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PREFACE 


This Final Report discusses the investigative results of a silicon-on-ceramic 
process for producing silicon sheet material. It was performed under the 
JPL Large-Area Silicon Sheet Task of the Low-Cost Solar Array Project. 

The objective of the program was to investigate the technical and economic 
feasibility of producing solar-cell sheet silicon by coating one surface of 
carbonized ceramic substrates with a thin layer of large-grain polycrystal- 
line silicon from the mejlt. The effort was divided into several areas of 
investigation in „rder to most efficiently meet the goals of the program. 
These areas include: 1) dip-coating; 2) continuous coating; 3) material 
characterization; 4) cell fabrication; and 5) theoretical analysis. The 
work in the continuous-coating technology development area was funded di- 
rectly by Honeywell Inc. from 1 January 1979 to 30 September 1979. Additional 
support to certain phases of the continuous-coating effort was contributed 
by the Solar Energy Research Institute (SERI) under the title of "Supported 
Growth of Polycrys taliine Silicon Sheet on Low-Cost Ceramic, Carbon, or 
Reusable Substrates," SERI contract No. XS-9-8119-1. Results from each of 
these separately funded efforts are included in this report for completeness. 

The work was performed during the period 21 October 1975 to 31 December 1980 
at the Honeywell Corporate Technology Center in Bloomington, Minnesota, 
under JPL contract No. 954356. 

During the course of the program, the Project Managers were P.W. Chapman 
and J.D. Zook, and the Principal Investigators were J.D. Zook, J.D. Heaps, 
and B.L. Gi-ung. Others who participated in this program are listed in another 
section of this report. 
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SUMMARY 


The objective of this Research and Development program was to investigate 
the technical feasibility of producing solar-cell-quality sheet silicon to 
meet the Department of Energy (DOE) 1988 overall price goal of $0. 70/watt. 
With the si 1 icon-on-ceramic (SOC) approach, a low-cost ceramic substrate 
is coated with large-grain polycrystalline silicon by unidirectional 
soldif ication of molten silicon. This R&D effort was divided into several 
areas of investigation in order to most efficiently meet the goals of 
the program. These areas Include: 1) dip-coating; 2) continuous ccfing- 
designated SCIM-coat ing, an acronym for Silicon Coating by an Inverted 
Meniscus; 3) material characterization; 4) cell fabrication and evaluation; 
and 5) theoretical analysis. 

Both coating approacbea were successful in producing thin layers of large- 
grain, tolar-cell-quality silicon. The dip-coating approach was initially 
investigated and considerable effort was given to this technique. The 
SCIM technique was adopted later because of its scale-up potentiality and 
its capability to more conveniently produce large areas of SOC. 

Large (10 cm x 100 cm) ceramic substrates were developed by Coors Porcelain 
Co. for the purpose of this investigation. Tie substrates contained slot 
openings for making electrical contact to the bottom surface of the silicon 
coating. Mullite was the ceramic selected. It was composit ionally modified 
to be thermally compatible with silicon. 

By virtue of unidirectional solidification, the thin (50 to 250um) silicon 
layers contain columnar grains as large as 5rara in width and several centi- 
meters long. The single grains are heavily twinned with their boundaries, 
as well as the normal boundaries, generally normal to the surface of the 
layer. Thus, the layers are, in general, only one crystal thick. The 
silicon-to-ceramic bond is mechanical in nature and has proven to be re- 
liably secure in practice. 

It was determined that substrate slot openings which are oriented perpendi- 
cular to the growth direction perturb the silicon structure to a lesser ex- 
tent than those which are parallel to the direction of growth. In either 
case, coverage with silicon is satisfactory. The thickness of the layer is 
inversely related to the melt temperature and is also proportional to the 
inverse of the coating speed squared, when the heat of fusion is removed 
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from the growing layer by radiation alone. It was experimentally demon- 
strated, however, that the growth rate can be substantially enhanced by the 
addition of convective cooling, and coating rates commensurate with project's 
economic goals can be achieved in this manner. 

The mullite substrates are slightly soluble in molten silicon, causing the 
melt to become progressively contaminated using the dip-coating process. 

This is minimized with the SCIU technique oy reducing the contact time and 
area the substrate has with the melt. 

The surfaces of SOC coatings are smooth and shiny and need no further pre- 
paration other than a standard clearup procedure prior to forming the p-n 
junct ion . 

SOC solar cells can be <;producibly produced with an average conversion 
efficiency of 9.6% (AMl, AR). The slot openings in the ceramic substrates, 
however, must cover approximately 50% of the total substrate area to 
achieve this level of cell performance. The need for this large slot open- 
ing ratio is to minimize internal series resistance in the cell, a fact 
that resulted from a supporting analysis of parasitic resistance. The high- 
est conversion efficiency for a dip-coated SOC cell was 10.54% (AMI, AH) and 
the highest efficiency for a SCIM-coated SOC cell was 7.6% (AMI, AR). The 
higher efficiencies are produced with a PH 3 diffusion at 850°C, followed 
by a slow cooldown at about 5°C/min. 

Theoretical modeling of SOC solar cells shows that present cell performance 
is limited by diffusion length, and that 13% SOC cells are possible with 
a diffusion Length of 50pm. 

The average diffusion length, L, in SOC material is 15 to 25pm. The dif- 
fusion i- r.gth within good grains is three times the average diffusion 
length. Closely spaced grain boundaries, high dislocation densities, and 
subsurface grain boundaries have been identified as structural causes of 
the reduction in L. 

2 2 

Cost estimates of the SOC process are $13/m added value, or $17. 25/m in- 
cluding silicon at $14/kg. Based on 11% module efficiency, this corresponds 
to a sheet price of $0.125/W added value, and $0.166/W , including silicon. 

r r 

The analysis shows that the cost of ceramic is the largest cost driver in 
the SOC process. 
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INTRODUCTION 


This RStD program began 21 October 1975. Its purpose was to investigate the 
technical and economic feasibility of producing solar-cell-quality shoot 
silicon by coating Inexpensive coramic substrates with a thin layer of 
polycrystalline silicon. The coating methods developed were directed toward 
a minimum-cost process for producing solar colls with a terrestrial con- 
version efficiency of 11% or greater. At this period in time, sheet 
silicon, for the most part, was produced by growing cylindrical single- 
crystal ingots and slicing them into round wafers. This slicing effort not 
only requires labor and time, but is accomplished by a costly diamond-edged 
saw which in turn produces u kerf loss of expensive silicon. The sliced 
wafers are then mechanically polished and subsequently chemically etched to 
remove surface damage. For the purpose of safe handling, those wafers are 
made thicker than is necessary to produce solar cells of a reasonable con- 
version efficiency. Once this type of solar cell is fabricated, much less 
than half of the starting silicon material plays a role in the photo- 
voltaic conversion of solar energy. The numerous processes required are 
labor-intensive and the round geometry of the final cell does not contribute 
to an area-efficient solar module. 

The s i 1 icon-on-ceramic approach, on the other hand, eliminates man 'f the 
above costly features. With the SOC process, a graphite couting is applied 
to one surface of an inexpensive ceramic substrate, causing molten 
silicon to wet, upon contact, only that suface of the substrate to produce 
a uniform, thin layer of large-grain polycrystal 1 iue silicon. With this 
process, nearly all of the silicon consumed is used in energy conversion. 

The surfaces of these thin layers are smooth and shiny and need no further 
preparation othor than a cleanup before being processed into solar cells. 
Such layers are coated onto large areas of rectangular substrates, thereby 
maximizing the active-area/total area ratio of a final solar module. 

During the course of this program, two silicon coating techniques were 
invest igated : 

• Dip-coating 

• Continuous SCIM-coatlng (SCIM is an acronym for Silicon Coating by In- 
verted Meniscus) 




Figure L illustrates tho principle of dip-coat Inn. A carbon or graphite 
coating Is applied to one face of a ceramic substrate. Tho substrate is then 
dipped Into and pulled from a crucible of molten silicon. As the substrate 
is withdrawn from tho melt, unidirectional solidification occurs on tho 
carbon-coated face of the substrate, forming a thin layer having large 
columnar silicon grains. 


CARBON "COATED 



Figure 1. Silicon dip-coating principle. 

The continuous SC1M technique is Illustrated in Fig. 2. Molten silicon is 
contained in a fusod-sillcu crucible having an attached trough as shown. 

A plunger is lowered into the crucible, displacing molten silicon along the 
trough to form a raised meniscus. As the ceramic substrate passes over this 
trough, tho meniscus contacts the carbon-coated side of tne substrate, 
forming a thin silicon layer which is virtually identical to that of the 
dip-coating technique. Since only the carbon-coated side of tho ceramic 
comes into contact with the melt, there is, in principle, no direct contact 
between the ceramic and the melt. Thus, contamination of the melt is greatly 
reduced as compared with that obtained from the dip-coating technique. 
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TECHNICAL DISCUSSION 


SHEET SILICON COATING FACILITIES 


Introduction 


The silicon dip-coating technique, discussed in the introductory section of 
this report, initially provided the basis for this R&D program. 

To satisfy the earlier needs of tho program, two dip-coating facilities were 
designed, constructed, and successfully used in pursuit of the work. As the 
program progressed, other sheet silicon coating techniques were considered 
and the SCIM technique (also previously discussed) was adopted and success- 
fully developed. In pursuit of this work, two SCIM-coating facilities were 
also designed aud constructed. 

This portion of the report addresses the purpose and function along with a 
detailed description of these four coating facilities. 

Dip-Coater No. 1 

To perform systematically an investigation of growth parameters which affect 
grain size, dislocation density, and layer purity, a versatile and reliable 
dip-coating facility was assembled. To the greatest extent possible, the 
facility was designed and constructed in such a way that its characteristics 
did not further complicate the problems to be solved. Since modifications 
were expected from time to time, sufficient design flexibility permitted 
such modifications to be performed with minimum time and expense. Once the 
necessary modifications were made, however, which permitted this coater to 
routinely produce good reproducible sheet silicon coatings, its design was 
for the most part frozen and it was then used for providing sheet silicon 
for the solar cell development effort. 

A sectional drawing of the facility is shown in Fig. 3 and a photograph of 
the assembled facility is shown in Fig. 4. 
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(1) Substrate Thermocouple Lead Wires 
' 2 ) Dip-Pull Quartz Tube 
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(15) Melt Thermocouple Positioning Rod 

(16) NichromeWire 

(17) Vacuum Gate Valve 

(18) Melt Thermocouple Weight 

(19) Water-Cooled "After Heater" Base Plate 

(20) Graphite Heat Shields ("After Heater") 

(21) Carbon Felt Insulation ("After Heater") 

(22) Graphite "After Heater" 

(23) Crucible Holder (Graphite) 

(24) Graphite Silicon Melt Heater 

(25) Carbon Felt Insulation (Silicon Melt Heater) 

(26) Graphite Heat Shield (Silicon Melt Heater) 

(27) Water-Cooled Silicon Melt Heater Base Plate 

(28) Quartz-Covered Melt Thermocouple 

(29) Quartz Crucible 

(30) Dip-Coating Facility Main Frame 

(31) Silicon Melt Heater Control Thermocouple 

(32) Crucible Positioning Rod 

(33) Water-Cooled Heater Electrical Leads 

(34) Melt-Dipping Chamber Pumping Port 

(35) Ceramic Insulation 

(36) Water-Cooled Surfaces 

(b) Nomenclature. 

Figure 3. Dip-coater No. 1 (concluded). 
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Figure 4. Dip-coater No. 1 and associated equipment. 


The coater included the following features: 

• 4 1 melt crucible (item No. 29) capable of dip-coating 6.5cm x 5.0cra 

substrates. This capability permitted us to produce sheet silicon 

2 

coatings in excess of the 20cm JPL requirement. 

• A purified graphite silicon melt heating element (item No. 24) powered 
by a temperature-controlled 25kW power supply. The heater assembly is 
mounted on a water-cooled copper base (item No. 27) and has only puri- 
fied graphite cylinders item No. 26 and carbon felt (item No. 25) for 
insulation. This reducf s melt chamber chemical contamination to a 
minimum. 
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• A Honeywell three-mode, precision temperature controller. 

• A water-cooled, stainless-steel melt chamber (item No. 36) and viewing 
port assembly. This water cooling is needed to offset the lack of 
heater-element insulation. This feature also contributes to a low 
level of chemical contamination. The cleanable viewing port is used 
for observing the substrate as it is entering and leaving the melt. 

• A gas-lock chamber (item No. 4) which permits the operator to load un- 
coated substrates and remove coated ones without shutting down the 
melt power. This feature is made possible by a large vacuum gate valve 
(item No. 5) which is located between the dip-coating chamber and the 
gas- lock chamber. The gas- lock chamber is readily purged with argon 
prior to opening the gate valve. 

• A programmable, variable-speed dipping and pulling mechanism. 

• A two-pen Honeywell electronic-type chart recorder for recording sub- 
strate position and temperature, both with reference to the melt. 

To efficiently characterize the optimum grov/th parameters it was 
mandatory that the growth conditions be accurately recorded for each 
sample substrate that was dip-coated. By doing so, as the physical 
properties of each coated sample were analyzed with respect to their 
particular growth conditions, positive direction could be given to 
the investigation. To provide this capability, the dipping and pulling 
mechanism was designed to provide an electrical readout which, when 
combined with the recorder's chart speed, provided dip-pull rate in- 
format ion . 

• An afterheater (item No. 22), if needed, located immediately above the 
silicon melt. This heater, however, was never used for layers grown 

at normal (3 to 6 cm/min) pull rates, but was used later in the 
program for fast (> 12 cm/min) growth studies. 

• A capability for orienting tie substrate at an angle while performing 
the dip-pull operation. There were no well-defined theories which 
suggested that angle dipping may enhance grain size, but intuitively 
we expected that dipping an angle-oriented substrate would modify 
the liquid-solid interface, perhaps in a beneficial way. Future ex- 
periments, however, demonstrated that this was not the case. 
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• Two adjustable gas flowmeters, one for regulating and monitoring the 
gas flow through the melt chamber and another for purging the gas-lock 
assembly . 

• A two-way gas valve for transferring from the argon gas cylinder to 
another cylinder should another type gas be desirable during a 
particular run. 

• A water failure shutdown switch. 

• A control panel for operating the facility. 

Dip Coater No, 2 

As the program progressed and dip-coater No. 1 was producing excellent 
large-grain sheet silicon coatings, we became reluctant to further modify 
No. 1 for experimental purposes. Dip-coater No. 2 was then designed and 
constructed to perform future experimental functions. The major task for 
this coater was to define the optimum conditions necessary to grow reason- 
ably thick (i.e., lOOum) films at high rates (i.e., 9 to 18 cra/rain). The 
design of this experimental dip-coater closely resembled that of dip-coater 
No. 1 with a few major exceptions. Dip-coater No. 2 was larger to allow 
more access around the liquid-solid interface and to allow special after- 
heaters to be used. The initial attempts to improve the growth rate while 
maintaining a film of reasonable thickness would involve the use of cold 
fingers to increase the temperature gradient at the liquid-solid interface. 
Dip-coater No. 2 had more access ports to accommodate the extra heating 
and/or cooling equipment necessary to grow at higher rates. The main ex- 
perimental approach used to increase growth rates is to remove the latent 
heat of solidification by forced convective and radiative cooling through 
the use of cooling shoes placed near the liquid-solid interface. As will be 
discussed below, when this is done, an afterheater must also be used to 
prevent fracture of the coated substrate during cooling. 

A photograph of dip-coater No. 2 is shown in Fig. 5. There are four access 
ports in the central portion of the chamber. Two are aligned at 30° to the 
geometrical center of the parting plane between the upper and lower 
chambers. There are four access ports in the top of the chamber that also 
point to the center of the parting plane and are inclined at an angle of 
fc>0°. During a run, the melt surface is typically level with the parting 
plane . 
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The heater assembly is shown in Fig. 6. A silicon charge Is shown inside a 
fused-sllica crucible in a graphite crucible holder. The crucible holder in 
turn is situated inside a graphite resistance heater which is surrounded by 
a heat shield consisting of two concentric graphite cylinders spaced with 
graphite wool. The heat shield, heater, and crucible holder can be inde- 
pendently positioned. The glass tube on the right of the heater is used to 
blow argon gas over the melt surface during some runs to minimize SiO 
formation on the silicon coating when cooling shoes are not use The heater 


is a picket-fence design considered to 
neater used in the other coater. 



Figure 5. 

Experimental dip-coator. 


be more rugged than the thin-wall 



Figure 6. 

Heater assembly in experimental 
dip-coater . 


The position of front and oack cooling shoes with respect to the substrate 
and heater assembly is shown in Fig. 7. The cooling shoes are machined from 
nickel and contain passages for both water cooling and gas flow. The cooling 
ard gas flow in a cooling shoe is shown in Fig. 8. The gas exits through 
holes drilled in the end of the shoe. This was not an acceptable situation 
because the nonuniform cooling on the silicon results in coatings with ve - *- 
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tlcal strlatlons. Examples will be shown later. On the other hand, nonunl- 
forra gas coverage clearly shows the effects of gas cooling on the growth 
pattern during dip-coating. Future modifications of the cooling shoes pro- 
vided more even gas coverage over the surface of the substrate. 


When forced convective cooling Is provided at the liquid-solid Interface 
during dip-coating, the thermal shock Is almost always sufficient to frac- 
ture the substrate particularly at higher pull rates. To minimize thermal 
stresses In the coated substrate, an afterheater must be used. Early 
attempts with "passive'' afterheaters (l.e., those heated only by conduction 
or radiation from the melt zone) were not satisfactory and an independently 
powered graphite resistance afterheater was eventually Installed. 



Figure 7. 

Photograph showing front and back 
cooling shoes in place in experi- 
mental dip-coater. 


GAS 


t I t f 



and gas flow in cooling 
shoes . 


The heater is a thin-wall, cylindrical, split-resistance element similar to 
those used early in the program as crucible heaters. The heater is powered 
by unregulated a-c current through two transformers in series. The relation- 
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ship of th« crucible, cooling, and afterheater la shown in Fig. 9. The cool- 
ing shoe* are inclined at an angle of 20° to the melt surface. The arrows 
on the sketch Indicate the direction of gas flow In each shoe. Note that In 
the shoe facing the uncoated r.lde of the substrate, the holes are drilled 
parallel to the axis of the i;hoe, while the holes In the other shoe are 
drilled perpendicular to the substrati surface. This arrangement allows the 
back side of the substrate to be cooled slightly below the liquid-solid In- 
terface, while the front side Is cooled above the interface. Typically, the 
gas flow rate on the back side of the substrate Is four times that through 
the front shoe. 


AFTER HEATER 



Figure 9. Schematic showing relationship of 
afterheater, cooling shoes, and 
crucible assembly in experimental 
dip-coater . 





SCIM-Coater No. 1 (SCIM-I) 


It was recognized early in the program that even though the dip-coating 
technique was successfully producing thin, smooth sheet silicon coatings of 
solar-cell quality, the method, nevertheless, was not easily adaptable to 
a continuous coating technique. In aldition t this, experimentation had 
shown that the silicon melt was being progressively c< 'dominated with each 
successive substrate that was silicon-coated. In an effort to overcome 
these disadvantages, the SCIM technique (described in the Introduction cf 
this report) was adopted as a promising cont inuous-cout ing method. The 
principal advantages of the SCIM technique are: 

• Its coating principle vastly reduces the area of, and the time that 
the substrate is in contact with, the silicon melt. This feature 
lessens the possibility of the melt being progressively contaminated 
by the substrate. 

• This technique permits the use of much larger substrates, which should 
allow the individual c. ystalline grains to develop larger in size. 

• The scale-up coating throughput of this method is potentially superior 
because it eliminates the dippitg and soaking time requirea with the 
dip-coating technique. 

• The thermal stability of the solidification .tone of a continuous 
coater should also be superior to a dip-coater, since the melt temper- 
ature would not be constantly upset by the immersion of unheated 
substrates . 

Figure 10 shows photographs of the SCIM-I facility. The coater is divided 
into five sections. (1) a compartment for substrate loading; (2) a preheater 
section for bringing the subs ate up to temperature; (3) a silicon coating 
compartment; (4) a postheater section for cooling the substrate; and (5) 
a compartment for substrate removal. During a coating run, these compart- 
ments are flushed ith an inert gas, most generally argon. The coater is 
powered by four temperature-controlled power supplies, two of which power 
the upper and lower preheaters. The third supply powers the main melt fur- 
nace and *he fourth supply the coating trough furnace. 
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(b) SCIM center No. 1 showing bus bars from 
remotely located power supplies. 


Figure 10. Photographs of SC’IM-l. 


( 
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Three of those four heaters are shown In FIr. 11. Tho preheater assembly 
cons i sled of u long Rruphlte tunnel throuRh which tho substrates pass. A 
Ioiir 1 1 udlnn 1 tomporaturo Rrudlont was created aloiiR tho tunnel by huuliuR 
om> oiul of tho tunnel with the upper and lower preheater elements. This 
tunnel assembly was subsequently surrounded by a dual-walled insulated 
housttiR. also shown til FIr, 11. A cross-sucl tonal view of tho preheater as- 
sembly is illustrated in FIr. 12. The postheater Is similar to the pre- 
heater except that it is sllRlitly loiiRor and passively heated by t lie 
coatiiiR t i'ourIi and preheater. FIrui'oh ll and 12 show the loitRitudinul ther- 
mal profiles which existed with this heat ln« conf tRurat ion . Initially, 

SIMM- 1 ’ s transport mechanism consisted of stain Loss-steel chain conveyors 
which carried the substrates to, throuRh, and from the irrowth cl, amber. The 
chain conveyors were driven by a "D i r 1 lok"-con t ro t led , variable-speed drive 
system which was adjustable for Rfowtlt velocities of tenths of em/miii to 
many cm/mtn. Given a relatively flat substrate, the fore and after chain 
conveyors were spaced to allow the substrate’s center of Rfavily to carry it 
from one conveyor to t lie other. Special silicone rubber end Rates permitted 
substrates to pass throuRh the coaler with a minimum of run leukuRO. While 

the facility was primarily dost Rued to hor t /.out a 1 1 y coat substrates, pro- 

vision was made for tiltiiiR t ho entire facility throuRh about * 20°. This 
proved to be a valuable feature, since it was never possible to successfully 
c,.at substrates in tho hori/.ontul position. 

The isotherms shown in FIrs. 14 and 1T> Illustrate that tin* transverse 
temperature uniformity near t he coatiiiR t roiiRh was very bad. Likewise, with 
this furnace con f i Rurut i on , the lotiRiludlnal temperature profile shown 
in FIr. tl did not indicate a sharp enouRh temperature drop in the solidifi- 
cation /.one of the confer. To correct both of these thermal inadequacies, 

a new type of t rouRh-heat .or assembly was installed in the center. 

This heater was fabricated in such a way that the Rraphtte l i-ourI) holder 
itself was an active resistive element ot the assembly ( see FIr. 111). The 
rectaiiRular Iors of this heat Iiir assembly served the final substrate huatiUR 
requirement. The power density of each Ior was Rrnduutod in such a way that 
the substrat«* tomporuturc never exceeded the melt (hr point of silicon. 
Contrary to the oriRinul approach (see FIr. 11). the heat Iiir /.one now ter- 
minated abruptly at the downstream edRO of the t rotiRti . 


1 H 


TEMPERATURE ('C, 



Figure 11. SCIM-1 preheater longitudinal thermal profiles. 
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Ft nu re 12. Cross-sect tonal view looking into 
SCIM-I preheater compartment as 
viewed from the coating chamber. 
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Figure 14 . 

SCIM-I preheater isotherms 
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Figure 15. SCIM-I postheater isotherms. 




Figure 16. 

New trough/substrate heater 
(also showing quartz crucible/ 
trough piece). 
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As shown in Fig. 17, the upper substrate heater was shortened, thereby 
shifting the peak temperature in an upstream direction. 


UPPf* 

<- KtAtEP 



Figure 17. SCIM-I coating chamber configuration 
showing combination trough/substrate 
heater . 


Following the above modifications, the liquid-solid interface, if desired, 
could now be shifted until solidification occurred at the edge of the 
trough. This is the dip-coating equivalent of having the surface of the 
melt freeze around the substrate as it is being withdrawn from the melt. 

The temperature of this shorter heater was now a critical coating parameter. 
If too hot, solidification, as usual, occurred too far downstream. If not 
hot enough, the substrate would freeze to the melt. In spite of this fact, 
the solidification position could now be controlled. 

The original chain-conveyor transport mechanism proved to be inadequate. 
Therefore, while the above furnace modifications were being implemented, a 
multiple-roller-type mechanism was installed in the coater. The lower 
roller was driven, whereas the upper roller served as an idler to maintain 
the alignment of the substrate. The rollers were made of graphite, since 
two sets of them operated at temperatues as high as 1200°C. The addition of 
these transport rollers required a modification to be made in the pre- 
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hoater and postbeater assemblies, also shown in Fig. 17. It was with this 
coating configuration that ths first succsssful continuous coating was 
performed. 

SCIM-Coater No. 2 (3CIM-II) 

In ths last ysar of this contract, a sscond continuous-coating system (de- 

slgnated SCIM-II) was dssignsd, constructsd, and mads opsrational. Its 

primary function was to dsmonstrats tbs capability of producing silicon 

2 

sheet at a throughput rate of 350 cm /min. 

To achieve this goal, the coater was dssignsd as shown in Fig. 18 to simul- 
taneously coat (side-by-side) two 12.5cm-wlde x 100cm- long mullite sub- 
strates The coating experience acquired using SCIM-I was a valuable asset 
in the design of SCIM-II. The design of this new coater is shown schemati- 
cally in Fig. 19. 



Figure 18. SCIM-II exterior showing two side-by-side 
substrates entering the coater. 


There are several key features which are different from the SCIM-I system. 
They are: 

1) All assemblies within the coater are attached to a common water-cooled 
chassis (see Fig. 20). This assures us that proper alignment of key 
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Figure 19. Schematic cross section of SC1M-II. 


parts will be maintained irrespective of operating temperature. 

2) All metallic parts within the coater are water-cooled to improve the 
purity of the silicon layer, and the majority of the remaining parts 
are fabricated from graphite which can be purified prior to assembly. 

3) The melt crucible and trough are contained in a graphite cradle as- 
sembly which permits convenient adjustment of trough-to-substrate 
distance during operation. 

4) Side heaters powered by a separate power supply were installed tc 
improve the transverse temperature uniformity in the substrate pre- 
heater. These side heaters, however, have proven so far not to be a 
necessary feature. 
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Figure 20. SCIM-I I interior showing water-cooled base 
and substrates emerging from preheater. 


5) The substrate transport mechanism consists of aluminum rollers which 
are fitted with several silicone rubber "0"-rings. The upper, as 
well as the lower roller, is driven by a loosely coupled spur gear. 

The entrance and exit rollers are positioned external to the cover 
containing the inert atmosphere needed for growth. In operation, they 
have proven to be very effective in performing the transport function. 
With this method of substrate transport, however, the substrate used 
must be at least as long as the distance (97cm) between the rollers. 
This can be a problem if a substrate breaks during coating. In this 
respect, the multiple-roller method used in SCIM-I is superior. 
However, the numerous rollers in SCIM-I, some of which operate at 
temperatures as high as 1200°C, proved to be only marginally 
dependable. 
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6) The water-cooled cover on SCIM-II serves only the function of con- 
taining the inert gas. As such, when this cover is removed, the 
inner walls are easy to clean. Likewise, when removed, the entire 
inner-workings of the coater are easily accessible for maintenance. 

The preheater and postheater design used in SCIM-II was fashioned after 
that of SCIM-I. This design was used because these functions .'n SCIM-I 
performed well. The purpose of the preheater is to bring the substrate from 
room temperature to the melting point of silicon in a manner which neither 
fractures nor warps the substrate. This is done by using a narrow graphite 
tunnel through which the substrate passes. One end of the tunnel is 
heated by a graphite resistance heater, and thermal conduction along the 
tunnel creates the needed temperature gradient. The postheater is similarly 
constructed. It, however, is a passive element deriving its heat from the 
preheater and the coating trough. This coating configuration produces a 
temperature profile with a rapid fall in temperature immediately downstream 
from the coating trough, thus creating the required thermal environment 
for solidification. To ensure that the coated substrate will emerge from 
the machine at a temperature of less than 232°C (the maximum operating 
temperature of the silicone rubber "0"-rings), a water-cooled heat sink is 
attached to the downstream end of the postheater. The thermal modeling 
contributing to this design is discussed later in this report. In order to 
achieve a uniform transverse temperature, the preheater and postheater are 
substantially wider than the width of the side-by-side substrates being 
coated. 

The trough heater and the fused-silica trough-crucible assembly are also 
fashioned after the design used in SCIM-I. The quartz trough nestles in an 
electrically active trough heater. The cross section of this heater is 
reduced at the electrode end in order that additional power can be generated 
to offset losses by conduction to the water-cooled copper electrodes and 
by radiation through the viewing port. The trough heater is also made pur- 
posely long for temperature uniformity reasons. 

As with SCIM-T, SCIM-II can also be tilted to coating angles up to 20° to 
provide the necessary meniscus stability. To date, lOcm-wide substrates 
have been coated at angles of 10° and 20° these results are included else- 
where in this report. 

The most obvious problem during the initial operation of SCIM-II was the 
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buckling and/or breakage of the coated substrates. Experiments were per- 
formed without silicon in the system to clarify the situation. Initial 
experiments quickly showed that at high speed (15 cra/min) there was only 
minor buckling, but at lower speeds, the buckling became more severe, and 
below 6 cm/min the substrates broke. The results were qualitatively the 
same regardless of the silicon coating and regardless of whether one side 
of the substrate was carbon-coated. 


The original design of the entrance and exit tunnels was based on a design 
speed of 15 cm/min. The most interesting result is that at this high speed 
there has been little warping and no cracking of any substrates. This is 
true of the original design as well as all subsequent modifications. 

The problem is that without forced convective cooling or use of an asym- 
metric growth mode, there is no chance of getting the desired silicon 
thickness at high speed. The thickness is needed to demonstrate slot 
coverage and cell performance. Thus, we wanted to do initial coatings in 
the range of 4 to 6 cm/min. At these lower speeds, the substrates have 
warped, buckled, or cracked. The proposed explanation is that the longitu- 
dinal temperature profile must be the cause of the problem. At high speeds, 
the longitudinal profile is smoother than at low speeds because the heat 
capacity of the substrate causes a temperature lag. Thus, the purpose of 
the changes in tunnel design has been to linearize the temperature profile. 
Fortunately, this has reduced the problem at low speeds, and has not 
created any problem at high speeds. 


The explanation in terms of thermal stresses caused by longitudinal temper- 
atuie gradients is the same as that proposed by Surek for edge-defined, 
film-fed ribbon growth or EFG.* The origin of thermal stress and the 
strategy to avoid it are discussed in SOC Quarterly Report No. 15, pub- 
lished 31 July 1980. 


As discussed by Surek, 


J E b 


xx 



the equation for the stress is 

2 „ 


( 1 ) 


in the case where the second derivative of the temperature is significant. 
Here, a is the expansion coefficient, E is young's modulus, and b is one- 
half the substrate width. The solution is positive at the center (y - 0) and 


^A. D. Morrison, et al., "Large-Area Silicon Sheet by EFG," Mobile Tyco 
Energy Corporation Annual Progress Report, September 15, 1976. 


27 


negative at the edges (y • t b), Indicating tension and compression, respec- 
tively, when d^ T/dx^ is positive. At the ends of the substrate, the solu- 
tion is not valid, since o - 0. Thus, in the vicinity of the end, there 

are <j and o stresses which can also be derived analytically, 

yy xy 


It is obvious that the temperature profile cannot be entirely linear in 
any substrate which is heuted and then cooled and that both types of non- 
linearities must be encountered, as shown in Fig. 21. Thus, thermal stress 
must always exist. However, there is a strategy which, at least in prin- 
ciple, can reduce the residual thermal stress to zero. 1 This strategy is 
detailed in Fig. 21. The idea is to have the concave portions of the pro- 
file in regions with high yield stress. Thus, thermal stresses do exist, but 

they result only in elastic deformation so that when the sample is again 
2 2 

in a d T/dx ■ 0 region, the residual thermal stress is zero. 



Figure 21. Strategy to eliminate effects of thermal stress. 

In the high-temperature regions, the ceramic can 
deform without buckling or cracking, and stress 
in the silicon is annealed as it is generated. 

At low temperatures, only elastic deformation 
occurs. The intermediate temperature ranges are 
critical, and the temperature profile must be 
1 inear . 


At sufficiently high temperatures, the yield strength of both the ceramic 
and the silicon is low so that the sample can undergo plastic deformation 
with no buckling or cracking. Stress relaxation studies in silicon 1 show 
that at 1000°C the stress relaxation time is a few seconds Thus, in the 
vicinity of the trough and especially at the liquid-solid interface, the 
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presence of large non 1 Inear 1 1 ies in the temperature gradient will not have 
harmful effects. 

In summary, the strategy is to consider the problems in terms of three 
temperature zones. The following strategy applies to both the heating and 
cooling zones: 

a) The high-temperature zone (T > 1000°C) can tolerate large nonlineari- 
ties in the temperature profile. 

b) The intermediate-temperature zone ( 1000°C ^ T > ")0°C) must have a 
very linear profile. 

c) The low-temperature zone (500°C l T) can tolerate large nonlinearities 
such that the thermal stress does not exceed the elastic limit of 
silicon on ceramic. 

According to the design strategy based on stiess analysis (Fig. 21), a 
linear longitudinal temperature profile is sought for the temperature 
range of low yield strength, from 500 to 1000°C. The linear profile is 
necessary in both heating and cooling in order to eliminate or at least 
greatly reduce residual stresses in the coated ceramic panels. A second 
requirement is to provide rapid initial cooling of the coated substrate, as 
prescribed by heat balance at the silicon liquid-solid interface. The 
second requirement will not be in conflict with the first one if rapid 
cooling is restricted to the plastic deformation region, as indicated in 
Fig. 21. 

It is expected that the thermal design of a production coater can be opti- 
mized for a given fixed feed rate. In the experimental growth system, 
thermal design is complicated by the need to accommodate more than one 
feed rate; an ideal temperature profile is difficult to maintain over a 
wide rmge of feed rates. At low speed, the substrate is likely to heat up 
over-rapidly on first entering the heating chamber. At higher speeds, this 
is moderated by the thermal inertial time lag. On the other hand, high- 
speed panels tend to emerge from the cooling chamber insufficiently cooled. 
The first problem can be eliminated by heat-sinking at the entrance to the 
heating chamber, whereas the second problem is relieved by heat-sinking 
at the exit from the cooling chamber. A more flexible means of accommo- 
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dating different speeds is to vary the amount of insulation placed in the 
respective chambers. 

A great deal of ins'ght into these and other aspects of the tnerreal design 
problem has been provided by a mathematical analysis of the heating and 
cooling chambers. A brief thermal description of the mathematical models 
fol lows . 

The substrate panels are transported at ambient temperature to the heating 
chamber, which is shown in schematic cros3 section in the left half of 
Fig. 22. This is made up of pairs of active heaters, carbon guides (a'.so 
functioning as the passive heaters), and insulating layers. The lmm-thick 
ceramic substrates move left to right between the guides, ucquiring heat by 
radiation from the guides, so that they reach 1693K when they arrive at 
the right end of the tunnel. At this temperature, the melting point of 
silicon, they are ready for silicon coating in the chamber immediately to 
the right. It will be noted that the moving substrates extract a total of 
pCwtAT Joules of heat per centimeter of length, where p and C are the 
density and specific heat, w and t are panel width and thickness, and AT 
is the temperature rise. If the panel feed rate is v cm/sec, the active 
heaters then must supply energy at a rate greater than pCwtAT watts. The 
"greater than” may be made to approach "equal to" if wasted heat is mini- 
mized; that is, the quantity and placement of insulation are such as to 
prevent significant losses to the water-cooled case. 
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Figure 22. Schematic cross section of SCIM-II showing 
simplified arrangement of various thermal 
elements used in modeling, h * -esistive 
heater; c * carbon (guide); i * insulation; 
w * heat sink (when present); s * substrate. 
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After the silicon coating step, the substrates traverse a cooling chamber 
(right half of Fig. 22) of similar geometry to the entrance tunnel. The im- 
portant difference is the omission of the active heaters. However, the 
water-cooled case and variable insulation are essential. Inasmuch as cool- 
ing is passive with the case being the heat sink, the carbon guides are a 
hindrance to cooling though needed for mechanical support. The guides act 
like heat shields, interrupting the radiative transfer from silicon-coated 
panels to the case. The important heat escape path is vertical from sub- 
strate to guide to case, with a lesser contribution by way of longitudinal 
conduction through the guides. 

The need for design flexibility is reflected in the variable structure of 
the mathematical models; provisions are included also for experimental 
calibration : 

• The amount of insulation used in the heating and cooling chambers 
can be varied. Adjusting insulation may bo the quickest and easiest 
way to set up different substrate speeds. 

• The far ends of the cooling chamber guides have tended to reach ex- 
cessively high temperatures in both theory and practice. The theore- 
tical model includes an optional heat sink at these ends. 

• "Standby" versions of the models have been provided. These represent 
heating and cooling chamber with heaters active but no substrate 
p'esent. The main purpose for these versions is to calibrate the 
models against temperature measurements, which are normally taken 
with the chambers empty. 

Only longitudinal temperature profiles are analyzed here, since these are 
believed to be associated with warping and breakage. Transverse temperature 
variations al»o occur and are discussed elsewhere. The main problem caused 
by transverse variations is coating nonuniformity, and this is being suc- 
cessfully attacked by heater redesign. 

The physical relationship upon which the thermal models are based will next 
be summarized. All temperature profiles are assumed to be one-dimensional; 
thus, variations in the transverse direction E.re ignored. Several thermal 
elements interact with one another (Fig. 22) so that as many as five coupled 
heat equations must be solved. Basically, the electrical input produces 
heat in the resistive heaters (h), and this heat spreads along both heating 
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chamber guides (c) by thermal conduction. In the steady state, the heated 
guides release part of this energy by radiation to the entering substrute(s) 
and the rest is lost by radiation and conduction to the ambient. Thu model 
Includes a convective heat transfer coefficient, but this has been set 
equal to zero in the calculations reported here. A significant portion of 
the heat loss takes place from insulated surfaces due to the limited 
thicknesses of carbon felt. The substrate(s) acquires heat from the guides 
so that its temperature Increases from approximately ambient at the left 
to the coating temperature (1420°C). 

In the cooling chamber, the process is essentially reversed. The substrate 
is progressively cooled by radiation to the cooling chamber guides and by 
conduction in the downstream direction. Because of the heat-shield effect, 
the cooling rate is lower than it would be in the case of direct radiation 
t.o the 3Q0K ambient. The cooling problem is further aggravated at the far 
(right) end because thermal conduction within the guides tends to elevate 
the temperature at that end. This effect can be greatly relieved by heat- 
sinking (w). 

End effects are found to be significant and are included both at the exter- 
nal end facing the cool environment, and the internal ends facing the 
coating chamber. Not only the actual end surfaces are affected but also 
parts of the internal surfaces close to the ends The moving substrate is 
not completely coolud when It emerges from the cooling chamber. It is 
important that enough cooling has taken place at this point to avoid a 
further precipitous temperature drop outside of the chamber. In the numeri- 
cal calculations, a practical approximation is required for the boundary 
condition at infinity. We use T - T (room temperature) at a point some 
distance (a few can) beyond the tunnel end. Over this distance, the radia- 
tion environment is assumed to be at temperature T . It is easily verified 
that the approximate boundary condition does not affect the numerical solu- 
tions for the temperature within the chamber. Figure 23 shows calibration 
of the model against measurements taken by using a thermocouple instead of 
the substrate. 

The finite-difference grid for the heating chamber represents 1cm incre- 
ments in longitudinal position, or 41 points for a 40cm tunnel. Since there 
are five coupled layers (two heater- insulat ion layers, two guides, and one 
substrate), there are 5 x 41 - 205 nodul equations which are solved by a 
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Figure 23. 


Calculated versus experimental 
exit tunnel temperature with 
no substrate in chamber. 


direct- inversion method for sparse matrices. The matrix has a bandwidth of 
11 and is asymmetric due to the moving substrate. Since the radiation law 
is nonlinear and thermal conductivities are temperature-dependent, the 
equations are resoived iteratively with updated coefficients for the nodal 
equations. Convergence occurs in about 20 iterations. For the cooling 
chamber, a finer grid (0.33cm) is used, but heater- insulat ion layers arc 
not required. 

Heating Chamber (Entrance Tunnel) Results - Calculated temperature profiles 
are shown in Fig. 24 for two substrate speeds. The solid curve gives the 
substrate temperature and the dashed curve gives the mean guide temperature. 
(Iho two guides may differ up to 20°C from each other at any x position. ) 

The third curve consisting of dnshes and dots, represents mean guide temper- 
ature under standby condi t ions- the some electrical input but no substrate. 
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The high-speed substrate profile is nearly linear, indicating optimum de- 
sign for high speed. Electrical input is 135 W/cm of width, of which 92 
W/cm heats the substrate, the rest being absorbed by the ambient environ- 
ment. The standby profile is considerably elevated because the entire 135 
watts is wasted to the environment, thus requiring a higher average radia- 
tion temperature. 

The low-speed profile demonstrates some of the difficulties encountered 
in trying to optimize the basic chamber design for a range of substrate 
speeds. With all the insulation in place, the initial temperature rise in 
the substrate is extremely rapid. This is moderated by removing all of 
the insulation, giving the Indicated profile, which is still somewhat too 
steep and has excessive curvature. The curvature is initially convex but 
becomes concave in the middle due to the large radiation losses from the 
guides. It becomes convex once again at the end, where the active heaters 
are located. The electrical input, 198 W/cm, is almost totally wasted, with 
only 12 W/cm required to heat the substrate. As discussed in SOC Quarterly 
Report No. 13, published 15 February 1980, the theoretical profile for 
the low-speed case is best improved by increasing the thermal resistance of 
the guides - for example, reducing their thickness. A compromise solution 
is being investigated experimentally; carbon has been removed in transverse 
grooves at the ends of the guides. 

Cooling Chamber (Exit Tunnel) Results - The substrate temperature distri- 
bution in the cooling chamber should satisfy two requirements: (1) a large 
initial gradient to carry the latent heat of fusion out of the grown 
silicon layer; (2) a nearly constant gradient from 1000°C to 500°C . As 
previously explained, these requirements are difficult to satisfy simul- 
taneously under all conditions. The location of the insulation is critical; 
since the cooling curve has the greatest curvature at high temperatures, 
the insulation should be placed fairly close to the beginning of the 
chamber. On the other hand, the first requirement demands a fast initial 
cooling, and this dictates leaving a finite uninsulated gap at the very 
beginning. A 1cm gap was found to be sufficient to retain the full initial 
cooling rate. The next approximately 18cm of both guides are insulated for 
present experimental work, and the same insulation interval (1 to 19cm) was 
used in the numerical calculations plotted in Fig. 25. The effects of this 
insulation are readily observed by comparison with Fig. 26, which repre- 
sents the uninsulated case. 
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Figure 25. Calculated exit tunnel profiles at different 
speeds with insulation located as indicated 
(carbon guides are not heat-sunk). The solid 
profile corresponds to guides under operating 
conditions, and the dashed and dashed— and— 
dotted prof lies to guides in standby conditions 
different electrical inputs. 







Figure 26. Calculated exit tunnel profiles at different 
speeds without insulation. Carbon guides are 
not heat-sunk. The solid profile corresponds 
to substrate; the dashed profile corresponds 
to guides under operating conditions; and the 
dashed-and-dotted profile corresponds to 
guides in standby conditions. 
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Figure 27. Calculated exit tunnel profiles 
with heat sink at exit end. The 
solid profile corresponds to 
substrate; the dashed profile 
corresponds to guides under oper- 
ating conditions; and the dashed- 
and-dotted profile corresponds 
to guides in standby conditions. 
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As discussed above, the high exit temperatures are intolerable. A water- 
cooled heat sink was therefore added at each guide end, and, according to 
calculations (Fig. 27), a substantial lowering of the exit temperature was 
expected. As discussed below, this was indeed the result. 

To eliminate this substrate warpage and breakage in practice, a large 
number of modifications were made in the system: 

• Entrance tunnel heaters were simplified and moved forward. 

• The side heaters were found to be unnecessary with the new entrance 
tunnel heaters. 

• The method of measuring temperature profiles was improved by relo- 
cating the thermocouples. 

• Insulation was added or removed as required to obtain the desired 
profile . 

• The trough heater was designed to compensate for unavoidable con- 
duction and radiation heat loss nonuniformities. 

• The exit tunnel was connected to a water-cooled heat sink as indicated 
by modeling and by measurements. 

• Closed-loop control for the temperature in the trough was implemented. 

• The angle of the quartz trough was changed so that it tilts backwards 
in operation. 

These modifications gave much improved operation. 


Temperature Profile Measurement - Careful measurement shows that when the 
thermocouple is inserted from the entrance port, the reading is different 
from when the same thermocouple is at the same position when inserted from 
the exit port of the coater. The discrepancy is opposite to that expected 
on the basis of conduction of heat along the thermocouple. For example, 
the temperature above the trough is 30 degrees hotter when measured from 
the exit port, where heat is conducted away from the probe. The only rea- 
sonable explanation for the discrepancy is that the thermocouple is re- 
sponding to the direct radiative environment and that the presence of the 
ceramic probe blocks the radiation. Thus, in the above example, the probe 
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pointed toward the hot entrance tunnel reads hotter than whon it is 
pointed at the cold exit tunnel. A different method of reading the temper- 
ature profiles was adopted. A series of thermocouples was installed by 
drilling holes into the half-inch plates and cementing a thermocouple into 
each hole. The results of the measurements which followed are discussed 
below. 

The temperature of the exit-tunnel plates as indicated by the type-K thermo- 
couples shows a sharp discontinuity at the edge of the insulation. This 
result also follows from the thermal modeling (see p ig. 27). Figure 28 
plots this profile and sketches the insulation position. The function of 
these carbon plates is to take the sample from a high temperature, where the 



Figure 28. Temperature profile of partially insulated 
exit tunnel. The horizontal scale is mea- 
sured in cm from the center of the trough. 
The data were taken with thermocouples im- 
bedded in the carbon plates to avoid radi- 
active coupling problems. 


ceramic is soft, smoothly down to a low temperature where the sample is 
strong enough to withstand the strains (Fig. 21). There is a change in slope 
at the end of the insulation. In this case, breakage of substrates oc- 
curred presumably due to the sizeable second derivative in the profile. 
Figure 29 shows the case where there was no insulation at all, and the pro- 
file is quite linear. Obviously, no insualtion is not very energy efficient. 
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Figure 29. Measured profile of uninsulated exit tunnel. 

It also prevents bringing the leading edge of the exit tunnel up to the 
1000°C temperature range required by the mullite. We then covered the exit 
tunnel with a uniform layer of graphite felt from the leading edge of the 
carbon all the way back to overlap the water-cooled copper plates. Figure 
30 shows the profile of the present exit tunnel and shows the position of 



Figure 30. Measured profile of exit tunnel with 
insulation along full length. 
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the insulation and water-cooling. A similar set of thermocouples was 
mounted in the entrance tunnel to provide a profile of the entire SCIM-II 
coater. Figure 31 is a plot of the profile of the entrance and exit cunnuls 
at two different input power levels. Considerable crosstalk between the 
trough heater temperature and the entrance tunnel heaters is indicated by 
the fact that the trough heater alone cannot reach 1400 U C. Also note the 
900°C temperature at the leading edge of the exit tunnel corresponding to 
400 A. Under these conditions, thermal stresses were excessive. We have 
concluded that at 900°C the exit tunnel is too low to begin cooldown on the 
substrate, whereas at 1000°C or more the 10cm by 100cm samples do not 
buckle at low or high speeds. 



Figure 31. Measured profile at two heater powers. 


The data in Figs. 28 through 31 all have been taken without samples or 
silicon in the trough. Figure 32 is the operating profile of SCIM-II, in 
eluding all of the following cases: 

• Substrates going through from 3 to 10 cm/min 

• Mo substrates at all 
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Silicon in the trough 


• No silicon in the trough 

• Coating SOC at 5 cm/min 

The height of the bars shows the total variation in temperature under these 
conditions. This variation was only t20°C at the most critical point (at 
the leading edge of the exit tunnel). With these conditions, the system can 
coat silicon on 10cm x 100cm ceramic substrates without buckling or crack- 
ing due to thermal stress. 



Figure 32. Measured profile during a coating run. 


Heater Insulation - Resistive carbon heaters are designed and built to sur- 
round the crucible or trough as much as possible. For the entrance tunnel, 
the heaters are adjacent to the carbon plate to be heated. In these cases, 
only one side of the heater sees the material to be heated, whereas the 
other side (outside, backside) radiates into some free-standing insulated 
structure which should act as a radiation shield. As the element to be 
heated approaches operating temperature, e.g., 1400°C in our system, it no 
longer requires a great deal of power from the heater, whereas the box 
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between the heater and the insulation that ie holding the insulation con- 
tinues to drain power from the system. 

Throe of our four heaters have been modified to minimize this loss and 
focus the radiation in the direction intended, by simply attaching the in- 
sulation directly to the electrically active heater element. The structure 
to hold the insulation is now the heater itself and no longer need be the 
elaborate structure previously used. The hottest part of the heater is now 
in contact with the insulation with reduced radiation, convection, or con- 
duction losses on that side. The bare portion of the heater still radiates 
directly into the intended heat sink. 

The insulation on the top entrance tunnel heater is simply a straight flat 
blanket cemented on top of the heater. The bottom heater is identical ex- 
cept for the insulation being underneath. The trougn heater is somewhat 
more elaborate in that three sides of the heater are insulated at the end 
where the current connections are made, two sides are insulated where the 
heater reaches into the crucible container and nearly touches the crucible 
heater. There is a small stretch of three-side insulation between the cru- 
cible box and the coating area as well. All along the trough heater the in- 
sulation is divided, since there is a voltage drop between the .wo sides. 
The 2mm gap in the heater is sufficient to maintain separation of both 
heater and insulation. Figure 33 gives a cross-sectional view of heaters, 
plates, and insulation. 
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PH INSULATION E3 PASSIVE PLATES 


Fig ? 33. Schematic cross section of 
coating zone, including 
heaters, insulation, and 
passive plates. 
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Coating results using this thermal configuration are reported in the Sheet 
Silicon Growth section of thin report. 

CERAMIC SUBSTRATE DEVELOPMENT EFFORT 

Introduction 


Early in 1975, when silicon dip-coating was first conceived, carbon appeared 
to be a logical choice as a substrate material. Small substrates were 
machined from graphite and dip-coated with silicon. Thin, smooth silicon 
coatings were produced having columnar grains similar in structure to the 
most recent SOC layers. Samples of Grafoil (a graphite foil) were also 
dip-coated at that time with surprisingly good results. A promising char- 
acteristic of carbon or graphite was the fact that it is electrically 
conductive and could conceivably form the back contact to the silicon 
layer. A preliminary investigation, however, into the economics of carbon 
substrates proved to be unfavorable and our attention turned to ceramic 
as a potentially low-cost substrate material. Sections of mullite tubing 
were then cut up into small substrates and successfully dip-coated with 
silicon. In the early phase of this program, Honeywell's Ceramic Center was 
asked to assist us in selecting a suitable ceramic. 

Later in the program, we learned of Coors Porcelain Company's expertise in 
large-area substrates and they were consulted and subsequently very much 
impacted the SOC effort. The results of these efforts are discussed below. 

Ceramic Material Selection 


The Honeywell Ceramics Center played an important role in helping us in the 
selection of a suitable ceramic material. This selection had to be made on 
Che basis of thermal expansion coefficient, thermal shock resistance, 
chemical purity and, most importantly, cost. 

The Ceramics Center, in addition to purchasing or developing substrate 

materials, spent time examining several processing approaches, which could 

be scaled-up to produce square-raeter-size substrates. These approaches 

included rolling, casting, doctor-blading, and hot-pressing. Rolling of a 

high-plasticity, mullite material was examined and substrates as large as 
2 

500cm were formed and fired. Rolling and doctor-blading of a calcium- 
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aluruinate bound alumina material were shown to be feasible. Casting and 

doctor-b' ading of silica substrates were demonstrated. Hot-pressing of 

2 

fiberglass products into large, 500cm sheets was also accomplished. 

Standard 5c.n x 7cm test samples for this program were prepared or received 
from vendors as indicated in Table 1. Each of these materials is discussed 
in the paragraphs that follow. 

Table 1. Substrate materials investigated. 


Material 

T rade Name 

Source 

Fabrication 

Process 

Firing 

1 emperature 

Mullite 

MV-20 

McDanel Corp. 

Rolling 

1640°C 20% glass 

MulUte 

MV-30 

McDanel Corp. 

Caettng 

15% glass 

Mullite 

Flberfrax 

Caroorundum 

Mot pressing 

MOOT 

Mullite 


American Lava 

Pressed 


t'ordiert’e 

Alaimag 701 

American l.avn 

(purchased) 

-- 

( ordiertte 

A-3171 

Du-Co 

(samples) 

-- 

Alumina 

Alolmag 614 

American Lava 

(purchased) 

-- 96% Al J 0 1 

Alumina 

Alaimag 798 

American Lava 

(purchased) 

- H5% Al 2 U 3 

Polygranular SIO^ 

GP-31 

Glasrock Products, Inc. 

lasting 

1160 't' 
1400 J C 
1640 J C 

Calcium Aluminate 

CA-25 

Alcoa 

Pressing and rolling 

60 “C 
1160 °C 
164 0' 1 1' 

/.ircotila 

Alaimag 475 

American Lava 



Sapphire 

To be aequirf . 





Mullito - Mullite is an aluminum silicate compound (SAlgOg • 2Si0 2 ) which 
is available in many forms. The crystalline material has a melting point of 
about 1850°C, but up to 30% glassy phase is commonly found in this material, 
which lowers the softening temperature. The thermal expansion and conduc- 
tivity of mullite is about 5 x 10 -6 /°C and 3.5 Btu/hr-ft°F, respectively, 
and its theoretical density is 3.26 gro/cc. Four types of mullite were 
initially obtained for this program: 

« Honeywell rolled MV-20 

• McDanel cast MV- 30 
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• Honeywell Hot-pressed Fiberfrax 


• American Lava pressed material 

McDanel mullite composition MV-20 was procured in a plastic extrudablo 
form. Test substrates (5 x 7 x 0.2cm thick) were produced by rolling the 
material over a smooth, oiled surface and controlled drying and firing to 
16i‘0°C. After rolling, the sheet material was cut to oversized 5cm x 7cm 
coupons. The coupons were dried betwoen flat plaster of paris blocks for 20 
hours at 60°C. They were then bisque-fired on mullite sagger plates to 
600°C and finally high-temperature-fired (on mullite) to 1640°C. Several 
larger 7cm x 12cra substrates and 20cm x 22cm substrates were successfully 
prepared of this material by the rolling methods. 

A second McDanel mullite composition, MV-20, was procured in fired tubular 
form. Eight Lost coupons were cut from this material. The MV-20 and MV-30 
materials typically contain 20 and 15% glassy phase after firing, re- 
spect i vely . 

Additional test samples of experimentally pressed mullite substrates were 
received from American Lava Division of 3M Inc. Rigid sheets of mullite 
fibers were also formed by pressing double or triple layers of Fiberfrax 
(1.2mm-thick Caroorundum sheet) under a load at 1400°C. Flat, very porous, 
rigid sheets lmra to 1.5mm thick with a density of 0.5 gm/cc were formed 
by this process. 

Alumina - Alumina substrates are the most common type of flat ceramic mate- 
rial used for electronic curcuitry. Pure A^O^ has a melting point of about 
2000°C and a thermal expansion and conductivity 6.7 x 10 - ^/°C and 10 Btu/ 
hr-ft-°F, respectively. The theoretical density of this material is 3.98 
gm/cc. 

Two types of alumina substrates were obtained from American Lava for this 
program, namely pressed 96% and 85% alumina. These are designated as 
Alsimag 614 and 798, respectively. A third type of alumina was fabricated 
at Honeywell by rolling, pressing, and doctor -bl ading a ca i c l um-uluminate 
bound alumina material. Table 2 gives four compositions that were formu- 
lated from Alcoa CA-25 calcium aluminate cement and tabular alumina filter. 
The consistency of the mixture was varied by addition of water based on the 
method of forming. 
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Table 2. Honeywell Ceramic Center formulation data. 


l ormulatton 
No, 

Materials 

Proportion* 
(‘ r > hv wt ) 

Water 
t hv wt ) 

harming 

Method* 

Drvtng 
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ttulil (*> 

II. T, * 
hiring 

<1160 i ) 

m. r, 

l iring 
<1400 t 1 

\\. V. 
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* 1040 i ) 
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\U'oa i‘A-25 

ftft 

20 

Pasting 

\ 
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30 
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Smooth compress surfaces with no sticking were obtained at 5000-psi pres- 
sure. No apparent warpage occurred during firing. Rolled substrates had 
some surface roughness, and slight warpage occurred .luring high-temperature 
firing. These materials were still somewhat porous but still had a density 
over 3 gm/cc. 

Zircon - Zircon is a zirconia silicate (Zr0 2 • Si0 2 ) compound with a melt- 
ing point of 1530°C and theoretical density of 4.68 gm/cc. The thermal 
expansion and conductivity of this material is 6.0 x 10 ®/°C and 3.8 Btu/ 
hr-ft-°F, respectively. 

Cordler lte - Cordierito is a magnesium aluminum silicate (2M^0*2Al 2 0 2 '5SiO 2 ) 
with a melting point of 1200°C; however, slightly different modifications 
may be as high as 1400°C. This material has a thermal expansion of 1.5 x 
10 _(3 /°C. 

Sample Cordierite substrates were obtained from the following two suppliers: 

• Minnesota Mining and Manufacturing Co. (3M - Alsimag 701) 

• Du-Co Ceramics Co. (A-3171) 
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A specimen of each material was fired at 1300°C with the results shown in 
Fig. 34. Both materials show leaching out of a lower melting-point con- 
stituent. The 3M material retained its shape, but the Du-Co material 
melted completely. The Du-Co material was therefore not evaluated for the 
silicon dip-coating application. 



Figure 34. Cordierite substrate samples. 


Silica - Silica in powdered fora was procured from Glasrock Products. This 
Glasrock grade "I" material (99% SiOg) was of two particle sizes designated 
GP31 (Fisher sieve 2.0) and GP71 (Fisher sieve 4.5). Two formulations were 
made (see Table 2). In each case, a slip was prepared by high-shear mixing 
of the solid ingredients and water. The material was then cast on a flat 
porous substrate. The castings were dried at 60°C, and subsequently high- 
temperature-f lred at either of two higher temperatures (see Table 1). 

A slight warpage of the cast sheets occurred after high-temperature firing. 
This was attributed to a nonunifora particle 3ize distribution through the 
sheet thickness. The fired substrates appeared to have adequate structural 
strength for subsequent processing. 
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With the exception of the cordierite substrate, all of the ceramics pre- 
viously discussed were dip-coated with varying degrees of success. All of 
the materials with the exception of American Lava's pressed mullite with- 
stood the thermal shock of being immersed into molten silicon with only a 
relatively short preheat cycle prior to their immersions. Of those listed, 
the VlcDanel MV-20 mullite plastic composition, which was rolled, dried, 
and fired at 1640°C by Honeywell's Ceramics Center, stood out as a strong 
favorite. This ceramic, as a substrate material, appeared to satisfy the 
important criterion of being defined a suitable material for the dip- 
coating process. Its cost-effective qualification also appeared to satisfy 
the long-range economic goals of JPL's Low-Cost Silicon Solar Array 
Project, and its thermal expansion coefficient (5 x 10 -6 °C) was reasonably 
compatible with that of silicon. No obvious stress or strain was evident, 
in either the silicon coating or the substrate. 

At this point in the program, the rolled, dried, and fired McDanel MV-20 
mullite substrates were adopted for further SOC studies. 

Mullite Substrate Development 

Later mullite substrates of similar composition to MV-20 were obtained 
from Coors Porcelain Co. They were designated as a S1SI composition and 
were produced by a cold-pressing and sintering process. 

The success of the coating experiments run with MV-20 mullite prompted an 
effort to see if the properties of the substrate could be optimized by 
compositional modifications. Based on discussions with M. Leipold of JPL 
and engineers at Coors Porcelain Co., substrates with experimental composi- 
tions based on variations of the 55 wt. % A1 2 0 3 composition were identified 
and were produced by Coors. The experimental compositions are listed in 
Table 3 and briefly discussed below. 

The data for the McDanel ceramic were obtained from a company brochure. The 
Coors analysis was based on analyses of the starting materials. The only 
notable difference between the two is the higher titanium in the S1SI. The 
high-mullite and high-silicon modifications were produced to determine 


♦For additional information, refer to reports titled "Development of Mullite 
Substrates and Containers" and "Silicon Sheet Growth Development for the 
Large-Area Silicon Sheet Task of the Low-Cost Silicon Solar Array Project" 
by D. G. Wirth and J. D. Sibold, Coors Porcelain Company Golden, Colorado. 
Work Performed Under Contract No. NAS-7-100-954878 


50 


I 


whether the mullite or the glass phase dominates the thermal behavior of the 
substrate during the coating operation. The boric acid modification was 
intended to be a body in which the glass phase is doped with B 0 °3 to 
produce a low-expansion glass and possibly act as a source of boron in the 
silicon coating. The open porosity modification is a body intentionally 
produced with higher porosity. The higher porosity will enhance the thermal 
shock resistance and will allow deeper penetration of the silicon coating 
into the substrate. The high-purity modification is a body produced from 
starting materials selected for their low-impurity content. As noted in 
Table 3, the iron, sodium, and titanium are noticeably low in this material. 
The electrically fused mullite is produced from a mixture of electrically 
fused mullite and silica powders. The impurity levels in this material are 
understandably low. The reducing- f ire modification is a standard S1SI 
mullite body that has been fired in a reducing atmosphere. The intent was 
to produce a material with hi B ii porosity that contained residual carbon 
in the form of unoxidized binder materials. 

A better idea of the phase relationships existing in fireclay refractories 
can be obtained from alumina-silica equilibrium phase diagram of Aksay 
and Pask shown in Fig. 35. Mullite is a compound that exists over a mea- 
surable range of compositions, and, according to Aksay and Pask, melts 
incongruent ly . Disagreements on the melting behavior of mullite still exist 
in the literature. Referring to Fig. 35, the equilibrium structure at room 
temperature of a composition containing about 55 wt . % A1 2°3 consists ot 
mullite crystallites in a eutectic matrix consisting of a finely divided 
mixture of Si0 2 and mullite. In reality, the liquid present at the firing 
temperature (approx. 1600 C) does not crystallize at the eutectic tempera- 
ture and instead forms a continuous glass matrix around the mullite 
crystallites. It should be further mentioned that the clay starting mate- 
rials do not contain any mullite but decompose during firing into a high- 
silica liquid into which the mullite crystallizes as firing proceeds. 

Structure and Properties of Ceramic Substrates - The structures and prop- 
erties of the substrates were measured at Honeywell and at Coors Porcelain 
Co. Microstructure, density, and modulus of rupture (MOR) determinations 
were made at Honeywell's Corporate Materials Science Center (CMSC), and 
density, firing, shrinkage, thermal expansion, and x-ray percent of mullite 
were determined at Coors. The Honeywell Ceramics Center provided firing 


2 A. Pask, Science 183, 69 (1974). 
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Code 


Material 


A1„0„ 
2 ■: 


sio„ 


mm 

UcDanel MV-20 

55.4 

lb 

to 

o 

■ 


52.6 

43.7 

B 

Coors SIS I 

57.6 

38.9 

ii 

UcDanel llV-30 

60.0 

38.0 



58.3 

38.9 






B 

High-mul 1 1 te 

67.2 

31.4 

C 

High -silica 

52.2 

44 . 7 

D 

Boric acid 

56 . 7 

38.3 

E 

Open-porosity 

57.6 

38.9 


High-pur 1 ty 

59.5 

39.1 


Electrically-fused 

77.5 

21.9 

H 

Reduclng-f Ire 

57.6 

38.9 


H 


ition of mullite-based ceramic substrates (in wt.%) 


Fe 2°3 


Na 2 0 T10 2 


WnO 


Coninerclal Compositions 



Esa 

0.7 

ES| 

0.5 

QQ 

0.74 

EM 

0,82 

IS 

- 

0.38 

1.30 

KM 

0.92 

0.25 

1 . 12 

0.2 

0.5 

0.08 

0.5 

0.19 

0.48 

0.10 

C 83 

0. 33 

- 

0. 13 

0.90 






Nominal UcDanel 




Coors analysis 

0.07 

0.01 

0.07 

Emission spectroanalysls 
Coors analysis 
Nominal UcDanel 
Coors analysis 

0.012 

0.01 

0.03 

Emission spectroanalysls 




Experimental Compositions from Coors Porcelain Co. 


.18 0.54 

.22 1.01 

0.13 0.24 0.90 0.24 1.10 1.50 

0.13 0.25 0.92 0.25 1.12 

0.21 0.13 0.34 0 05 0.30 

0.00 0.00 0.00 0.35 0.05 

0.13 0.25 0.92 0.25 1.12 



Coors 

analysis 

Coors 

analysis 

Coors 

analysis 

Coors 

analysis 

Coors 

analysis 

Coors 

analysis 

Coors 

analysis 







































A1 2 0 3 (WT PERCENT) 

0 10 20 30 40 50 60 70 80 90 



A1 2 0 3 (MOLE PERCENT) 


Figure 35. Alumina-silica phase diagram. Also shown 
are the compositions of the various sub- 
strate materials currently being examined 
(the code letters are explained in Table 3). 

shrinkage data on the MV-20 and MV-30 substrates. The properties measured 
to date are summarized in Table 4. The thermal expansion of these materials 
is in accordance with the good thermal shock resistance of fireclay refrac- 
tories. As shown in Fig. 36, the thermal expansion of the mullites is 
slightly greater than silicon. Thus, after cooling to room temperature, the 
silicon coating should be under a compressive status. Preliminary x-ray 
measurements (discussed later) of residual stress in the silicon have in- 
dicated this to be the case. 

The strengths of the MV-20 and boric acid modifications were the lowest of 
the group. The strongest samples were those containing greater amounts of 
muilite (i.e., S1SI and the high-mullite modification and the material 
fired in a reducing atmosphere). Both the strength and thermal expansion 
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Table 4. Properties of mulli te-based ceramic substrates. 


Code 

Material 

Density 

(gm/cc) 

% 

Porosity 

Thickness 

Shrinkage 

l 

Length 

Shrinkage 

(%) 

X-Ray 
% Mullite 

Equilibrium 
% llullite 

Thermal 
Expansion 
KT to 1000°C 
(%) 

MGR 
(kpsi ) 

| Commercial Compositions | 

• 

HcDanel MV-20 

2.50 

- 

14.0 

10.2 

- 

74 

0.495 

15.3 : 1.5 

KS 

Coore SIS! 

2.70 

0.32 

12.1 

16.7 

76.5 

81 

0.520 

23.7 i 1.4 

n 

McDanel MV-30 

- 

- 

- 

- 

- 

82 

0.470 

- 



Experimental Compositions from Coors Porcelain 

Co. 




B 

High-mul lite 

2.86 

0.22 

12.1 

12.7 

70.2 

95 

0.485 

20.9 1 1.2 

C 

High-si 1 ica 

- 

- 

- 

- 

56. 1 

74 

0.480 

- 

D 

Boric acid 

2.60 

0.19 

11.7 

13.7 

66.8 

80 

0.535 

15.6 i 0.7 

E 

Open-poros i ty 

2 . ? 5 

7.27 

11.6 

13.8 

66.4 

81 

0.525 

18.8 : 1.8 

mm 

High-purity 

7.70 

0.53 

14 .9 

19.5 

66.8 

84 

0.510 

20.1 1 1.0 

■ 

Electrically-f used 

2.6 to 
2.8 

- 

10.4 

11.3 

84 . 2 

87 

0.530 

19.6 -• 1.4 

Ll 

Heducing-f ire 

2.7 to 
2.9 

- 

12.5 

18.1 

66.4 

bi 

0 . d!5 

■ 

23.8 i 1.6 






f 








Figure 36. Thermal expansion of mullite 
compared with silicon. 


affect the thermal shock resistance. Materials with low expansion and high 
strength, such as the high-mullite modifications, are expected to have 
good thermal shock resistance. 

The average thermal exoanoion for mullite from room temperature to 1000°C 

3 

is 0.526o. This value is at the upper limit of those listed in Table 4 
and, as discussed later, indicates that the thermal properties of the glass 
phase rather than the mullite are dominating the properties of the sub- 
strates . 

The microstructures of McDanel MV-20 and Coors S1SI mullite refractories 
are shown in Fig. 37. The light regions are mullite crystallites and the 
dark phase is the glass matrix. The S1SI structure has more mullite and 
less porosity. This is consistent with the data in Table 4; the S1SI has a 
higher density and contains more AlgO^, which corresponds to a higher 
mullite content. The microstructures of most of the experimental composi- 
tions were similar to those shown in Fig. 38. The major differences were 
in the amounts of the phases (i.e., glass, mullite, and porosity). The 
microstructure of the electrically fused mullite body is shown in Fig. 38. 
The light areas are presumably mullite and the dark regions are silica. 


3 

T. G. Godfrey, Jr., MS Thesis, Clemson College, Clemson, S.C. 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOG3MPH 



(a) McDanel MV-20, specimen MR-114. 



(b) Coors SIS I, Specimen MC-1. 

Figure 37. Optical micrographs of mul.Mte 
refractory specimens. 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



Figure 38. Optical micrograph of an 

electrically fused mulllte 
substrate (500x). 

Characterization of Acid-Leached Mulllte Refractory - Most of the Impurities 
In a mulllte refractory are In the glass phase. In an effort to eliminate 
some of the Impurities prior to coating, a number of pieces of McDanel 
MV-20 were leached In hydrofluoric acid to dissolve the glass phase. A 
scanning electron micrograph of a fractured cross section of an MV-20 sub- 
strate Is shown In Fig. 39(a). This substrate was Immersed In hydrofluoric 
acid until about half the total cross section had been acld-leached . A 
higher magnification photo of the Interface region of the leached layer Is 
shown In Fig. 39(b), while a still higher magnification photo Is shown In 
Fig. 39(c) The aclcular morphology of the mulllte crystallites that pre- 
clpated during firing is readily apparent. The structure appears to contain 
only the continuous-glass phase and the mulllte needles. The results of 
EDAX analyses of regions In the center and at the edge of the leached 
sample show* d that, as expected, the aluminum/silicon ratio Increased and 
the potassium and titanium were lower in the leached region. The iron 
content was hardly altered by leaching. 

Although the leaching operation lowered the surface Impurity content, it 
also weakened the substrate. Leached substrates that were carbon-coated and 
dipped in silicon repeatedly fractured during dipping. 
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ORIGINAL PAGE 

BLACK ANO WHITE PHOTOGRAPH 




Figure 3i>. Scanning electron micrographs of a fractured surface 
of an MV-20 sample that had been Immersed In HF to 
leach the glass phase to a depth of approximately 
one-third the specimen thickness. 



• • 



The residual stress in the silicon layer can be estimated from an expression 

4 

used for stresses in glaze coatings applied to ceramics. For a thin silicon 
coating on an infinite substrate, the stress at temperature, T, is 

°S1 ■ E l T o - T l hi- a cerajnlc ) | l ' 3 * * ^ I < 2 > 

where E is Young's modulus and is assumed to be equal in the substrate 

and the coating, T is the temperature at which the substrate and coating 

are at a stress-free state, (a-.,- a , ) is the difference in linear 

Si ceramic 

thermal expansion coefficients of the coating and substrate, and j is the 
ratio of coating to substrate thickness. 

For the coatings applied in this work, the last term in the equation is 
negligible. Silicon is assumed to deform above 600°C (i.e., T ■ 600°C) and 

Q O 

E is approximately 15 x 10 psi. Toe room temperature stresses in the 
coatings estimated from the above range from about 4500 to 7500 psi, de- 
pending on the substrate material. A compressive stress in the coating is, 
of course, beneficial, since it raises the fracture stress of the coating. 


Some x-ray measurements have been made which have detected the residual 

5 

stress in the silicon layer. In these measurements the Bond technique of 
precision lattice parameter measurement was used to detect any stress- 
induced change in lattic* spacing. Since the technique can be used only on 
stressed single crystals or samples with large grain sizes, only large- 
grained samples were examined. The x-ray measurements showed the (110) 
lattice spacing in the silicon layer to be less than the <110> spacing in 
a stress-free single crystal by 65 ppm. This confirms the sign of the re- 
sidual stress. The <110> direction in the grain examined lay 20 degrees off 
the normal to the specimen surface. As a result, the strain (and stress) 
distribution parallel to the specimen surface could not be calculated be- 
cause the orientation of the crystal axes of the silicon grain with respect 
to the specimen surface are not known. A very rough estimate of the residual 
stress in the coating corresponding to a strain of 65 ppm can be made as 
follows. If we assume that the strain normal to the surface is 65 ppm, then 


a = Ee 


15 x 10 6 X x 10" 6 
0.3 


3250 psi 


(3) 


where Poisson's ratio, v, is assumed to be 0.3. This stress is on the order 
of that arrived at earlier on the basis of the thermal expansion mismatch. 

4 

W. D. Kingery, H. K. Bower and D. R. Wilraan, Introduction to Ceramics, 2nd 
Ed., Wiley, New York (1976) 1, p. 609. 

5 W. L. Bond, Acta Cryst . 13 , 814 (1960). 
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Figure 40 shows the thermal expansion of the experimental substrate compo- 
sitions from Coors plotted as a function of mullite content as determined 
by x-ray analysis. Also plotted on the figure are the expansions of pure 
mullite and vitreous silica. Each data point is identified by the code 
letters defined in Table 4. The results can be discussed in te r ms of the 
glass content and impurity content in each material. The addition of im- 
purities to vitreous silica increases the thermal expansion. The total 
impurity content as determined from Table 3 increases steadily through the 
series of samples with similar glass contents as follows: B, F, H, E, D. 



PERCENT MULLITE (X-RAY MEASUREMENT) 

Figure 40. Thermal expansion from 25°C to 1000°C of 
experimental substrate materials. 
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The boric acid modification material (sample D) has the highest purity 
content of all and the largest expansion. Sample A has the snme impurity 
content as samples E and H but contains more mullite. Sample C has a higher 
impurity content than B but contains more glass. The conclusion here is 
that the thermal expansion of a mullite-bas refractory is determined not 
only by the relative amounts of glass and mu. e in the material but also 
by the impurity content which primarily resides in the glass and alters its 
expansion. Impurities such as transition metals may be deleterious to the 
ultimate performance of the photovoltaic coating. Thus, if substrates are 
to be designed to be more thermally compatible with silicon, the most likely 
approach will be to alter the glass content rather than the glass chemistry. 

The boric acid additions were to produce a low-expansion glass to lower the 
thermal expansion coefficient to a value closer to that of silicon. This 
was not found to be the case. This problem was resolved with other com- 
position modifications discussed later. The porous body was produced to im- 
prove the thermal shock resistance and to allow more penetration of sili- 
con into the ceramic. The results were marginal. Th9 material fired in a 
reducing atmosphere contained a high porosity and a carbonaceous residue 
from unoxidized binder. Both were to enhance silicon penetration into the 
ceramic. Again the results were marginal. 

Modifications B through F and H survived the dip-coating operation but too 
few were dipped to make a valid comparison of the effects of the composi- 
tional modifications. In fact, only one solar cell was fabricated on one 
of these substrates. The efficiency was only 3.41%. The G-Mod material was 
the only substrate not to survive the dip-coating operation. This is at- 
tributed to the higher total mullite content in the body. The thermal 
expansions of these materials are discussed in more detail later. 

Semiquantitative analyses of the impurity contents of a number of these 
materials have been run and are given in Table 5. 

Notable differences include the lower impurity content in the F-Mod material 
and the higher titanium and copper and lower iron in S1SI compared with 
MV-20. As noted above, there have not been enough solar cells made on the B 
through H modifications to evaluate the impurity effects. The mechanical 
properties and thermal shock resistance of S1SI and of the Coors substrate 
materials A through H which were extensively characterized are discussed in 
the paragraphs that follow. 
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Table 5. Seraiquantltative emission spectrochemical 
analysis of mullite substrates (in wt . 1 ) . 


Sample 

Tl 

Cu 

Mg 

Fe 

Ca 

V 

Ni 

Cr 

Mn 

MV20 

0. 78 

<0. 01 

0. 20 

0. 89 

0 . 11 

0. 031 

— 

— 

0. 041 

S1SI 

1 . 1 

0. 071 

0. 29 

0. 68 

0. 070 

0. 051 

<0. 03 

0. 030 

<0. 03 

Open- 

1. 45 

0. 038 

0. 29 

0. 57 

0. 062 

0. 036 

<0. 03 

0. 026 

<0. 03 

porosity 

modification 










Reducing- 
f ire 

1 . 0 

0. 27 

0. 27 

0. 52 

0. 062 

0. 042 

<0. 03 

0. 28 

<0. 03 

modification 










High- 

purity 

modification 

0. 27 

0. 047 

0. 14 

0. 45 

0. 080 

<0. 03 

<0. 03 

N. D. 

<0. 03 


The fracture strength and thermal shock resistance of these materials are 
listed in Table 6, The fracture strengths were measured in four-point 
bending at room temperature on bars with a nominal cross section of 0.1 x 
0.1 inch over a span of 1 inch. The fracture strengths of the Coors C-Mod 
and McDanel MV-20 material were about the same. As previously determined, 
the A^Og/SiOg ratios and presumably, the glass contents, are also similar 
With the exception of the open-porosity and electrically-fused ruullite 
substrates, the fracture strengths of the rest of the Coors materials were 
about 20 kpsi. 

The critical quench temperature refers to a thermal shock resistance mea- 

g 

surement attributed to Hasselman. In the technique the room temperature 
fracture strength (usually in bending) of samples quenched from elevated 
temperatures is measured as a function of the quench temperature. When the 
quenching stresses are sufficient to propagate localized surface flaws in 
the ceramic, the room temperature strength decreases abruptly. The quench 
temperature causing the decrease in strength is an indication of how 
much thermal shock the material can withstand. Higher critical quench 
temperatures imply greater thermal shock resistance. The dat.; '‘sted in 
Table 6 were obtained on bars with a nominal cross section of 0.08 x 0.08 
inch tested in four-point bending over a span of 0.75 inch. Tests were run 


g 

D. P. H. Hasselman, J. .Am. Ceramic Soc. 52, 600 (1969). 
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Table 6. Room-temperature fracture strength and critical 
quench temperature for mullite refractory 
substrates . 


Code 

Material 

Strength 

(Kpsi) 

Critical Quench 
Temperature (“C) 

7-50 

MV20 

15. 0 t 2. 7 

340 

7-300 

MV20 

13. 8 ± 2. 3 

340 

A 

SI SI 

21. 7 ± 1. 4 

# 

B 

High-mullite 

21. 0 ± 1. 6 

300 

C 

High-glass 

13. 8 ± 0. 8 

360 

D 

Glass -property 
modification 

21. 1 ± 1. 4 

325 

E 

Open-porosity 

16. 6 ± 1. 8 

325 

F 

High-purity 

19. 4 i 0. 9 

325 

G 

Electrically-fused 

16. 2 ± 6. 5 

315 

H 

Reducing-fire 

22. 3 ± 1. 5 

315 


Not measured to date. 


on S1SI samples but the bars were larger in cross section. Since there is 
a size effect in thermal shock testing, these data ueed to be repeated 
on smaller specimens. The critical quenching temperature increases for 
smaller cross sections. Representative data showing the fracture strength 
as a function of quench temperature for MV-20 mullite are shown in Fig. 41. 
With the exception of the MV-20 and C-Mod material, the critical quenching 
temperature is around 325°C. It is slightly less for materials with higher 
mullite contents (B and G-Mod) and for the material fired in a reducing 
atmosphere. It is slightly higher for materials containing more glass, 
presumably as a result of their slightly lower thermal expansion. The ther- 
mal expansion data will be reported later in this section. 

Fracture toughness measurements have been made on MV-20 and on the Coors 
material A through H, excluding the C-Mod material. Briefly, the fracture 
toughness is an indication of the resistance of a material to fast, catas- 
trophic crack propagation and is usually denoted by the critical stress 
intensity factor, K IC is a measure of the stress at a crack tip during 

fast fracture in terms of the crack tip and loading geometry, the crack 
size, and the remote applied stress according to K IC = Y Op a, where Y is a 
geometrical constant, c F is the fracture stress, and a is the crack length. 
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Figure 41. Room-temperature fracture strength 
of MV-20 mullite. 


Kj C is a material property and is determined by measuring the load required 
to fracture precracked specimens with known loading and crack geometries. 

In this work, the constant-moment modification of the double-cantilever 

7 

beam technique, devised by Freiraan et al., is being used to measure Kj C - 
An advantage of this technique is that the stress intensity factor is 
independent of crack length. Thus, the fracture toughness measurement can 
be made by simply loading a precracked specimen to failure. Crack length 
measurements are not necessary. A schematic of the specimen and loading 
geometry is shown in Fig. 42. A side groove is cut into the specimen, as 
shown in Fig. 42, to guide the crack. 

K IC is calculated from 

K ic " (4) 


7 

S. W. Freiman, D. 


R. Mulviile and P. W. Mast, J. Mater. Sci. 8, 1527 (1973). 
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P P 6 



Figure 42. Schematic showing specimen and loading 
arrangement used in the constant-moment 
test. 6 is the deflection of the point 
of load application during the test. 

where I is the moment of inertia of one-half of the specimen, L is the 
moment arm, and t is the thickness at the base of the crack-guiding notch 
shown in the cross-section diagram in Fig. 42. Neglecting the effect of the 

3 

reduced cross section at the crack, I * wh /12. Note that the general form 
of both these expressions is K JC ■ AP c , where A is a geometrical term and 
P c is the applied load at fracture. In this test, Kj is independent of the 
crack length, a, and the measurements are correspondingly simplified, as 
noted above. 

Fracture toughness data taken to date are summarized in Table 7. 

Table 7. Fracture toughness measurements 
on mullite substrate materials. 


C ode 

Material 

Kj C (MNm' 3/2 ) 

MV20 

MV20 

1. 98 i 0. 23 

A 

S1SI 

2. 39 ± 0. 34 

B 

High -mullite 

2. 22 ± 0. 06 

D 

Class property 
modification 

2. 37 ± 0. 29 

E 

Open- porosity 

2. 16 t 0. 25 

F 

High-purity 

2. 23 t 0. 25 

G 

E lec tr ic ally - f used 

2. 10 ± 0. 31 

H 

Reducing-fire 

2. 30 t 0. 34 
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The fracture toughness of the mullite refractories listed in Table 7 are 
about the same. The lowest of the Coors materials is the electrically- fused 
mullite fired with fused silica. To put the numbers in Table 7 in per- 
spective, values of the fracture toughness of ceramics range from 0.75 

-3/2 -3/2 

MNm ' for soda lime glass to greater than 6 MNm 1 for hot-pressed 

silicon nitride. Fine-grained aluminum oxide has a fracture toughness of 
-3/2 

about 5 MNm ' . The fracture toughness of metals is typically an order of 

magnitude higher. 

An estimate of the flaw size initiating fracture in these materials may be 

1/2 

obtained from K^, » J p YA ' , where J p is the fracture strength and Y is a 

constant depending on the flaw and loading geometry. For penny-shaped 

g 

surface flaws, Y * 1.26. 

Flaw sizes calculated from K._. and o„ measurements on S1SI and MV-20 

IG r 

•materials are 160 and 230um, respectively. These are relatively large 
flaws and may indicate most failures in these ceramics originate from pores 
or pore clusters. 

Before the thermal expansion data are presented, some recent work on new 
mullite compositions should be pointed out. One of the goals of the work on 
the first series of Coors substrates (i.e., B through H-taod) was to see 
ii the thermal expansion of mullite-based refractories could be altered by 
compositional control. This was the reason that boric acid was added to 
S1SI to form the D-Mod substrates. The resulting borosilicate glass was 
expected to lower the thermal-expansion coefficient of the material. As pre- 
viously noted, this was not the case and the expansion of the D-Mod sub- 
strates was actually slightly higher than that of S1SI. We realized that 
substrates with thermal expansions closer to that of silicon would probably 
have to be developed when slotted substrates were used because the stress 
concentrations in the substrates due to the slots might initiate frac- 
tures in the substrates and coatings. This we later found to be true. 

Since fused silica has a lower thermal expansion than silicon and mullite 
has a higher thermal expansion, in principal, the thermal expansion of the 
two-phase refractory should be a function of the relative amounts of each 
phase. These are controlled by the A^O^/SiC^ ratio according to the phase 
diagrams shown in Fig. 35. With this in mind, J. Sibold and D. Wirth at 
Coors started formulating another experimental series of mullite composi- 

g 

G. R. Irwin and P. C. Paris, pp. 1-46 in Fracture , Vol. Ill, H. Liebowitz 
ed. , Academic Press, New York (1971) " 
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tlons by adding excess silica to the S1SI composition. This had alroady 
been done in the first series of compositions (i.e., the B and C-Uod mate- 
rials), but in the "high-silica" modification, only enough silica was 
added to S1SI to bring the Al 2 0 3 /Si0 2 ratio down to that of UV-20. In the 
next series, even more Si0 2 was added. The compositions of the second 
series of Coors experimental substrates are given in Table 8. 


Table 8. Composition of the second series of 
Coors experimental substrates. 


Code 

Material 

Compos i t i on 
wt. % A1 2 0 3 

Al 2 0 3 /S10 2 
Hat io 

A 

S1SI 

57.6 

1.48 

C 

"High-silica mod" 

52.2 

1.17 


S1SI ♦ 10 wt. % excess Si0 2 



I 

S1SI + 18 wc. % exciss Si0 2 

48.8 

1.01 

J 

81SI + 25 wt . % excess SiO^ 

45.8 

0.89 

K 

b IS I ^ J8 wt. % cxcoss SiC>2 

41.4 

0 . 74 


We see in what follows that the material with the highest silica content 
(i.e., K-Mod) has properties that make it much more applicable as a sub- 
strate for SOC than any others tested to date. 

The microstructures of the A, I, and K-Mod materials are shown in Fig. 43. 
The darker phase is the continuous glass phase. The amount of glass phase 
is seen to increase as the AlgO^/SiOg ratio decreases. The pore size also 
appears to increase with increasing silica additions. Measurements have 
not been made of the densities at this time. The mechanical properties and 
thermal shock resistance of these materials have also not been determined. 
Thermal expansion data have been provided by Coors and are discussed below. 

The thermal expansions of the substrate materials produced by adding excess 
silica to the S1SI composition are shown in Fig. 44. The data are plotted 
as the linear thermal expansion from room temperature to 700°C as a 
function of the A^O^/SiOg ratio. The residual stresses between the silicon 
layer and the substrate are presumed to be negligible above 700°C. The 




67 



ORIGINAL PAG? 

BLACK AND WHITE PHOTOGRAPH 




(a) S1SI 


(b) S1SI ♦ 18 wt. % excess SiO, 



(c) S1SI + 38 wt. 1 excess SiOg. 

Figure 43. Photomicrographs of various mullite bodies 
(magnification 150x). 
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data corresponding to S1SI. I and K modifications are indicated as such on 
the figure. The thermal expansion of silicon is also included in the figure 
Vote that as the silica content increases, the thermal expansion decreases 
m the manner expected and that the thermal expansion of a mullite re- 
fractory with an Al 2 0 3 /Si0 2 of approximately 0.45 should be equal to that 
of silicon. 



Figure 44. Thermal expansion of S1SI mullite 
containing excess silica. 


A silicon-coated substrate with this composition should not show any re- 
sidual stresses between the coating and the substrate. The residual stresses 
in the silicon coating can be estimated by 


EAT a 


Si 


- a 


ceramic 


(5) 


where Young's modulus, E, is assumed to be equal in the substrate and the 
coating, AT is the temperature difference when the composite cools below 
its stress-free state, and a Si and are the linear thermal-expansion 

coefficients of silicon and the substrate, respectively. If we assume E is 
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15 x 10 6 psi, a gl (RT - 700°C ) - 3.7 x 10" 6o C _l , and a stress-free state 
at 700°C (i.e., AT ^ 675°C), the residual stresses in the silicon layers on 
SIS I and the K-Mod substrates are: 

a (RT - 700°C ) 

. o 

Material (°C _1 ) (psi) 

S1SI 4.7 x 10~ 6 10,000 

K-Mod 4.0 x 10~ 6 3,040 

Note that with even a small thermal mismatch, significant thermal stresses 
can develop. These stresses can cause the slow growth of flaws. When sub- 
critical flaws reach critical size under a constant applied stress, large- 
scale fracture may occur. 

The K-Mod substrates proved to be very comfortable with the SOC process in 
every respect and innumerable satisfactory SOC layers, large and small, 
were achieved using them. Minor modifications to the K-Mod substrates were 
later made and designated as N and 0 mods. While no in-house characteriza- 
tion of these latest modifications was performed, *-hey too proved to be 
completely satisfactory for SOC purposes. 

Slotted and Large-Area Substrate Development - The ceramic substrates used 
in SOC solar cells will necessarily have to be produced by a continuous 
process in order to meet the economic goals of the program. To this end, 

J. Sibold and D. Wirth at Coors have been adapting their roll compaction 
process to the manufacture of mullito substrates for the SOC process. This 
effort was quite successful. Large, 10cm x 100cm S1SI substrates were 
produced from roll-compacted tape. These substrates were 1.0mm thick and 
are some of the largest single substrates produced by a tape process. 

Roll-compacted tape has also been used to produce slotted substrates for 
back-contacted SOC solar cells. In this p ocess, slotted substrates with 
dimensions shown in Fig. 45 and 46 were produced by punching the slot 
patterns in green, unfired tape with plastic die sets. The properties of 
silicon coatings put down on slotted substrates are discussed elsewhere 
in thi; report. 

A final comment is made pertaining to the mechanical integrity of slotted 
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Figure 46. 10cm x 100cm slotted substrate 
produced at Coors. 
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substrates. It was noted that slotted substratos have more of a tendency to 
crack after coating with silicon and during subsequent processing than un- 
slotted substrates. This was expected and results from the tendency for 
subcritical crack growth to occur in ceramics and from the stress concentra- 
tions associated with a slot. Since the thormal expansion of the substrate 
is greater than silicon, the subst ,te la in tension at room temperature. 

The stress concentration at the end of a slot can raise the local stress to 
values larger than those listed above. The fracture stress of a ceramic 
depends on the fractive toughness and the flaw size. If flaws in the sub- 
strate can slowly grow under a constant: residual stress, they may grow 
to a size where they become critical and initiate fast fracture. This is 
the origin of delayed failure in glasses and other ceramic materials and 
has been observed by us in silicon-coated mullite. This problem h<*s been 
minimized by using the K-Mod and G-Mod substrate materials, but is still oc- 
casionally observed. 

SUBSTRATE CARBON COATING EFFORT 

As discussed in the Introductory section of the report, the carbon or 
graphite coatings play a k ?y role making SOC possible. 

Initially, the substrate carbonization was performed by simply hand-scrub- 
bing high-purity graphite onto the substrate. It was originally believed 
that this graphite coating, in addition to making the substrate wettable to 
molten silicon, could perhaps serve as a back contact to silicon coating 
which otherwise was electrically insulated by the ceramic substrate. It was 
also not known to what extent, if any, the type or nature of carbon coating 
would ultimately affect the crystalline structure of the subsequent silicon 
layer . 

The scrub-on-coatings proved to do an excellent .job in producing good 
adhering, smooth, large-grain layers, but they did not appear to be cost- 
effective in a scale-up process. For this reason, early in the program, 
the following alternate carbon coating methods were investigated: 

• Applying carbon soot to one surface of the substrate, using an 
acetylene torch. 

• Applying black Aqua-Dag to one surface of the substrate and subse- 
quently firing it at 900°C in an argon atmosphere. 
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• Applying a liquid coating of pure acrylic which had been dissolved 
in 1-2 Dichloroethane (CgH^Clg) to one surface of the substrate 
followed by a 900°C firing in an argon atmosphere. 

• Pyrolyzing a layer of transparent 3M-brand Scotch tape in a 900°C 
furnace having a nitrogen atmosphere. 

Of these methods, only the rubbed-carbon technique produced i.niform coat- 
ings of silicon. Except for the Aqua-Dag which flaked off the substrate 
following the 900°C firing, the remaining methods produced only spotty 
wetting of the silicon. 

Later in the program, five different 1 iquid-carbon-coat j r,g compounds were 
purchased from Dylon Industries. These coatings had the viscosity of latex 
paint and were water-soluble. A t .al of 22 substrates were silicon coated 
that were first coated with the various Dylon compounds. It was found 
necessary to oven dry these coatings, as residual moisture caused them to 
blister or "boil" off when dipped in the molten silicon, as shown in 
Fig. 47. It was also noted that the granular nature of most of these coat- 
ings promoted nucleation which decreased the grain size, and that carbides 
formed throughout the silicon layer. However, in cases where the granular 
surface was smoothed by wiping off the excess graphite and the coating 
was oven-dried, the Dylon coatings gave large grain size similar to the 
hand-rubbed graphite coatings. Unfortunately, our chemical analysis of 
these coatings has revealed that they contain from 1 to 2 % impurities and 
were unsuited for our application. 

Another somewhat novel method of carbon coating was also investigated - 
plasma decomposition of propane gas in an RF sputtering apparatus. The 
carbon is deposited uniformly and its purity reflects that of the gas. The 
morphology and grain size of the silicon film formed on such substrates 
is comparable to those on hand-rubbed graphite coatings. 

Several methods of incorporating the carbonization with the ceramic pro- 
cessing were also examined. These were rolling graphite or silicon carbide 
particles into one surface of the substrate and then firing in an inert 
atmosphere or just firing the ceramic in an inert atmosphere so that the 
residue from the organic binders would form a carbon film. None of these 
approaches proved satisfactory for silicon coating. The plain ceramic fired 
in an inert atmosphere shattered in the temperature gradient above the 
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Figure 47. Blistered and spotty silicon film formation as tbe -esult 
of incomplete drying of Dylon carbon coatings (1.7x). 
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melt before it could be dipped. The substrate with silicon carbide de- 
veloped large blisters and it appeared that the silicon did not wet the 
surface. The substrate with the graphite particles appeared to have been 
wet by the silicon, but either there was not sufficient carbon or the 
particles did not stay embedded to permit the growth of a silicon layer. 

Some ’-ery promising experiments were conducted using colloidal carbon 
suspensions obtained from Acheson Colloids Co., Port Huron, Michigan. Four 
different mixtures were tested: Dag 154; Electrodag 154; Electrodag 155; 
and Electrodag 502. The 154 mixtures are graphite in isopropyl alcohol 
with a cellulosic binder. Eloctrodag 155 is a mixture of graphite in 
trichlorethylene with a cellulosic binder. Electrodag 502 is a graphite 
and carbon-black mixture in methyl ethyl ketone (2-butanone) with a 
f luroelastomer binder. Drying posed no problems and coating thickness was 
relatively easy to control. 

Spectrographic analysis shows that the dags should be suitable for solar 
cell application. The results of the analysis are shown in Table 9. 


Table 9. Spectrographic analysis results of colloidal 
carbon suspensions from Acheson Colloid Co. 
Values are in ppm. 


Element 

Dag 154* 

Dag 154 

E 'dag 502 

E 'dag 154 

E 'dag 155 

ar 

88 

ND < 10 

110 

20 

33 

Mo 

100 

ND < 33 

34 

ND < 33 

ND < 33 

Sn 

52 

ND < 10 

ND < 10 

ND < 10 

ND < 10 

Ca 

180 

43 

68 

39 

56 

V 

33 

ND < 33 

ND < 33 

ND < 33 

ND < 33 

Ti 

24 

8 

53 

7 

8 

Nt 

15 

2 

2 

2 

2 

Si 

220 

35 

72 

73 

51 

Pb 

35 

ND < 33 

ND < 33 

ND < 33 

ND < 33 

Mg 

4 

2 

14 

2 

3 

B 

3 

ND < 1 

37 

ND < 1 

3 

Mn 

4 

2 

3 

2 

3 

Na 

200 

ND < 100 

270 

ND < 100 

ND < 100 

| Fe 

1500 

260 

370 

285 

290 


Jjt 

Older stock. 


/ 5 
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The growth morphology of the silicon was not adversely affected by the 
dags. Metal lographlc examination of cross sections of SOC which were 
carbonized with Dag 154 showed excellent wetting of and bondi.'ig to the 
ceramic surface by the silicon, as shown in Fig. 48. 



Figure 48. Silicon-ceramic interface 
carbonized with dag 154 
(500x) . 


Vitreous carbon coatings were applied to several mullite substrates by the 
Tylan Corporation, Torrance, California. Their process is proprietary, but 
they report the thickness of the coating to be approximately 7pm. While 
our primary purpose for these coatings was to prevent substrate dissolution, 
upon dip-coating witu silicon, a secondary interest was related to the 
possibility that these coatings may provide a conductive silicon-ceramic 
interface. With this pu’-t se in mind for all of the carbon coatings, the 
electrical resistivity oi several carbon coatings was measured. 

There are several important aspects of the electrical resistivity of carbon 
layers. First, the electrical resistance can be used to evaluate the 
uniformity ami reproducibility of the carbon-coating process. Secondly, the 
sheet resistance of the carbon may affect the resistivity measurement of the 
silicon. Thirdly, it would be desirable to use the carbon as a back contact 
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to the silicon to decrease parasitic series resistance of solar cells. 
Table 10 summarizes our measurements on sheet resistance of carbon layers. 

It is clear from Table 10 that with the rubbed and the dag-borosilicate 
coatings, graphite is almost entirely removed by the melt. The dag only 
layer is partially removed, but has the appearance of being continuous ani 
intact after removal of the SOC by etching or spalling. 


Table 10. Typical carbon s eet resistivities. 


Carbon Type 

As Deposited 
(ohms /square) 

After CP-4 Etching to 
Remove SOC Layer 
(ohms /square) 

Graphite rubbed on rolled 
raullite (MV-20 ) a 

10 to 25 

250 to 600 

Graphite rubbed on pressed 
raullite (American Lava) 

100 to 150 

IK to 60K 

Vitreous carbon (Vitregraf) 
on MV-20 

2.5 to 5.0 

2.5 to 5.0 

Da; (colloidal carbon in alcohol) 
on MV-20 

2 OK to 40K 
(un fired) 

120 to 600 


20 to 40 
( fired) 

70 to 450 

(where SOC spalled off) 

Dag and borosilicate glass (1:1) 

- 

5K to 1M 


a Th.e thickness of the layer based on weight gain is 7um, which 
would give a sheet resistivity of about 2.5 ohms/square. 


The vitreous carbon layers, on the other hand, even after 1 hour >f contact 
with the melt, are almost wholly intact. Since the sheet resistivity of the 
vitreous carbon is only 2 to 5 ohms/square, it is much lower than the 
sheet resistivity of the p-type silicon layer (50 to 500 ohras/square) and 
could seriously affect the sheet resistivity measurement if there is low- 
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resistance contact between p-silicon and vitreous graphite. The nature of 
this contact was therefore investigated. 

Samples were mounted in a metallurgical mount, and cut and polished to 
show che cross section. The microscope showed that there was an interface 
layer that was harder than either the carbon or the silicon - a strong 
suggestion of SiC formation. 

The samples were then removed from the potting compound and the cross 
section was examined using scanning Auger microscopy at Physical Elec- 
tronics, Inc. Some of the results are shown in Fig. 49. A complete Auger 
spectrum of the interfacial material showed only silicon and carbon, 
and no impurities (at the 0.1% level). It is interesting, however, that 
the shape of the carbon peak and the silicon-to-carbon stoichiometry 
varied from that of a polycrystalline silicon-carbide standard. It is not 
known whether these discrepancies are significant, but they suggest that. 
the form of silicon carbide is different from the usual form. 

At any rate, several comments on rue silicon-carbide interface are in 
order. The SiC appeared thickest at the top of sample MR-78, which had 
been, in the melt for an hour. Since it was much thicker than the carbon 
coating, the SiC must be formed from carbon particles which float to the 
top of the melt. There was no detectable SiC at the top of sample MF-106 

t 

which had been in the melt a minimum amount of time. Away from the surface, 
the SiC If .ckness apDeared about the same in the two samples, varying 
from 2 to 10pm. Typically, the SiC was somewhat thinner than the vitreous 
carbon coating remaining. 

The samples were etched in Sirtl etch after polishing, which revealed the 
presence of cracks. The region shown in Fig. 49 was selected because 
the SiC was of maximum thickness and because it showed some cracks in the 
silicon and at the silicon-SiC interface. These cracks may be responsible 
for the nonuniformity of the electrical contact, or it may be due to 
variations in SiC thickness. 

A goal was to improve the conductance (increase thj thickness) of the 
carbon layer used to wet the ceramic and to obtain ohmic contact between 
the carbon and the p-type silicon so that the carbon contributes signifi- 
cantly to decreasing the series resistance of the solar cells. Two addi- 
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Figure 49. Scanning Auger micrographs of the silicon- 
vitreous graphite interface at 2500x. Photo 
at left is a secondary electron image; the 
center photo is the carbon image; at right 
is the silicon image of a section of sample 
MR-106. Silicon is at the top, carbon is at 
the bottom, and the silicon carbide is in 
between. 
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tional goals were: (1) to investigate the use of doped carbons as diffusion 
sources for doping of the back surface of silicon to provide a BSF effect, 
and (2) to investigate the use of the glassy carbon (GC) as a protective 
barrier against unwanted impurity diffusions from the substrate into the 
silicon layer. 

An important part of this effort was to have control of the carbon layer. 
For this reason, we attempted to develop a method of producing doped GC 
layers. We had previously demonstrated that "Vitregraph" GC coatings made 
by Tylan Corporation showed very little erosion by the molten silicon. 

This is important if the vitreous carbon is to be used as a protective 
barrier . 

In Table 11, sil icon-to-carbon contacting is indicated in the dipped GC- 
coated substrates by a lower measurement of 8 to 40 ohms/square. In 
contrast, the hand-rubbed Ultra-Carbon control sample is 200 ohms/square. 

On some of the dipped samples, the silicon was etched off with CP4 etch and 
considerable erosion o' the vitreous carbon layer was observed. There 
was evidence of the boron-doped carbon doping the melt. This is indicated 
by contrasting the p y of the silicon layer on the doped-carbon samples 
and the p g of the silicon on the hand-rubbed Ultra-Carbon sample dipped 
after the doped-carbon samples. 


Table 11. Sheet resistivity of silicon coatings 
(first dip run). 


Substrate 

Silicon Layer 

o* <n/n) 

Hand-rubbed Ultra-Carbon control before 

200 

( 00° fired GC (4 samples) 

8 to 40 

900° fired boron-doped GC (dipped) near 
end of run 

to 

Hand-rubbed Ultra-Carbon dipped after 
boron-doped samples 

17 


In a second dip run, eight GC-coated substrates with thinner carbon coat- 
ings (average p^ = 52 ohms/square were dipped. The silicon coated all of the 
substrates very well. Sheet resistivity measurements for the sequence are 
summarized in Table 12. All samples were dipped with the same pull-rate, 
resulting in similar silicon coating thicknesses. 
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Table 12. Sheet resistivity of silicon coatings 
(second dip run). 



Silicon Layer 

Substrate 

e d (nr-) 

Hand-rubbed Ultra-Carbon control before 

120 

Unheat-treated GC (3 samples) 

avg. = 130 

Boron-doped GC (5 samples) 

avg 3. 4 

Hand-rubbed Ultra-Carbon control after 
GC samples 

2. 1 

I 


Comparison of the first two lines of Table 12 shows no evidence of electri- 
cal contact between the silicon and carbon. This is possibly due to the 
high erosion of the thinner carbon layer by silicon. The bottom two lines 
indicate boron doping of the melt by the GC for the same reasons as pre- 
viously mentioned. Diodes were made from this dip run and results are 
listed in Table 13. 

Table 13. Summary of SOC with glassy carbons. 


Diode 

No. of 
Diodes 

Avg» Active 
2 

Are* (cm ) 

V 

oc 

(V) 

.1 

SC 

(mA/cm ) 

Fill 

Factor 

B 

Max 

Eff 

(%> 

Comments 

120-lA-Qll 

5 

0. 04 6 

0 47 

22 2 

0. *568 

7. 2 

8. 2 

Hand-rubbed 
carbon control 

120-3B-QUA 

6 

0 070 

0. 50 

21.6 

0 679 

7. 25 

8. 4 

GC 

120-3B-Q11B 

6 

0. 085 

0 50 

22 0 

0. 672 

7 4 

7. 7 

GC 

120-3A-Q1! 

! 

1 05 

0 50 

19. 5 

0. 619 

6. 0 

6. 0 

GC 


The cell measurements show that the GC coatings can produce performance at 
least as good as that produced with hand-rubbed carbon. 

It was expected that the erosion resistance and the impermeability of the 
GC layers would be increased by high-temperature firing. For this purpose, 
an induction furnace for high-temperature firing of vitreous carbon 
coatings was set up. Some problems arose with the first few samples which 
were fired. The first coating fired successfully, but succeeding coatings 
tended to flake off. This problem was later rectified by thinning the 
polymer and laying down a much thinner coat for the second and succeeding 
coatings . 
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In summary, GC layers have been demonstrated to act effectively as wetting 
agents ani as diffusion sources for doping silicon layers. Firing at 
900°C does not appear to give adequate erosion resistance and impermeabil- 
ity. although such films are adequate to produce good cell performance. 
High-ter.iperature firing is needed to reduce the erosion effect of the 
molten silicon. 

At this point in the work with glass-carbon coatings, the work was tempo- 
rarily discontinued due to the loss of personnel who were performing the 
work . 

With respect to the behavior of a carbon coating as a wetting agent, it 
appeared that any means of applying a smooth, uniform, dense coating of 
carbon of suitable thickness and purity will promote the coating of ceramic 
with silicon. Because of the ease of application of the dag mixtures, they 
were eventually adopted as a final method for applying the carbon coating. 

SHEET SILICON GROWTH EFFORT 

Introduction 


Sheet silicon growth was obtained during the investigation, by two very 

different coating techniques. These techniques were, as previously discussed, 

dip-coating and SCIM-coating methods. In this section of the report it will 

be shown, nowever, that each o, these techniques produced a final product 

which was structurally very similar. In an effort to fully understand the 

growth mechanisms of these two growth processes, considerable analytical < 

support was given to the experimental effort. Since day-to-day experi- 

f 

mental progress can only be judged through analysis of results, consider- I 

! 

able material characterization was also performed. This subsection of j 

the report summarizes the various phases of analytical and experimental 
growth studies along with the numerous characterization studies required ’ j 

to pursue this investigation. j 

f 

Thermal Analysis of Sheet Silicon Growth Processe s 

Introduction - The growth rate of sheet silicon appears to be limited by ^ 

heat-extraction efficiency. Our economic analyses indicate that a through- 
put speed of about 14 cra/min is required to meet the price goals, and, 
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to date, our best solar cells have been grown at one-fourth this speed. We 
have dedicated a large amount of effort to increasing our understanding 
of the thermal mechanism of our growth process with tho aim of attaining 
the necessary high throughput. 

The present SOC growth appears to proceed from a liquid-solid interface 
which is roughly perpendicular to the pull direction (see Fig. 50). We 
have been using the word "symmetric" to refer to this type of growth. 




Figure 50. Vertical growth of sheet silicon, (a) Unsupported 
growth, (b) Supported growth. In the pi .jsent dis- 
cussion, the effect of a die to stabilize the 
unsupported growth is neglected. 
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Symmetric growth may be unsupported, as lu the edge-defined, film-fed 
9—13 14 

growth ( EFG) or the web-dendritic ribbon processes, or supported on 

15 1 0 

a substrate, as in the silicon-on-ceramic (SOC) ' or r ibbon-agaiast- 
17 

drop (RAD) methods. In either case, the attainable growth rate is 
limited by heat extraction efficiency. The steady-state analysis of unsup- 
ported silicon growth is greatly simplified by the symmetry of the struc- 

1 0 

ture, leading to an uncomplicated solution to the heat equations. Im- 
portant by-products of the analysis are expressions for pull speed versus 
sheet thickness and for the shape of the crystallization or freezing front. 
This simplicity does not exist in the case of silicon grown on ceramic, 
and numerical methods must be chosen to solve the heat equations. The 
first part of this subsection is an analysis of supported symmetric growth, 
for the SOC case in particular. It is found from t.he symmetric growth 
studies that the economic throughput target of 12 to 18 cm/min will prob- 
ably net be met in symmetric growth with only passive (radiative) cooling. 
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We have also begun a theoretical study of the so-called wedge-shaped or 
asymmetric growth as an alternative to the more conventional, symmetric 
growth. Asymmetric growth offers high throughput potential and possibly 
other advantages such as sweeping impurities out of the active solar-ceil 
region. However, very little is known about this growth mode, particularly 
the thermal and mechanical design factors which must be incorporated in 
the growth apparatus to produce and maintain a stable wedge-shaped growth 
region. These factors are analyzed in the s icond part of this subsection. 
Only unsupported asymmetric growth has been considered here. The exten- 
sion to substrate supported growth is left for a future investigation. 

Part I: Supported Symmetric Growth - The thermal analysis of supported 
symmetric growth does not assume any particular orientation in the gravi- 
tational field but a vertical orientation is elected for purposes of 
illustration (Fig. 51). The important considerations are the rate of re- 
lease of the latent heat of fusion and the thermal properties of the SOC 
beginning at the location of the liquid-solid interface and extending down- 
stream. These consist of the neat loss mechanism, assumed to bo primarily 
radiative, at the front silicon and back ceramic surfaces, and conductive 
as well as convective heat transport through the silicon and ceramic 
materials. Except for possible gaseous convection, which is not considered 
here, these properties are independent of orientation in the gravitational 
field. The analysis applies to both dip-coated and SCIM-coated material, 
except for the details of the mensicus shape which do not affect the main 
result of the analysis - the growth velocity versus silicon thickness 
tradeoff. The dip-coating geometry was selected to represent the upstream 
boundary conditions in our thermal model of the process. The essential 
features of dip-coated SOC are shown in Fig. 51. Ceramic substrates are 
pulled vertically from the silicon melt at constant speed, v. The freezing 
front appears at a height y„ above the melt. The value of y„ at zero 

" * ig 

pulling speed may be calculated from hydrostatic considerations. Moreover, 

this appears to remain unchanged in actual steady-state growth for pulling 

19 

speeds well into the millimeter-per-second range. Above y the tempera- 
tures in the structure begin to fall toward some effective radiation 
temperature, T , of the environment. 
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Figure 51. Cross-sectional view of SOC growth. The meniscus 
shape is based on hydrostatic equilibrium. Radi- 
ative background temperatures Tq and Ti are seen 
above and below the liquid-solia interface. 


A possibly different radiation temperature, T^, is seen by the portion of 
the structure below the freezing front. In the present study, only radiative 
heat removal was considered aid the temperatures and T q were assumed 
uniform over 0 < y < yp and y^ < y < «, respectively. The simple two-part 
background is used to give a crude approximation to the effect of a heat 
shield located at y p . More general heat loss functions are easily incorpo- 
rated and may be investigated later. As was true in the unsupported case, 
the rate of heat extraction by radiation must balance the rate of removal 
of sensible heat from the melt (by convection) plus the rate of release of 
latent heat at the freezing front. In the present case, most of the sensible 
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heat is carried in the ceramic owing to the greater thermal mass. 


Hea t Equations — Two zones are identified because of the different thermal 
properties of molten and solid silicon. In the zone <_ y < «, the steady- 

state heat equations in the ceramic portion follow from two-dimensional 

20 

heat balance, including convective terms: 
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The two-dimensional model is appropriate because the width of the structure 
(z-dimension ) is assumed large compared with the thickness, t„ + t . Phy- 
sical parameters are identified in Tabl? 14. For the silicon portion, a 
one-dimensional temperature profile, U (y), is used, equivalent to 
assuming the isotherms are horizontal. This assumption can be checked by 
estimating the inclination of the isotherms in the solid silicon portion. 
The steady-state equations in the solid silicon (y > y F ) are 


t 


o 



I <*U 
k -3-2 
1 o dy 


p C v 
o o 



f Q (y) 


(9) 


f Q (y) 


e o °{[ V *>] 4 




T (t c ,y) 


( 10 ) 


U o ( y) 3 T(t c ,y) (11) 

U o (y F ) " T F (12) 

The extraction term, f (y), includes both radiation at the external surface 
and conduction across the internal surface into the ceramic. 


In the zone 0 < y < y F> the equations for the ceramic portion are similar 
to Equations (6) and (7), except for the radiation temperature, T^ , in 
place of T q r ig. 50) and a new boundary condition at y = 0 instead of 


20 

H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids . 2nd Ed. , 
Oxford Univ. Press (1959). 


87 


C-'L 


Table IT. Nomenclature and numerical data. 


Symbol 

Description 

Value/ Units 

C c 

Specific heat of ceramic 

1.04 (J/g °IO 

C o C 1 

Specific heat of silicon 

0. 954 (J/g K) 

.1 % <y> 

o 

Silicon heat flux density (y "» y^,) 

(W/cm 2 ) 

.ij(y) 

Silicon heat flux density (y < y^) 

(W/em 2 ) 

k 

c 

Heat conductivity of ceramic 

0.041 (VV/cm °K) 

k 

o 

Heat conductivity of solid silicon 

343/ I J (W/cm °K) 
o 

k i 

Heat conductivity of liquid silicon 

0. 6 (W/cm 'K) 

i. 

Lutent heat of fusion of silicon 

1802 (J/g) 

T(x.v) 

Ceramic temperature 

( 'K) 

T 

o 

Downstream (y y^,) ambient temperature 

(K) 

T i 

Upstream (y < Vj.) ambient temperature 

( ’K) 

T M 

Melt temperature 

(°K) 

T 
1 F 

Silicon freezing temperature 

1685 ( K) 

t 

c 

Ceramic thickness 

(cm) 

t 

o 

Solid-silicon layer thickness 

(cm) 

t t (y) 

Liquid- silicon layer thickness 

(cm) 


Silicon temperature (y * y^) 

C’K) 

U 1 <y> 

Silicon temperature (y s y^,) 

(°K) 

V 

Pull speed 

(cm/ sec) 

y F 

Total meniscus height above melt 

0. 5 or 0. 69 (cm) 

5 

Numerical grid spacing 

0. 02 (cm, 

p c 

Ceramic density 

2. 25 (g/cm 3 ) 

p o' °1 

Silicon density 

2. 33 (g/cm 3 ) 




. 3T 

3 

k c Tx 

* 77 


Eq, 

(15): 


h 

3T\ r „ 3T . o 

7y j - P c C c V 77 0 

(13) 


[T(0,y)] 4 - Tj 4 } 

(14) 



(15) 


T( x , 0 ) - T u 

m the liquid-silicon portion, the variable thickness ^ must be taken into 
account. The following is derived based on the assumption of horizontal 
iso the ms see next subsection for derivation of Eq. (16) 


. ; dU i\ dU i 

37 | c i k i it) " p o vt o c i 37" “ 

f x ( y ) 

(16) 


[ u i.(y )] 4 - T i 4 j " k c h T(t c y) 

(17) 

U 1 (: ) - T(t c .y) 


(18) 

u i^ y f ^ 


(19) 


A final condition expresses tie fact that the freezing front is the 
source of a heat flux density (vertical component) equal to P o vL: 


-k. 


d 

3y 


( y p ) 


p vL 
o 


- -k 


d 

37 


( y p ) 


( 20 ) 


Derivation of Eq . (16) — Refer to Fig. 62. We assume horizontal isotherms 
so that temperature can he denoted by U^y). Vertical mass flow, inde- 
pendent of y, is equal to P 0 vt Q . If t^y) denotes local thickness, vertical 
heat flux. H(y), per unit width is 


H(y ) - p o vt Q C l U 1 - t 1 k 1 dU l /dy ^ 

Denoting by f^y) the heat flux density lost by radiation at the curved 
surface and by conduction to the ceramic, we then have 


dli/dy f x - 0 


or 
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Figure 52. Diagram of meniscus region, 


dU, 


dU, 


h ri k i3T ' p o vt o c i W m f i (y) 


For t^(y), we use the result from Harrill, et al. . 


14 


( 22 ) 


( t i - t Q )/K - cos 


h -1 ( 2K/y ) - cos h -1 (2K/y p ) + 


4 - (y F /K) 2 l 


1/2 


-[4 - (y/K) k y 


(23) 


where K/y, - [2(1 - sin 9)]" 1/2 - 0.79 if the contact angle, 0. is 11°. 
I* 


Finite-Difference Mod el — A finite two-dimensional array of nodes was 
selected (Fig. 53) as the framework to approximate the solution of the 
above equations by finite differences. The ceramic thickness, t c> is 
represented by seven nodes, which are equally spaced except for the edge 
pairs, where one-half the standard spacing is used. Hence, the standard 
spacing, 5 , is equal to t c /5. The standard spacing is also used vertically 
except for the end pairs. An integer, m F , labels the approximate location, 
y , of the freezing front, separating the liquid and solid silicon phases. 
I f ' m * 0 denotes the bottom row of nodes, the correspondence between m and 

y is . 
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y - (m - 1/2) 6, m ■ 1,2, ... 


(24) 


and mp Is the largest Integer contained in 1 + y^/i. 

The right-hand column of nodes is sufficient to represent the entire silicon 
thickness due to our assumed horizontal isotherms. There are no discon- 
tinuities in the two-dimensional temperature profile; hence, the silicon 
temperatures are equal to the ceramic temperatures at the silicon- ceramic 
Interface. 




• 101 
• 100 

• 99 

• 98 


K« -I 


• M f + 2 

• Mp + 1 



• M p - 1 

• 4 


• 3 

• 2 

• 1 
• 0 


H 


Figure 53. Finite-difference array for calculating SOC 
temperature profile. Right-hand column of 
nodes corresponds to silicon layer (one- 
dimensional approximation) and right-hand 
edge of substrate. The rest of the nodes 
represent locations within and at the left- 
hand edge of the substrate. 
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The details of the solution technique, a modified over-relaxation method 
21 

(Gauss-Seidel), are omitted here. Tho equations are nonlinear due to the 
Stefan-Boltzmann relations and the temperature-dependent k Q . For the latter, 
an approximate Inverse temperature dependence, k Q - 343/l' o , Is inferred 
from recommended data over the range BOOK < U () < 1685K. These nonlinear It ies 
are handled by ir^ans of coefficient updates at each Gauss-Seidel iteration. 

The steady-state temperature profile depends on the thicknesses of ceramic 
and silicon, the pull ?pw.ed, tho melt temperature, and the background 
temperatures, T and T, . Of course, T„ is fixed at 1685K. Not all of the 

O 1 “ 

foregoing can be independently specified, since the problem would be 
overdet mined and a solution would, in general, bo non-existent. As a 
mathematical convenience, we have chosen the melt temperature to be a de- 
pendent variable, and this is therefore an output of the calculation along 
with the temperature profile. Al<»o, we specify the temperature differences 
T M - T^ and T y - rather than T and T^ directly (e.g., T y - T^ - 1385K 
will result in T q being a little above room temperature because T^ will 
generally be a few degrees above 1685K. 

Interpretation -- At the present time, the complete temperature profiles 
in ceramic and silicon are of only casual interest. Isotherms obtained by 
interpolation of output data of two computer runs are shown in Fig. 54. 

The sparse, sharply-bent isotherms below the freezing front, and the 
sudden Increase in isotherm density above the point of latent heat release, 
are features found in all the runs. 

On the other hand, the silicon part of the profile is crucial and gives 

the nucleus of information for our interpretation of thickness versus growth 

16 

speed. As in the case of unsupported silicon sheet growth, a major as- 
sumption here is that stable growth of SOC requires the liquid meniscus to 
be at or above the freezing temperature; that is, heat flow in the liquid 
must not be away from the freezing front. A condition for the limiting 

stable growth rate is that the temperature gradient on the liquid side of 

18 

the interface approach zero. According to Brice, maximum silicon thickness 
at a given velocity may be obtained oy adjusting the melt temperature to a 
point just above the value where the interface becomes unstable and den- 
dritic growth develops. In a computer calculation, this would correspond 
to zero temperature gradient in the liquid at y » y F - 


21 


L. Fox, Numerical Solution of Ordinary and Partial Differential Equations , 

PergammotrTre^s~TT^S5X P^ 
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Figure 54. Isotherms for T ■ Ty. T - Ty - 5n, 

n m l, 2 by two-dimensional 

interpolation; t c =* 0.1cm, y p * 0.5cra 
T M - T q * 1385K, T M - Tj - 30OK. 

Left: t Q = 0.003cm, v m 0.153 cm/sec. 
Right: t Q =• 0.03cm, v = 0.0685 cm/ .sec 



Fig. 55 shows calculated normalized heat flux, Jt/ovLT Q as a function of 
y, for three melt temperatures at the same pull speed, where 


J - J 1 - Jy* (liquid) 


( 25 ) 


■* - - * k o 37 (solld) (26) 

In Fig. 55, the middle case (II) gives maximum stable t Q because J 1 (y F ) - 0. 
Case I, with a h.gher melt temperature and lower t Q , is within the stability 
limit, but Case .II (lower melt temper! ture) is unstable. 



Figure 55. 


armalized thermal flux in liquid and solid parts 
f the silicon layer at height y above meit surface. 

- 0.098 cm/ sec for all three cases. The degree of 
tability is determined by the melt temperature, 
ase I (stable): T M - T F - 22.30; t D - 0.0067cm; 
i (yF) > 0. Case II (stable): T M - Ty - 15.1 , 
l -5; 0100cm, Ji (yp) “ 0 . Case III (unstable). 

?. _ Tf - 6.50; t 0 " 0.0143cm; Ji (yF> < °- 


Figure 56 gives further computed data assuming the same fixed v. J 1 (y p ) is 
seen to be a raonotonic increasing function of T M> and t Q is a monotonic 
decreasing function of T M - According to this figure, silicon thicknesses up 
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to 0.01cm can be stably grown at v ■ 0.098 cm/sec on 0.1cm ceramic, assum- 
ing that heat is removed by radiation to an environment at room temperature. 



0.3 


0.2 


0 1 c 


0 5 


H -0.1 


- 0.2 


Figure 50. Normalized (y p ) and silicon thickness as functions 
of melt temperature. Pull speed fixed at 0.098 cra/sec; 
t c ■ 0. 1cm; T U - T 0 ■ 1335K; Ty - Tj. ■ 50K; yp ■ 0.5cm. 


Figure 57 shows t and J, (y_)/p^vL as functions of pull speed when melt tem- 
perature is held constant. We note that t always decreases as v increases, 
but J^ypJ/p^L reaches a minimum and then increases with v. In general, 
the model iudicates that growth is more stable at both extremes of pull 
speed (fixed Ty) and tends to become less stable or unstable at intermediate 
pull speeds. If a lower Ty were used in these calculations, there would be 
a larger range of pull speeds corresponding to unstable growth, and stable 
growth would be possible at both higher and lower pull speeds outside of 
this range. 



STABLE 


STABLE 



Figure 57. Normalized Jj (yy) and silicon thickness as 
functions of pull speed. Ty - Ty fixed at 
15.1°; t c - 0.1cm; T M - T 0 - 1385K; Ty - Ty 
- 50K: y p - 0.5cm. 

Figure 58 gives conditions of limiting stable growth; in other words, v 

and Ty are simultaneously varied to keep J^y^,) " 0. Solid curves represent 

the effects of two substrate thicknesses. Dashed lines are limiting growth 

16 

conditions for unsupported silicon from the one-dimensional anai,ais. 

2 

That analysis predicts a constant v t for each background tenrperature , 

giving straight lines with a slope of -2 on the logarithmic plot. On the 
2 

other hand, the v t product for SOC is seen to fall rapidly at higher pull 
speeds, as thermal convection in the ceramic begins to dominate the heat 
disposal problem. The experimental data shown were obtained with the "fast 
dipper" on O.lcm-thick ceramic with an undetermined background temperature. 

The Role of T^ and yp — The calculations discussed so far have been based 
on the values y p * 0.5cm and Ty - ■ 50°. To test the generality of the 

results, the calculations for the Fig. 58 curves were repeated for two 
different combinations of y p and T^t (1) y p * 0.69cm, Tj - T y ■ 50°; (2) 
y p * 0.5cm, T^ - Ty ■ 300°. It was found that these new combinations gave 


96 



the same limiting v vs. t . relationships as before. This is a surprising 
result and it indicates that neither of these parameters has a significant 
effect on the important Smiting relationship This also tend* to enhance 
the credibility of our ear’ ier analysis of the unsupporteu cu'-et silicon 
growth, in which was assumed <?qual to zero, and further details of the 
liquid meniscus, such as its shape and temperature distribution, were 
ignored. 



Figure 58. Silicon thickness versus pull speed for 
supported (1mm and 2mm substrates) and 
unsupported 1 ® growth. Maximum stable 
growth assumed LJj/yp) ■ 01 . T« - T Q ■ 

1385K. Experimental data (circles) 
obtained on lmra unslotted substrate. 

Isotherms Shapes Within the Solid Silicon — The isotherms were assumed 
flat in the silicon part of the structure. The error in this assumption may be 
estimated by analyzing heat flows up the silicon and laterally toward the 
radiation and ceramic surfaces. In the solid silicon portion, the greatest 
departure from flatness occurs near the freezing isotherm itself, so the 
discussion is restricted to that region. No attempt is made to examine the 
liquid portion. We proceed as follows. 
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Let J and J„ be flux density vectors at the corners *s shown In Fig. 59. 
in a column of cross-section dxdz, the latent heat released per unit time 

p vLdxdz . This is obviously the same quantity asjj,|dsd Z , hence, 

o 


is 



Figure 59. Analysis of freezing front shape by 
means of heat flux components. 


j - p Q vLdx/ds - p Q vLcos 


(27) 


Similarly , 


J 2 j “ P 0 vLcot a 2 


The lateral components are then 


- p vL sin a, cos a. 
J lx o 1 1 


1 „ - p vL sin x 0 cos a„ 

J 2x o 2 z 


(28) 


(29) 

(30) 


By continuity of lateral heat flow, J lx is also the heat conducted into the 
ceramic : 


j 


lx 


k c | 3T/ 3x j 


(31) 
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i 


( 32 ) 


and J 2x is the radiated heat 

J 2x * e o° ( T F 4 - V) 

- 1/2 sin -1 | 2k c |3T/3x | /d o vL | 

- 1/2 sin -1 j2e Q ojT F 4 - T o 4 ) /p q vL 


(33) 

(34) 


The angle a 2 may be evaluated directly, but for it is necessary to cal- 
culate 1 3 T / 3 x | . A three-point Interpolation formula -nay be used incorpo- 
rating the calculated temperatures at y ■ y F and the last three columns in 
the finite difference lattice, recalling the spacing as shown: 


U- a -JU-6/ 2-4 

T. T_ T r 


The interpolation formula yields 


3T 1 I 8 
TX " S T p 



at the right-hand edge. 


(35) 


Angles a 1 and a 2 are plotted against thickness, t Q , in Fig. 60 from data 
for the cold background case represented by the right-hand solid curve in 
Fig. 58. The rraximum angles occur at the largest t , which is 0.1cm. We 
restrict our attention to this condition, which is the severest test of our 
flatness assumption. 

In Fig. 61, that particular crystallization front is represented as a 
circular arc. From this geometry, the highest Ay spanned by the crysaliza- 
tion front is calculated to be 0.008cm. Since the finite difference incre- 
ment, < 5 , is 0.02cm, we conclude that the worst-case Ay is only about half 
the numerical lattice spacing, and this justifies the flatness assumption. 

It would be incorrect to conclude that most of the heat from the silicon 
is conducted into the substrate, although this is locally true near the 
interface, where > c» 2 . In fart, angle decreases and eventually be- 
comes negative with inc-easing height above the interface. In the balance, 
more heat flows from the ceramic into the silicon than the other way around. 
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Figure 60. 


Inclinations of ».«!« 1»* ^"^0^1*°°^ . 
_ _ raM/Mtm C11T 1 face . 


OL 


A . . 



Figure 61 . 


pproximate worst-case c 7® ta J;“^ tlon 
-* — / -nvnorcrerated scale). 


»«nn t orftometrv 


This is certainly to be expected in view 
the substrate to the total heat that must 


of the large contribution made by 
be radiated away. 


Part II ■ Unsupported Wedge-S h aped (Asymmetric) Growth - The growt 
supported sheet silicon or silicon ribbon from the melt can happen in 
P clPle, in either of two modes, as illustrated in Fig. 62. In the con- 
ventional mode (a), the liquid-solid Interface (LSI) is oriented appro*!- 
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raately perpendicular to the pull direction, as was the case in the supported 
growth discussed in Part I. The latent heat of fusion released at the LSI 
is carried in the pull direction and must be disposed of by radiation and 
convection to the relatively cool environment of the solid ribbon. The 
word "symmetric" fits this mode because of the plane of symmetry through the 
half-thickness points. The other mode (b) lacks a plane of symmetry, as 
diagrammed here, and may be called "asymmetric." Because the LSI is in- 
clined at a small angle, 6, to the pull direction, most of the latent heat 
of fusion is eliminated immediately by crossing transversely to the free 
solid surface. Only a small fraction of the latent heat remains to be 
disposed of downstream (upwards as represented in Fig. 62) from the solidi- 
fication wedge. It is debatable whether asymmetric growth has ever been 
successfully demonstrated, but is is critically important technology to 
pursue theoretically and experimentally. It is inherently much faster than 
symmetric growth and may have other advantages such as entrapment of 
impurities . 


x 

4 


-► y 


SOLID 




mTTB 


LIQUID 

(*) 





Figure 62. Silicon ribbon growth in (a) symmetric and (b; 

asymmetric modes. The letters v and t are pull 
speed and ribbon thickness, respectively. 


22 

A heat transfer model was developed for the wedge-shaped region in asym- 
metric growth but it does not address the problem of how to cool the ribbon 
as it travels away from this region. This report indicates that downstream 
cooling is too rapid for pure wedge-shaped asymmetric growth if a 300K 


J. Zoutendyk, "Theoretical Analysis of Heat Flow in Horizontal Ribbon 
Urowth . " 
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;imb tout is assumed for the entire solid ribbon. The cooling can be retarded 
to the proper rate by including postheaters in that apparatus. The present 
analysis applies directly to unsupported ribbon growth. Extension to sup- 
ported asymmetric sheet growth, as in SOC, will be done it a later time. 

In both symmetric and asymmetric growth it is useful to consider lines of 
demarcation between the solidification zones and the fully developed solid, 
downstream regions. These are indicated as dashed lines in Pig. 62. The heat 
flux approaching the line from below, consisting of convected heat, heat 
flux from the liquid, and any latent heat of fusion not alreudy disposed 
of, must balance the heat flux away from the line into the downstream 
region. This simple conservation principle leads to a v vs. t tradeoff in 
the case of unsupported symmetric growth, as reviewed briefly below. In the 
asymmetric growth model, the conservation of heat flux does not imply any 
direct tradeoff limitation, but it does dictate the use of postheaters. 

Review of Symmetric Growth Model — The conservation principle for sym- 
metric growth can be expressed as 

t ( pcvTp + J Q + pLv) - t ( kq pcvTp) • (36) 

where p and c are the density and heat capacity of solid silicon at T ^ T„, 

r 

Tp Is the melting temperature, the x-coraoonent of heat flux density in 
the liquid, L is the latent heat of fusion, k is the heat conductivity of 
solid silicon (independent of temperature in this analysis), and q - -dT/dx, 
the downstream component of the temperature gradient at the dotted line 
in Fig. 62a. Values used for these parameters are given in Table 15. The 
first term in Eq. (36) represents convected heat flux, the second term is 
the heat flux from the liquid silicon, and the third term is heat flux due 
to latent heat of fusion. The tkq term is heat flux by conduction into the 
downstream solid region. Note that convected heat cancels in Eq. (36). If 
J is assumed equal to zero (maximum throughput ) , the equation becomes 

pLv/k ■ q I T=T (37) 

F 

The value of q is derived Independent of the solidification zone parameters 
by integrating the heat equation for the downstream silicon portion, subject 
to appropriate initial condition. Since symmetric growth is very nearly one- 
dimensional in character, the following heat equation is convenient to use: 1(5 
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Table 15. Values for silicon parameters 
for Part 1 1 . 


Parameter 

Value 

P 

2.33 g/cm 3 

c 

0.954 J/g deg 

k 

0.216 W cm/deg 

L 

1802 J/g 

t f 

J o 

1085K 
0.0 W/cm 2 


t(k q j}$ ♦ t*SV q| - F(T , ) (38) 

where F is the cooling function depending on the ribbon temperature, T, 
and effective ambient temperature, T A . This equation is equivalent to the 
standard one-dimensional heat equation in which the coordinate x is the 
independent variable. The transformation was made by substituting -q for 
dT/dx. If heat disposal is by radiation from both surfaces to an effective 
300K blackbody ambient, Eq . (38) becomes 

t(kq + ocv q) ■ 2eo (T 4 - 300 4 ) (39) 

The initial condition is that the gradient, q, tends to zero as T approaches 
ambient; i.e., q(300) - 0. The solid silicon eraissivity is taken to be e 
* 0.46, and Eq. (39) is integrated numerically from T ■ T^ to T « Tp, since 
q(Tp) is the value required in Eq. (37). 

The q obtained from the integration depends upon t, and to a lesser extent 
upon v due to the convective term in Eq. (39). Figure 63 is a plot of q 
vs. t for various pull speeds, along with pLv/k (dashed lines). The inter- 
sections of the two sets, marked with heavy dots, reveal the v vs. t trade- 
off and illustrate the flux conservation principle. 

Ribbon Thickness in Asymmetric Growth — A thermal model of asymmetric 
growth is independent of orientation, and the pull direction has been taken 
vertically (Fig. 64) for comparison with symmetric growth. We will use a 
horizontal orientation in the later analysis. Mechanical factors could 
dictate other orientations, such as the one shown in Fig. 64. In this con- 
figuration, point A, which presumably defines the downstream limit of the 
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t (um) 


Figure 63. Heat fluxes from solidification region (dashed 
lines) and into downstream region (solid lines). 
Balance is at the intersection points (heavy 
dots), establishing the v versus t tradeoff for 
symmetric growth. 



Figure 64. Nonvertical implementation of 
asymmetric growth. 
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LSI, is fixed by meniscus stability consideration."^** The wedge angle, 0, 
is given by 

'* 3 ‘■•j | Tp“* - T^^J/pvL radian (40) 

it heat flux into the LSI from the melt is neglected.*”*" If the top— side 
ambient T is 300K, this becomes 

0 => 0 . 005/ v 

and for pull speeds of the order of 0.1 cra/soc or higher, the wedge angle 
is less than 3°. Although t) is established by thermal considerations in- 
volving the pull speed, the thickness, t, is not, since 


and the length AB is an adjustable parameter bused on the design of the 
apparatus. The position of B might be fixed, for example, by heating the 
melt surface to the left or B and allowing it to cool on the right. Equations 
(40) and (41) show that an Indirect tradeoff exists between v and t; given 
a fixed wedge length AB, t is inversely proportional to v. For the subse- 
quent calculations, this dimension is arbitrarily assumed to be 0.8cm. 

l.? H ) l >e raturu Profile in the Solidification Region -- Subject to a number of 
assumptions which will be examined, the two-diraens ionul heat equation yields 
a nearly linear temperature distribution in the solid portion of the region 
x ■ 0 to x ■ t (see Fig. 65): 

T(x.y) » T p - A • (x/t - y/t) ( 42 ) 

w i t h 


A * 



t /k 


(43) 


and 


t/e 


0 * to T* 

r 



(44) 


23 „ 

C. A. Rhodes, M. M. Sarrof and C. H. Liu, ’’Invest! 
Stability in Horizontal Crystal Ribbon Growth ” J 
00-101 (1980). 


gat ion of the Meniscus 
Crystal Growth 50, pp. 
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Figure 65. Linear temperature distribution in 
wedje-shuped solidification region. 


The validity of the assumptions rests on their self-consistency,, which we 
will attempt to establish: 

1) We assume the LSI is indeed a straight line and an Isotherm at 

temperature T„. 

r 

2) Given (1), the linear temperature distribution obviously satisfies 
the steady-state heat equation: 

k(3 2 T/3x 2 + 3 2 T/3y 2 ) - p C v 3T/3x - 0 ( 45 ' 

provided that the pcv term can be neglected, Zoutendyk 22 gives 
the following condition for neglecting the pcv term: 

vt << k/pc - 0.097 (46) 

We will return to this condition. 

3) A heat flux emerges from the LSI due to the latent heat of fusion. 

If there are no temperature gradients in the liquid, this flux is 
normal to the LSI and of constant density. That is, 

J - pvL sint) (4?) 

because the LSI is a straight line of length t/sin 0 and the total 
flux is PvLT. From Eqs . (42) and (47) we see that 
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(.48) 


J x ■ kA/1 » pvLsin^0 


and 


Jy *■ kA/t - pvL sin0 cose (49) 

4) The radiated flux density from the surface y ■ 0 is approximately 
constant, provided that A is small, and equal to eo jTp 4 - T q 4 J. 
Equating this to the vertical component J y in Eq. (49) gives 



which establishes Eq.(43)and confirms that A is small (A - 4.9° if 
t ■ 0. 1cm) . 


5) Combining Eqs . (43) and (49) also gives 


sine cose a 6 => ea |T F 4 - T o 4 |/° vL 


)/, 


(51) 


which establishes Eq. (49) . Finally , since e a t/l, Eq. (51) may be 
rewritten as 


vt = eo 1 Tp 4 - T q 4 • it/pL 


( 02 ) 


or vt = 0.005 if £ = lcra. (It has been already noted that l is 
somewhat arbitrary.) Since 0.005 << 0.097, the condition Eq. (46) 
for neglecting the pcv term is satisfied. 


Matching the Temperature Distribution to the Downstream Region -- As ob- 
served previously, the component of latent heat carried downstream is much 
smaller in asymmetric growth. In fact, since 


J x * pvL sin 9 = 


ca Tp - T q 


4 j 2 


/pvL 


from Eqs. (44) and (48), the conservation principle for the wedge model 
becomes 


/ 4 

E0 1 t f 


/ PvLk - q Tvr 


(53) 
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For a cooling environment of 300K along the whole silicon ribbon, q was 

2 4 

found to be in the 10 to 10 deg/cm range (Fig. 63), whereas the left-hand 
side of Eq. (53) is much smaller (e.g. , about 5 deg/cra fo - . v - 0. 1 cra/sec). 
It is obvious that the cold downstream environment is not compatible with 
the wedge model. We will see in more detail later (Fig. 70) what happens 
when we try to force the cold downstream environment on the wedge model. 

It is not so obvious that the cold environment precludes any kind of 
asymmetric growth. It is conceivable that some kind of curved LSI is possi- 
ble with a portion of it nearly parallel to the surface, for higher 
throughput, and another portion at a large inclination to the surface in 
order to match the high q. 

The approach taken here is an attempt to stabilize the wedge pattern by 
extending the linear distribution for some distance into the downstream re- 
gion. This will automatically satisfy continuity of not only 3T/3x but 
3T/3y as well. It is assumed that afterheaters can be constructed to pro- 
vide an arbitrary distribution, T^(x), of effective ambient temperature on 
the meniscus side of the ribbon. The other side will remain exposed to 
the cold ambient. A restriction on this approach is that T A (x) (i.e., tho 
design of the afterheaters) must depend on v and t. A change in either of 
these parameters would force a redesign of the afterheaters. The required 
T a (x) is estimated using the one-dimensional heat equation, Eq. (38), 
and then "verified" by applying it to a computer program which (numerically) 
solves the two-dimensional heat equation for the entire silicon ribbon. 

The method is demonstrated for the particular case, l ■ 0.8cm, v - 0.1 
cm/ sec, and t * 0.04cm. 


To begin, let F(T, T A ) in Eq. (38) be written 


F[T, T a (T)] - f ( T , 300 ) + f[T,T A (T)] (54) 

where f(T,U) - eo(T 4 - U 4 ) and hence F(T,300) - 2f(T,300) - 2ea 

4 4 

(T - 300 ). The procedure will be to find a T A (T) to suit our purposes and 
then convert to T A (x) using 


x 



dT/q(T) 


(55) 
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where q is obtained from Eq. (33) with T A • T^(x). Let q o (T) be the inte- 
gral of Eq. (38), with ■ 300K. This is plotted in Fig. 06 over the range 

| T ^ 1085K to T ■ 1200K. For the desired distribution, q(T), let the linear 

i temperature range (constant q * A/ 8,) be extended to x ■ 1.21 - 0.90cm, and 

let and Tg denote the half-thickness temperatures at x - 1 and x - 1.21, 
| as in Fig. 67. The values of T^ and T,, are Indicated in Fig. 66. Next, let 

q(T) increase as the temperature decreases from T^ and join it smoothly 
, to q (T) at some T„ , as shown by the dashed curve in Fig. 66. "Smoothly'' 

I means that q ■ q Q and dq/dT ■ dq Q /dT »t T ■ Tj. 
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Figure 66. Original (q Q ) and afterheater-modified (q) 
temperature gradients for the case v ■ 0.1 
cm/ sec, t ■ 0.04cm. 
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X a o »*l *=l,2t 


Figure 67. Linear temperature 
distribution. 


Some observations may be inferred from the approach as outlined so far. It 
will be noted that q and q Q satisfy 

t(k q dq/dT + pcv q) - f[t,T a (T)J (56) 


and 


t(k q o dq 0 /dT + PCV q Q ) - F(T . 300 ) (57) 

Inspection of Fig. 66 reveals that q < q Q and qdq/dT < q Q dq /dT. It follows 
that F[T,T A (T)J < F(T, 300) , and, therefore, T A (T) > 300. In°other words, 

T /\ ( T ) is realiabl e with heaters. Second, since both q, q Q and their deriva- 
tives match at T - Tg, the value of T A (T) falls to 300K at T - T 3 and re- 
mains constant at lower ribbon temperatures. Finally, to the same approxi- 
mations that were discussed earlier, continuity of both dT/dy and dT/dx is 
assured at x - l. That is, the linear temperature distribution is extended 
past x - l: (a) Since dT/dx is constant by design throughout the extension, 
dT/dy is also constant throughout the extension because d 2 T/dx 2 d 2 T/dy 2 - 
0 (neglecting the pcv term), (b) Because the ambient seen by the ribbon 
surface opposite the meniscus is 300K, it follows that 

dT/dy = (l/k)f(T p ,300) (58) 

throughout the extension, as it is also in the wedge region. This is the 
primary rationale for placing the afterheaters on the meniscus side. 

A simple function having the required properties is 

q(T> " + a i (T 2 ~ T > + a 2 (T 2 * T > 2 * where T 2 - T 2 T 3 (59) 
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with ^ and * 2 selected so that q(T 3 ) - q Q (T 3 ) *nd q'(T 3 ) “ q 0 '(T 3 ). If 
T 3 - 1300K, these are q(T 3 ) - 672.17 and q'(T 3 ) - 1027.40, and tha coeffi- 
cients are found to ba a^^ ■ 5.0)97 and a 2 • -0.0086849. Now T^(T) ia easily 
determined algebraically from Eq. (38) over the range > T > T 2 using 
q - A/£ and q' - 0, and over the range T 2 > T > T 3 using q and q' from 
Eq. (59). Then T^(x) may be found by integrating as in Eq. (55) to obtain 
the relation between x and T. 

Two-Dimensional Verification — The two-dimensional steady-state heat equa 
tion, Eq. (45), was solved numerically on a 7 x 300 grid representing a 
rectangular region t • 0.04cm wide and about 2.4cm long. The first 0.8cm 
(100 grid points) corresponds to the wedge portion, the next 0.16cm (20 grid 
points) represents the linear extension region T^ > T > T 2 , and the re- 
mainder covers the rest of the afterheater distribution, including a very 
short final pcvtion (approximately 0.04cm in length) where T^ ■ 300K. 

Boundary conditions were set up as indicated in Fig. 68. Temperatures were 
specified on segments 1, 2, and 5, whereas derivative conditions were used 
on segments 3 and 4 corresponding to the nonlinear radiation law. An ap- 
proximately linear temperature distribution was set up in the wedge region 
by treating the liquid as an extension of the solid. The temperature on the 
edge segments 1 and 2 wa t, set to increase linearly from Tj, at the LSI 
corners to T p + A at (o,t). 



Figure 68. Boundary conditions for numerical integration of 
Eq. (57). Segment 1: T ■ T F + (y/t )A. Segment 2: 

T * (Tp + A) - (x/i) A. Segment 3: |3T/3yJ ■ (1/k) 
f (T, 300). Segment 4: |3T/3y| - (1/k) f [T,T A (x)l 
Segment 5: T - 1275K. 
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The 2100-point system was solved by means of a successive over-relaxation 
method. The main features of the calculated temperature distribution are 
plotted in Fig. 69 along with the effective ambient, T A (x). The partials 
jT/Jx and 3T/3y are of particular interest and are seen to remain essen- 
tially constant through the solidification region, ending at the point 
marked with an arrow. Ir. other words, the linear profile remains intact 
throughout this region. The required T^ is approximately constant at just 
over 2000K to a distance of 1cm from the LSI-meniscus Juncture and then 
falls rapidly to 300K over the next 0.55cm. 

For comparison, Eq. (45) was solved with T A (x) sot equal to 300K, corres- 
ponding to the omission of afterheaters. The boundary conditions on seg- 
ments 1, 2, and 3 were the same as given in the caption for Fig. 68 except 
for segment 4, 3T/3y - ( l/k)f (T,300) . The results are plotted in Fig. 70. 

Note that the rapid downstream cooling causes drastic variations in both 
3T/3x and 3T/3y near toe LSI-meniscus Juncture. It is clear that tho linear 
temperature distribution is completely destroyed. 

Summary — 

Part I: Supported Symmetric Growth — For supported synanetric growth, a 
steady-state vortical growth model of dip-coated SOC yields a system of 
nonxinear, coupled differential equations for the temperature p vfiles in 
ceramic and silicon portions of the structure. Cooling of the structure 
is assumed to be by radiation to a specified environment, and tho cooling 
rate must balance the rates cf diffusion of latent heat of fusion and 
convected heat from the melt. Liquid meniscus height is fixed by hydro- 
static considerations and does not vary over the normal range of growth 
speed. Because of the complicated geometry and essential conlinearity of 
the system, numerical solution of the equations is the method of choice. 

Due to the overall heat balance requirement, only certain combinations of 
ceramic thickness (T c ), silicon coating thickness (T q ), pull speed (v), and 
melt temperature (Ty) will yield steady-state temperature profiler. In 
actual growth, the coating thickness is not arbitrarily specified but is a 
dependent function of t c> v, and T u (>Tp). 

18 

A theoretical stability criterion places a further restriction on the 
allowable combinations of growth parameters. Growth is not stable if heat 
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Figure 69. Two-dimensional temperature distribution in silicon 
ribbon drawn through the heated environment, 

T,(x) (dashed curve). T and the partials, 3T/3x 
and 3T/3y, are averaged over the width of the 
ribbon. A well-established linear temperature dis- 
tribution is indicated by the near constancy of 
3T/3x and 3T/3y up to the end of the wedge (marked 
with an arrow). 
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Figure 70. Two-dimensional temperature distribution 
in silicon ribbon with no heated environ- 
ment. The linear temperature distribution 
is totally destroyed, indicating that a 
wedge-shaped solidification zone is not 
possible in a cold environment. 
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flows backwards from the crystallization fronts; in other words, stability 
requires that dy/dy U^J(yp) < 0. For given pull speed, t Q Increases mono- 
tonically as Ty is lowered toward Tp, but d/dy U^y^) decreases at the same 
time. Maximum coating thickness for stable growth occurs, therefore, when 
d/dy U^(yp) - 0. This limiting thickness and the corresponding melt temper- 
ature both depend on pull speed, ceramic thickness, and radiation environ- 
ment . 

Finally, the one-dimensional approximation in solid-silicon portion was 
critically examined by estimating the maximum departure, Ay, from flatnoss 
for the freezing isotherm. This turned out to bo of the order of 0.01cm, 
wtuch would appear to Justify the approximation. Also, the freezing isotherm 
was found to be bowed upward, making a large inclination to the horizontal 
at the ceramic surface than at the as-grown surface. Although most of the 
latent heat is carried vertically, this would indicate that the net lateral 
heat transfer in the neighborhood of the freezing front is into the 
ceramic as opposed to out from the free surface. This is true only locally, 
however. In fact, over the structure as a whole, net heat transfer across 
the silicon-ceramic interface is from ceramic to silicon, and not the other 
way around. 

Part II: Unsupported Wedge-Shaped (Asymmetric) Growth - In part II, we 
have attempted to establish some consistent thermal conditions for high- 
speed, asymmetric growth. It is theoretically possible to grow thicker 
layers of silicon at a given pull rate in asymmetric growth. One way to ex- 
plain this is in terms of the rapid removal of the latent heat of fusion. 
This takes place on a radiating surface very close to its source in the 
liquid-solid interface. In symmetric growth, heat tends to be transported 
farther in the solid material before being released (less efficiently) at a 
lower radiation temperature. A thermal analysis of the asymmetric growth 
zone, together with the cooling portion downstream, has shown that a post- 
heater is a necessary part of the growth apparatus if a strictly wedge- 
shaped growth zone is to be maintained. 

Dip-Coating Results 

Introduct ion - When a 7cm x 5cm carbon-coated mullite substrate is dipped 
into and subsequently pulled from a silicon melt, silicon coatings or 
layers of large elongated grains result which are aligned in the direction 
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of withdrawal. When growth condition are proper, the layers are smooth and 
shiny In appearance as shown In Fig. 71. To satisfy the needs of this R&D 
program, the number dipped In this manner was in the thousands during the 
course of the Investigation. 



Figure 71. Etched surface of SOC sample. Note the 
directionality of the structure and 
grain size. 


This subsection of the report describes rho experimental options investi- 
gated, along with growth parameters of the technique and the physical and 
electrical properties of the resulting layers. 

Angle Dip-Coating - In the interest of improving grain size and/or growth 
rate, a number of substrates were dip-coated, with the substrates held at 
an angle as opposed to being perpendicular to the melt. Examination of 
th.?se silicon coatings showed that angle dipping did not drastically alter 
the rate or manner in which the coatings solidified. The coatings looked 
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quite comparable to those dipped vertically, with coating thickness de- 
creasing with increasing pull rate. 

Silicon-Ceramic Bonding - Initially, it was assumed that the silicon re- 
acted with the carbon to form a thin silicon-carbide (SiC) film which sub- 
sequently wet the silicon. It was further presumed that this thin SiC film 
penetrated the pores of the ceramic to form a mechanical bond. Occasional 
layer separation from the ceramic substrate prompted a further study of this 
silicon-substrate bonding mechanism. 

It was observed that if a very thorough and heavy carbon coating was 
applied to the substrate, thicker (125 to 250um) silicon layers sometimes 
flake off if the total time the substrate remains in the melt is short. 

This conclusion is further supported by the fact that the flaking occurs 
more predominantly near the top of the substrate (that portion which re- 
mains in the melt the least time). Visual Microscopic examination of the 
bottom surface of separated layers, along with examination of that portion 
of the substrate from which it separated, strongly suggests that the 
formation of SiC, if any, plays an insignificant role in bonding the silicon 
layers to the substrate. 

One could conclude from the visual inspection, however, that silicon bonds 
directly to the pores of the ceramic only in regions where there are 
voids in the carbon coating. These voids could occur partially by an in- 
complete carbon coating and partially by the thinner regions of the carbon 
coating entering the silicon solution upon dipping. From the beginning, 
it was observed that when the substrate was dipped approximately 1/8 inch 
into the melt and allowed to remain there for a prolonged time, this portion 
of the substrate had a poor silicon coating. This small portion of the 
silicon-coated substrate also had a distinct absence of carbon remaining. 

The possibility exists that at this high temperature there is sufficient 
oxygen present from the disintegrating quartz crucible to permit carbon 
monoxide to form, but it is also possible that the carbon is merely going 
into the silicon solution. 

Further examination of coated substrates revealed that layers which have 
separated from the substrate all remained in the melt relatively short 
periods of time. On the other hand, if the carbon coating were too thin, 
the entire carbon film could be dissolved before the silicon could 
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intermittently wet and bond to the substrate, thus producing only spotty 
wetting. While the presence of a carbon coating is essential to producing 
a silicon coating on a ceramic substrate, it appears to act merely as a 
wetting agent. 

Figure 72 shows the silicon-ceramic interface of substrate MR-8, unetched. 
The silicon-ceramic bond is apparently mechanical in nature, formed when 
the molten silicon penetrates the porous substrate. This intimate contact 
is achieved over less than 10% of the surface but is generally sufficient to 
provide excellent adherence of the film to the substrate. The black areas 
between the silicon and ceramic are believed to be holes where unreacted 
carbon was pulled out during polishing. Figure 73 shows photographs of the 



Figure 72. Silicon-ceramic interface (.unetched) of 
substrate MR-8. 


back side of a portion of a separated silicon layer and that portion of the 
substrate from which it separated. These photographs further illustrate 
that the silicon bonds directly to the substrate only in regions that are 
void of the carbon coating. 
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Cleaned Interface 


Bare Silicon 


(b) Separated layer ( 30 x ) . 

Figure 73. Substrate MR-32 and separated layer 
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Figure 74 shows photographs of s i 1 ieon-aluminu and silicon-pressed mullite 
interfaces, respectively, where the ceramics are not as porous as the rolled 
mullite and not penetrated by the molten silicon. It was noted that the 
silicon coating did not adhere as well to these denser coramics as it did 
to the moro porous mullite. However, it was also noted that if an excep- 
tionally heavy carbon coating is supplied to the mullite, the silicon does 
not adhere well to it either. 

Seeded Growt h - In an effort to promote large-grain growth, four substrates 
were dipped with single-crystal silicon seeds attached. Thin wafers of 
single-crystal silicon were simply clamped to the substrate at the top edge 
of the carbon coating. The results can be seen in Fig. 75, which shows the 
interface between the single-crystal seed and the ensuing silicon film. 
Although the desired result, single-crystal film growth, did not occur, the 
procedure shows promise. The seed-film junction is continuous, indicating 
that with a suitably oriented seed, large grain growth may be induced. 

Growth Paramolurs - it was observed early in the program tnat tee layer 
thickness was dependent on such growth parameters as melt temperature and 
growth rate (pulling speed in the dip-coating teehnigue). The layer thick- 
ness was inversely related to both melt temperature and growth rate. This 
observation, along with other unexplanable observations, gave rise to the 
previously discussed growth, modeling study in support of our experimental 

effort. Velocity thickness data are summarized in Fig. 76. The straight 

2 

lines denote constant values of v t product. We observe that at a given 
crucible position the maximum thickness at a given velocity tends to a con- 
stant v 1 2 t relationship. When the crucible is low with respect to the 
heater, the v t product also is lower. 

Early in the program we were not sure where solidification occurred along 
the SOC layer. Now, however, we believe that solidification occurs at the 
top of the meniscus. Direct evidence of this is a bright line right at 
the upper edge of the meniscus. The increased brightness is due to the 
higher emmissivity of the solid. The narrowness of the line is because of 
the steep temperature gradient in the solid. Our observations can be sum- 
marized as follows: 

1) Solidification takes place at the edge of the meniscus. This is the 

same as in Czochralski growth and ribbon growth. 
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Figure 75. Seeded silicon sheet growth. The seed is 
the trapezoidal shape at the top of the 
photo (Wright etch, 2.5x). 

2) The meniscus height is constant over the range of velocities normally 
used: 2.4 to 6.0 cm/min. 

3) The contact angle between the liquid and solid is 11 as in other 
types of growth from the melt except the r ibbon-to-r ibbon case. 

4) Starting from good conditions of growth, an increase in the melt 
temperature, T , always causes a decrease in thickness. Films can be 
grown in a melt temperature range of at least 20 to 30K. 

5) At a given velocity, T^ can be decreased to a minimum value, T mo> 
corresponding to maximum thickness. Under these conditions, a 

few dendrites form at the top of the layer, but these decay within 
a few mil’imeters growth. 

6) If T is further decreased below T , dendrites will form and will 
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Figure 76. SOC layer thickness, t , as a function of growth 
speed, v, for a number of dip-coating runs under 
different conditions. Triangles and circles are 
used to indicate series of coating runs in which 
attempts were made to hold the conditions constant. 
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not decay. The entire range of melt temperatures between dendritic 
growth to complete disappearance of dendrites is on the order of 3°K. 
The center of this range we consider to be the optimum growth tem- 
perature at that velocity. 

7) If the velocity is increased, the optimum growth temperature, T rao , 
is increased. Thus, to obtain maximum thickness at a given velocity, 
the melt temperature must be increased as the velocity is increased 
to prevent dendrite formation. 

8) Width instability has not been observed in SOC growth. This contrasts 
with the case of EFG where careful design and control are needed 

to maintain the width. 

9) SiC particle formation on the surface of SOC layers ht not been ob- 
served. At times, proturberances on the surface ".re observed, but 
these have always been found to be dendrites, formed under conditions 

where T was just below T . Again, this contrasts with the case of 
rn rao 

EFG ribbon growth, where care must be taken to avoid SiC particle 
format ion . 

The above observations, coupled with the results of our thermal analysis 

of growth, indicate that SOC growth is very ^similar to EFG, but has some 

important advantages. The thermal analysis, which applies equally to 

2 

SOC and other ribbon growth, shows that the v t relationship is very 
general. However, one restriction is that the heat current density (temper- 
ature gradient) on the liquid edge of the liquid-solid interface must be 
negligible compared with vL, the velocity times the latent heat per unit 
volume . 

The condition of nearly zero temperature gradient in the liquid is consis- 
tent with the observations listed regarding the optimum melt temperature, 
T mQ . Thus, the analyses and the experimental results lend credibility to 
each other. In fact, the value of the v 2 t product can be calculated quite 
simply, and agrees with the experimental results. 

2 

The theory basically states that the v t product is determined completely 
by the coating function. Thus, we can confidently predict that additional 
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radiative and convective cooling will enable us to increase the v t product 
by a significant factor. 

Slotted Substrate Coverage - To serve as solar cell material, sheet silicon 
conventionally must have both of its surfaces accessible for electrical 
contacts. With SOC, contact at the silicon-ceramic interface is only possi- 
ble via openings in the substrate. Initially, 1mm- diameter hole oponings 
were tried and good silicon coverage was achieved. It was observed, however, 
that there was considerable perturbation of the layer surface around each 
hole. A series of long slot openings was then placed in the substrates. In 
part of these substrates the slots were vertically oriented (parallel to 
the growth direction), in other substrates the slots were horizontally 
oriented. As discussed later in this report, in order to minimize series 
resistance in solar cells made from these layers, it was necessary that 
these slot openings cover 50% of the coated area. 

Horizontal slots were found to be more easily bridged by the silicon than 
the vertical slots. Horizontal and vertical slots that cover 50% of the 
area were successfully and routinely dip-coated. Figure 77 is a photograph 
of a vertically slotted substrate in which the slots are over 1mm wide and 
the ribs are less than 1mm wide. These layers wore grown at a speed of 3.6 
cm/rain. Figure 78 shows a substrate configuration in which the slots were 
horizontal with the same slot-to-rib ratio. We observed that there is less 
perturbation of the surface with the horizontal slots, and fo. 1 this reason 
this configuration was adopted for the remaining SOC effort. 

Back Surface Field Metallization Study - It is generally known that the 
application of a back surface field (BSF) to the back surface of a solar 
cell Improves its conversion efficiency. With SOC such application is 
difficult. Nevertheless, it was decided that it may be worthwhile to 
attempt a BSF application through the slot openings in the substrate. 

Different materials were evaluated as oack contact materials that would 
produce the BSF effect when used with the slotted SOC material. 

The first of these new materials was an aluminum paste formation. The use of 

this for BSF contacts was described in a paper by J. A. Thornhill of 
24 

Spectro Lab. We prepared two paste mixtures which are very similar to the 
24 

.1. A. Thornhill, "Automated Fabrication of BSF Silicon Solar Cells with 
Screen-Printed Wraparound Contacts," Spectro Lab, NASA Contract 
NA53-20029 , Final Report, pp. 11-15 (Aug st 1977). 
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Figure 77. Vertically slotted substrate 


(a) Front surface 


(b) Back surface 


Figure 78. Horizontally slotted substrate. The back has 
been coated with solder for a back contact 
The lower part of the front surface has been 
AH coated. 




Table 16. Single-crystal diodes with aluminum back contacts. 


vVafer 

No. 

Type of 
Aluminum 

Furnace 
Temp. (°C) 

Furnace 

Time 

Max. V tK 
(V) 

Comments 

6 

Paste 

P03 

60 sec 

0. 56 

Showed shunt and 
series resistance 

10 

Paste 

916 

40 sec 

0. 55 

Showed shunt and 
series resistance. 

11 

Paste 

~900 

25 sec 

0. 56 

Showed shunt and 
series resistance. 

12 

Sputtered 

876 

30 sec 

0. 54 

Showed shunt and 
series resistance. 

13 

Sputtered 

825 

15 min 

0. 54 

Showed series 
resistance, but the 
shunt resistance was 
much better than 
observed on the 
other samples 


Table 17. Slotted SOC diodes with sintered aluminum paste back 
contacts. 


Diode No. 

Furnace 
Temp. (°C) 

Furnace 

Time 

V oc 

(V) 

Fill 

Factor 

Comments 

1 -13A-2 

823 

5 min 

0. 515 

0. 740 


139-13A-4 

823 

5 min 

0. 507 

0. 744 


139-13A-1 

823 

5 min 

0. 321 

— 

Diode located at the 
e :ge of the mullite. 

139-13A-3 

823 

5 min 

0. 445 

•• • 

Crack In the silicon 
extending toward the 
diode. 
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in-house mixture described by Thornhi-1, one with an average particle size 
of 9um and the other an average particle size of 6ym. The second material 
was also aluminum, but was sputtered on the back surface rather than applied 
as a paste. Contact resistivity measurements and diodes were made on both 
paste and sputtered aluminum. The contact resistivity of the fired paste 
was 0.16 ohm/cm, while the sputtered aluminum was 0.010 ohm/cm. The value 
for the paste was higher than expected. This may have been due to process 
conditions . 

Diodes were made on both single-crystal and SOC silicon. The results for 
single-crystal silicon are listed in Table 16 and for the SOC diodes in 
Table 17. Aluminum samples were sintered in an oxygen/argon atmosphere 
and all the aluminum paste samples were air-baked at 210 to 250°C for 15 
minutes prior to sintering. This paste was made with Reynolds Aluminum 
1-131 with an average particle size of 6um. The improved shunt and series 
resistance on the SOC diodes in contrast to the shunting effect and high 
series resistance of the single-crystal diodes may be the result of lower- 
temperature firing. 

Tne results of this work at this point were inconclusive and it was dis- 
continued due to the loss of personnel involved. 

Fast Growth Studies - The 30C economic analysis reveals that a substantial 
cost advantage is achieved if lOOym-thick layers could be coated at a 
lineal velocity of 14 cm/min. Theoretical and experimental studies to date 
had shown that such a coating velocity was not possible with radiative 
cooling alone. For this reason a second, more versatile dip-coater was 
constructed to investigate fast growth. To overcome early problems with 
substrate breakage due to thermal shock upon cooling, an afterheater as 
shown in Fig. 9 was added to the coati r and encouraging progress was made. 

A series of runs was made with argon gas flowing through the cooling shoes 
onto unslotted substrates. In all cases, uneven, vertically striped coat- 
ings were obtained. The stripes corresponded to the gas openings in the 
cooling shoes. In some cases the coatings were partially dendritic. The 
thickest portions of the smooth regions were measured and are shown plotted 
(as crosses) as a function of growth speed in Fig. 79. 

The line drawn on the far left side in Fig. 79 represents a large amount 

of data taken on dip-coater No. 1 with radiativ' cooling only. The data 

2 9 

follow the relationship v ' t = 0.059. The circles are data taken in 
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Figure 79. Summary of data obtained with dip-coater No. 2 
using an "active" afterheater and helium-gas 
cooling. The open circles are taken on substrates 
coated without gas cooling, and the line is a 
summary of data taken on the other dip coater under 
the same experimental conditions. 


dip-coater No. 2 with radiative cooling only. The agreement between the two 
coaters is good. The crosses are the data taken with argon flowing through 
the cooling shoes. The argon data lie on a line roughly parallel to the 
other and is displaced toward higher v and t. With argon cooling, a smooth 
coating thickness of 300pm was achieved at a growth speed of 6.0 cm/min, 
but all coatings produced at 12 cm/min were less than 100pm thick. 


A series of runs was then made with helium flowing through the cooling 
shoes in an effort to got thicker coatings (the thermal conductivity of 
helium is approximately eight times that of argon). These results are shown 
as the triangles plotted in Fig. 79. The appearance of the coatings pro- 
duced under helium was the same as observed previously (i.e., vertical 
striations corresponding to individual gas jets plus some dendritic 
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regions). The major difference is that the coatings grown under helium are 
thicker. Coatings exhibiting smooth regions with thicknesses of 200 and 
300um have been grown at 12 cm/min. An example of a sample with a 200i.m- 
thick coating produced on a slotted substrate at 12 cm/min is shown in 
Fig. 8o The smooth portion of the coating has thicknesses ranging from 
25 to 200um. Smooth coatings were also grown at 14 and 18 cm/min under 
helium and the data are plotted in Fig. 79. These data must be considered 
preliminary because the melt temperatures have not been optimized for the 
higher-speed growth. As indicated in Fig. 79, a goal of lOOum thick at 
a growth speed of 14 cm/min appears to be feasible. 



(a) Front Side 


(b) Back Side 



Figure 80. Photograph of a slotted substrate coated at 12 cm/ 
rain w)th helium gas directed at (a) the front and 
(b) the back sides of the substrate. The thicker 
regions are 200um thick. 
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Physical Propart lea of Dip-Coated Material - When ceramic substrates are 
carbon coated on one side and then dipped Into molten silicon, directional 
solidification occurs upon withdrawal on the carbon-coated side. With di- 
rectional solidification, the primary crystallization originates from the 
liquid-solid interface as opposed to proceeding from the surface of the 
substrate. Thus, continued nucleation results from previously grown silicon, 
thereby permitting single crystallites, after a few millimeters of growth, 
to enlarge in width as well as in length. The films are reasonably uniform 
over the area and their thickness, as previously mentioned, can be con- 
trolled by varying the melt temperature and withdrawal rate. Through the 
course of this investigation, layer thicknesses from approximately 10um 
to 300pm have been grown. 

Since the final selection of mullite as a substrate material was made on 
the basis of thermal expansion coefficient and thermal shock, little, if 
any, cracking occurs in the substrates or silicon layers. 

The crystallographic texture of these SOC layers was determined by x-ray 

nc i n 

diffraction. The structure is quite similar to that of ribbon. ’ The 
predominant structure is heavily twinned with jlll( twin planes perpendi- 
cular to the surface and parallel to the <211> growth direction; the sur- 
faces are close to the {lio} planes. A photomicrograph of an etched surface 
of an SOC sample showing twin boundaries is shown in Fig. 81. Note that 
some areas between twin boundaries are populated with dislocations and 

others are not. The grains are heavily twinned and the dislocation densities 

(3 g o 27 

are high, 6 x 10 to 8 x 10 cm - . According to some theories, the 

twinned structure is necessary to permit rapid growth. 


25 

J. D. Zook, R. B. Maciolek and J. D. Heaps, "Silicon-on-Ceramic for Low- 
Cost Solar Cells," IEDM Tech. Digest , IEEE Press (1977), p. 218. 

26 L. C. Garone, C. V. Hari Rao, A. D. Morrison, 1 Surek and K. V. Ravi, 
"Orientation Dependence of Defect Structure in £FG Silicon Ribbons," Appl. 
Phys. Lett . 29 511 (1976). 

27 

M. H. Leipold and R. J. DeAngelis, Proc. of Photovoltaic Solar Energy 
Conf. , Luxembourg, (1977) p. 872. 
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Figure 81. Photomicrograph, of an etched surface of an 
30C sample showing twin boundaries and 
dislocation distribution. 


A photomicrograph of a metallographic cross section taken perpendicular to 
the growth direction is shown in Fig. 82. The twin boundary and dislocation 
distributions coincide with that observed on the surface. Note that the 
layer is only one crystal thick and that there is a contiguous interface 
between the silicon and the raullite substrate. 

In the seeded growth experiments, where single-crystal and EFG ribbon 
material were used to enhance grain growth, it was determined that at the 
accelerated growth rates of interest in this investigation, seed orientation 
was a critical factor. Successful seeding was achieved only with seeds 
having the same orientation as the texture of the rapidly grown films, i.e., 
J 110 | <211> . A photograph of an SOC film grown at 2.4 cm/min and seeded 
using single-crystal silicon ribbon material |110}<211> is shown in Fig. 75. 
The surface was etched to reveal how the structure of the seed extends into 
the SCC layer. The use of other seed orientations resulted in the growth 
being influenced for only a few millimeters before the interface became 
unstable and nucleation took place, with growth proceeding as if it had not 
been seeded, i.e., {110} <211>. 
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millimeters before the Interface became unstable and nucleatlon took place, 
with growth proceeding as if it has not been seeded, i.e., |110| <211>. 



Figure 82. Photomicrograph of a raetallographic cross 

section of an SOC sample taken perpendicular 
to the growth direction. Note that the layer 
is only one crystal thick and that a con- 
tinuous interface exists between the silicon 
and ceramic. 


Electrical Properties of Dip-Coating Material - Early in the program, the 

15 -3 

silicon melt was doped with boron to a level of 5.65 x 10 atoms cm , 
corresponding to a resistivity of 2.5 ohms-cm based on a single-crystal 
mobility. Invariably, however, dip-coated layers ended up with resistivities 
much lower than 2.5 ohms-cm. It was suspected that the molten silicon was 
dissolving the mullite substrate, thereby freeing aluminum atoms to contri- 
bute to the p-type doping level. 

The primary means for determining impurity doping levels was using a four- 
point resistivity probe. The four-point probe technique gives an accurate 
measure of resistivity if the thickness of the layer is known, but the 
carrier concentration is obtained by assuming that the mobility in the poly- 
crystalline layers is the same as j.n single-crystal silicon. This assump- 
tion, in general, is not valid; therefore, three Hall samples were prepared 
from two dip-coated substrates (MR-1 and MR-7). Of the three samples 
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prepared from dip-coated layers, two were prepared from substrate MR-7 
which had the melt surface solidify around it. This required it to soak in 
the melt for 9 minutes before remelting was accomplished and temperature 
equilibrium could be established. In a sense, the samples prepared from 
MR-7 thus represent a worst-case example as far as time allowed for the 
substrate to contaminate the silicon coating. The remaining substrate was 
prepared from MR-1, which represents a more typical time in which the 
substr'.le is in contact with the molten silicon. 

Measurements were made on these three dipped layers to determine if the 
resistivity measurements could be used to determine the carrier con- 
cent rat ion . 

Hall samples were prepared by masking the sample configuration and sand- 
blasting the remaining silicon from the substrate. The samples were about 
lcm long and 1 to 2rara wide, with the voltage legs about 5mm apart. Contacts 
were made by evaporating aluminum and alloying at 750°C for 10 minutes. 
Electroless nickel was plated over the alloyed aluminum and indium 
soldered to the nickel. The measurements were made in a standard manner 
using a d-c field to 2 kilogauss. 

The room temperature measurements are summarized in Table 18. The sample 
designations MR-7-P and MR-7-X refer to the crystalline orientation in the 
Hall sample, with "P" meaning the long axis of the grains were parallel to 
the direction of current, and "X" meaning the long axis of the grains were 
oriented across the sample, perpendicular to the direction of current. As 


Table 18. Sumr ary oi 300K electrical properties. 


Sample 

No. 

Carrier Cone. 
(cm" 3 ) 

Resistivity 
(ohm- cm) 

Mobility 
(cm 3 /V sec) 

MR-7- P 

6. 7 x 10 16 

0.45 

208 

MR-7-X 

7 x 10 16 

0. 6 

146 

MR-1 

7. 5 x 10 16 

0. 54 

153 

Bulk Poly 
(Melt- 8) 

1. 8 x 10 17 

0.24 

146 
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seen in Table 18, these two samples have slightly different nobilities, 

with the highest mobility being in sample MR-7-P. The mobilities in all 

o 

three samples are lower than the 250 to 300 cm /V-sec expected for single- 

17 3 

crystal material with carrier concentrations in the 0.6 to 2 x 10 /cm 
range. 

The temperature dependencies of carrier concentration, resistivity, and 

mobility for sample MR-7-X are shown in Figs. 83 and 84. In general, they 

behave as expected for silicon: the carrier concentration decreases and the 

resistivity goes through a minimum before increasing as the sample is 

cooled. The one exception occurs at 77K where the carrier concentration 

apparently increases over the 85K value. The temperature dependence of the 

mobility is compared with a single-crystal sample with a carrier concen- 
17 3 

tration of 2 x 10 /cm in Fig. 84. The notable difference is the sudden 
drop in mobility in the layers at 77K. No explanation for this decrease is 
postulated at this time. 

The main point to note from these measurements is that the mobility in 
polycrystalline samples will be lower than in a single-crystal sample with 
the same doping density. As a result, the resistivity measurements will 
indicate a higher purity by as much as a factor of 2. 

To investigate the source of this excess p-type conductivity in melt 8, a 
thin wafer was prepared from this melt and an infrared transmission mea- 
surement made at a temperature of 10K using an FTS-14 Fourier Transform 

Spectrometer. The absorption coefficient is shown in Fig. 85 as a function 

-1 4 -1 

of photon energy in cm X(pm) * 1 x 10 /K(cm ) . The prominent absorption 
lines are labelled according to the standard nomenclature for acceptors in 
silicon. At the lower energies, three lines due to boron are clearly seen 
and are labelled 1 through 3. At the higher energies, a number of lines 
characteristic of aluminum are seen, labelled XI, X3, 1, 2, and 4. The most 
prominent aluminum line (line 2) is totally absorbing in this particular 
sample because of the relatively thick sample used at this high carrier 
concentration . 

An estimate of the impurity concentrations of both boron and aluminum can 
be made from the intensity of the absorption peaks. These are (4 ± 1) x 
10 15 /cra 3 for the boron and (8 t 4) x 10 16 /cm 3 for the aluminum. The alumi- 
num concentration is more uncertain than our estimate of the boron concen- 
trations, since we did not run a standard with aluminum doping. We had to 
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Figure 83. Carrier concentration and resistivity vs. 1000/T 
(deg K) - * for sample MR-7-X. 



Figure 84. Mobility vs. temperature (K) for sample !'R-7-X. 


i 
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rely on aluminum spectra which have appeared in the literature. The boron 
concentration in the sample is approximately the amount loaded in the melt, 
but less than the aluminum. 



Figure 85. Absorbance spectrum of a polycrystalline silicon 
sample. 


The first and most obvious conclusion to be made from the above measurement 
is that raullite was, in fact, being dissolved by the molten silicon. This 
was probably true, but it would have been a mistake not to consider the 
possibility of some coincidental source caused by uncleanness in the sub- 
strate processing technique or the impurities in the carbon coating. Tylan 
Corporation of Torrance, California, was then requisitioned to coat all 
surfaces of several raullite substrates with their reportedly impervious 
glassy carbon coating. Discussions with Tylan were encouraging. Based on 
their experience with coating impure graphite crucibles, they felt that a 
Tylan coat in t would remain intact and that molten silicon would '.vet it. 
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The thermal expansion coefficient of mullite was compatible with Tylan's 
coating and the scale-up cost of the coating process should be modest. 

Upon receipt of the "Vitregraph" coated substrates, the extent to which 
the mullite dissolution could be inhibited was investigated. 

The Investigative experiments are summarized in Fig. 86, which gives re- 
sistivity as a function of time in the melt. Samples were prepared from the 
tops of solidified melts and resistivity measurements were made, where 
possible, within single grains of silicon. The top of the solidified melt 
was selected to be representative of the solidification rates involved in 
actual dip-coated layers of silicon and thus should more accurately simulate 
the effects of impurity segregation. 

An undoped silicon charge was melted and a 125um-thick silicon layer was 
dip-coated on a mullite substrate (MR-45). Following this, an uncarbonized 
substrate was dipped and soaked in the melt for a period of 63 minutes. 

This shaking procedure was intended to contaminate the melt. Following this, 
a third substrate (MR-47) was then also dip-coated with a 125ym-thick layer 
for the purpose of comparing its resistivity with that of layer MR-45. 

The respective resistivity levels of these two layers along with that of 
the melt are shown in Fig. 86. The above experiment with a new silicon 
charge was repeated with substrates which were fully coated with graphite 
on all sides and edges. The first dipped substrate (MR-48) was in the melt 
for 22 minutes. As shown in Fig. 86, the rubbed-on graphite coating neither 
helps nor hinders the final impurity level of the melt. 

Using another undoped silicon charge, five mullite substrates (MR- 101 
through MR- 105) were progressively dip-coated with silicon layers using 
identical growth conditions. Their respective resistivities, as the melt 
became progressively contaminated, are also shown in Fig. 86. 

Fig. 86 shows that the aluminum impurity level of the melt increases 
rapidly during the first 2 or 3 minutes of contact time with the substrate, 
and then increases at a lesser rate, becoming very nearly saturated after 
an hour of contact time. 

The initial rate of dissolution would be proportioned to the amount of 
exposed surface area, with the rate decreasing as the solution approached 
saturat ion . 
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TIME (MIN) 

Figure 86. Layer and melt resistivities V3. accumulated 
time that the undoped melts were exposed to 
mullite. 
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The effective surface area of the mullite depends on the surface roughness 
and on the degree of wetting of the substrate. The saturation solubility of 
the aluminum Is determined by the amount of oxygen In the melt. Because of 
the quartz crucible, the silicon Is saturated with oxygen, limiting the 
amount of dissolved aluminum to a concentration well below the solubility 
limit of aluminum itself. Since several chemical reactions are involved, 
the detailed kinetics are undoubtedly not simple, as indicated by the 
several slopes in different regions of Fig. 86. 

Figure 86 also shows the effect of the vitreous carbon ( "Vi t regraph" ) in 
inhibiting the mullite dissolution rate in molten silicon. Vitreous 
carbon-coated substrates MR-97-TC and MR-98-TC were successively dip-coated 
in an undoped silicon melt. Substrate MR-98-TC was allowed to soak in the 
melt for 1 hour before it was removed with a smooth, uniform silicon layer. 
MR-97-TC and MR-98-TC layer resistivities are also plotted in Fig. 86 for 
the purpose of comparing them with the more generally used rubbed-on , 
graphite-coated substrate silicon layers. 

Sectioning of layer MR-78-TC caused portions of the vitreous carbon 
silicon layer to separate from the substrate. This made it possible to 
measure the sheet resistance of the silicon layer as well as that of the 
attached vitreous carbon coating. Since the sheet resistance of the silicon 
layer was approximately 20 times that of the vitreous carbon, it is 
evident that some type of conductivity barrier exists between the two. 

The sheet resistance of vitreous carbon on the uncoated portion of the 
dip-coated layer was very similar to that attached to the silicon layer, 
giving further evidence that little, if any, of the vitreous carbon coating 
was dissolved by the molten silicon after soaking in it for a period of 
1 hour. 

Evaluation of SQC Material Using LBIC - During the course of this program, 
a technique was developed to evaluate in detail the various contributions 
to the short-circuit current of SOC solar cells. This technique, called 
LBIC ( light -beam- induced current) technique, resembles other light-beam 
scanning techniques but differs somewhat in the details. The technique a..d 
results have been described in various SOC Quarterly and Annual Reports and 
in several publications. Therefore, the purpose of this subsection is only 
to summarize the highlights. A more complete description is given in Appen- 
dix A, which is based on a 1980 Photovoltaic Solar Energy Conference at 
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Cannes, France. Summaries of LBIC results also were given in SQC Annual Re- 
ports Nos. 2, 4, and 5. The basic features of the LBIC experimental tech- 
nique used at Honeywell are as follows: 


e A high-pressure xenon arc lamp is used as a source. It provides peaks 
at convenient wavelengths in the infrared with a focused spot in- 
tensity of about 1 sun. 

e The focused spot can be a few micrometers in size, typically 3 x Hum, 
measured on the sample. 

e The absolute quantum efficiency, Q, is obtained from the ratio of the 
LBIC current in the sample to the LBIC current in a calibrated, highly 
uniform silicon photodiode inserted in place of the sample. 

e The minority carrier diffueion length, L, is obtained from the abso- 
lute quantum efficiency, Q, by solving 




i 


1 - R) alL 
+ °X L 


where a is the absorption coefficient at wavelength X, 


(60) 


e By averaging the values of L obtained for six different wavelengths, 
an average value of L is obtained. T^e value oi (1 - R) can be 
selected to minimize the mean square deviation, AL/L. 

e At the longest wavelengths (approximately l.Oum), the values of aL 

tend to be much less than unity, so that the LBIC scan gives a mapping 
of an effective diffusion length as a function of position, as seen 
by the above equation. 


e LBIC can also be used to evaluate material before processing, by using 
an electrolyte technique. The signal-to-noise ratio is not as good 
as in the d-c case, however, since chopped light must be used. 


The theory of LBIC response has also been developed, assuming that an ef- 
fective recombination velocity parameter, S, can be used to describe grain 
boundary recombination. The analytical results were published in Appl. Phys 
Lejt_t. 3_7, 223 (1980), as well as in SOC Annual Reports Nos. 4 and 5. Some 
comparisons with the theory are given in Appendix A, as well as in these 
A*"iual Reports. 
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The experimental results of LBIC measurements on a number of SOC samples can 
be summarized as follows: 

• There usually is a characteristic maximum diffusion length, L . This 

is evident from the good LBIC response at 0.98ym, and generally occurs 
in the largest visible grains (greater than 400um wide). Cross sec- 
tions of samples at lines where the LBIC scan is known show that the 
best diffusion lengths occur in grains having high densities of twin 
boundaries and few dislocations. 

• There usually is a characteristic minimum diffusion length, L min . 
characteristic of the most active grain boundaries. Surprisingly, the 
depth of the minimum LBIC response at 0.98yra does not depend strongly 
on the spot size or on obtaining exact focus of the light beam. The 
minima of the LBIC response tend to be constant, even if the substrate 
is bowed, so that part of the scan is out of focus. Varying the spot 
size from 5 to 15ym does not significantly affect the value of the 
minima. Thus, a "minimum diffusion length" is a fairly well-defined 
experimental quality for a given SOC cell. 

• Different large grains within the same sample may have different 
diffusion lengths. Even adjacent grains may have values of L differ- 
ing by 25%. 

• SOC material appears to consist of large-grain regions and small- 
grain regions. The small-grain regions may have a significantly 
smaller diffusion length than the large-grain regions, especially in 
processed cells. An example is shown by the LBIC scan in Fig. 87. 

• An average diffusion length, L, can be obtained from a single LBIC 
line scan by electronically integrating and averaging the LBIC 
signal during a scan. 

• The values of L max i L m i n > and L were found to depend significantly on 
the boron doping in the sample, with a strong decrease in L v^alue 
with doping, at least using the "fast cool" processing sequence. The 
results are shown in Fig. 88. 

• The value of L is also found to be a function jf processing. Using 
the electrolyte technique before cell processing and comparing the 
results with the p-n junction results after thermal diffusion using 
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Figure 87. Correlation of LBIC scan with structure. A solar cell was scanred and then 
cross-sectioned through the scan line. The LBIC scan (upper trace) is at 
the same magnification as the cross-section (middle photo) and the front 
surface view (bottom photo). The letters, such as LI and V, correspond to 
minima in the trace. The numbers correspond to maxima. The distance between 
the 10.0 and 8.0 markers is 2.0mni. I. the cross-sectional view, the ceramic 
is on top and the top surface below it is rougher because of inevitable 
chipping during polishing. 
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a phosphane gas source the value of L seems to be Increased by the 
diffusion. It is further increased by the "slow cool" post-diffusion, 
as shown in Fig. 88. The response within grains as well as the re- 
sponse at grain boundaries is improved. 



BOHON DOPING, N A (cm 3 ) 


Figure 88. Diffusion length as a function of doping and cell 
processing. 


• Hydrogen passivation performed at Sandia Labs is effective in reduc- 
ing grain boundary recombination. The details are given in Appendix A. 


• The effective diffusion length, or minority-carrier lifetime, is not a 
function of intensity, either within grains or at grain boundaries. 

The experiment is described in SOC Annual Report No. 5. Nevertheless, 
there are some discrepancies in the values of L measured on two 


Slow cooldown, following the phorphorus diffusion 
Solar Cell Effort portion of this report. 


is discussed in the 
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different systems. The second system illuminates the whole cell, but 
with much lower intensity. The latter technique gives somewhat higher 
values (50% higher) of L than the average scan LBIC technique. The 
discrepancy remains to be resolved. 

• There are some grains in which the value of L is almost as low as at 
the grain boundaries. These regions are largely responsible for 
lowering the average value of L below the value in the good grains. 

By correlating a metallographic section with the LBIC scan shown in 
Fig. 87, these regions can be seen to be high densities of disloca- 
tions. Thus, the cause of low L appears to be structural rather than 
chemical. The largest values of L occur in large grains that are 
heavily twinned and have few dislocations. 

In summary, the LBIC technique has been found to be quite useful in evalu- 
ating minority-carrier lifetime effects in SOC material. It is found that 
the present limitations in efficiency are due to both grain boundaries and 
defects within grains. The indications are that structural defects (grain 
boundaries and dislocations) are more important in limiting the present 
performance of SOC cells than are chemical defects within grains. On this 
basis we expect that SCIM-coated material with large, heavily twinned 
grains will produce improved efficiency. Hydrogen passivation, by reducing 
the electrical activity of grain boundaries and dislocations, appears to 
be another way to improve efficiency. 

SCIM-Coating Results 

Introduction - The first problem encountered when SCIM-I became operational 
was related to substrate breakage. As detailed in the Sheet Silicon Coating 
Facilities subsection of this report, this problem was corrected by adding 
appropriate pre- and post-substrate heaters and adjusting the system to 
achieve a proper longitudinal temperature profile and an adequate trans- 
verse temperature gradient along the entire substrate path. 

A second and more time-consuming problem at the outset of SCIM-coating was 
solving a meniscus instability problem which turned out to be related to 
horizontal coating. This problem is illustrated in Fig. 89. As the sub- 
strate passed over the trough, molten silicon wet the substrate and was 
drawn downstream of the trough before solidification occurred (see t^). As 
the substrate continued its motion, this molten layer of silicon thickened 
(see t 2 through t g ) until it eventually spilled over and annihilated the 
trough heating element located below. 
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Figure 89. Illustration of molten spillover problem 
in relation to coating length. 


For the most part, the solution of this problem resulted from a Honeywell- 

sponsored effort along with support from SERI Contract No. XS-9-8119-1 and 

28 

was reported in the January, 1981, Final Report to SERI. The solution was 
relatively simple. Meniscus stability was attained by tilting the coater 
to an angle greater than 10°. Additional modeling was done on growth stabil- 
ity and is discussed later in this subsection. 

The design of pre- and post-substrate heaters in SCIM-II borrows heavily on 
that which was finally successfully used in SCIM-I. The pre- and postheaters 
in SCIM-II proved to be adequate for fast (greater than 10-cm/min) coating 


J. D. Heaps, S. B. Schuldt , B. L. Grung and J. D. Zook, 'Supported Growth 
of Polycrystalline Silicon Sheet on Low-Cost Ceramic, Carbon, or Reusable 
Substrates," Final Report, SERI Contract XS-9-8119-1 (January 1981). 
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speeds. The thermal design analysis played an important role in achieving 
the proper longitudinal profiles needed to prevent substrate warpage and 
breakage. The actual modifications made to improve this profile are des- 
cribed in the Sheet Silicon Coating Facilities subsection. 

SCIM-Coated Silicon Layers - One of the better SCIM-coated layers to come 
out of SCIM-I is shown in relation to dip-coated samples in Fig. 90. The 
coating was made on an unslotted raullite ( K-roodi f icat ion ) substrate which 
was 5cm wide by 55cm long and 0.1cm thick. The layer, as with dip-coated 
layers, was relatively smooth and shiny, with very long columnar grains up 
to 5mm in width. The layer was grown at a coating angle of 20° and a coating 
rate of 3.6 cm/mln. Its nominal thickness was lOOura. 






Figure 90. Photograph of SOl sheet silicon showing dip-coated 
and SCIM-coated substrates. 


Later in the program, larger layers were coated onto 10cm x 100cm slotted 
substrates using SCIM-I I. A photograph of one of these coated substrates is 
shown in Fig. 91. This layer was also grown at a coating angle of 20° but at 
a lower coating rate of 3 cm/min and its nominal thickness was greater, 
being approximately 150ura. 

SCIM-coated layers, in general, have better developed columnar single- 
crystal grains which are somewhat wider and significantly longer than 
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dip-coated layers. This is understandable, since the shorter length of dip- 
coated substrates terminates the development of many grains which probably 
would have otherwise been similar to SCIM-coatings . 



[4 — ■ 100cm 


Figure 91. Photograph of a 10cm x 100cm SCIM- 
coated SOC layer. 


SCIM-Coating Growth Parameters - As in dip-coating with radiative heat of 

fusion removal only, SCIM-coated layers also approximately obey the constant 
2 

v t relationship. Thus, if faster coating rates are to be achieved, con- 

28 

vective cooling or some form of asymmetric cooling will have to be 
implemented . 

Figure 92 sketches the coating menisci. It shows that the trough mensicus, 
as well as the inverted meniscus, is controlled by the displacer. The in- 
verted meniscus is also affected by tilt angle, trough meniscus, and the 
thermal conditions which position the liquid-solid interface (LSI). 

When the tilt angle is 20° and the trough meniscus is horizontal, the in- 
verted meniscus comes up from the trough with its center of curvature on 
the exit tunnel side of the trough and the LSI is above point A in Fig. 92. 

At this same tilt angle, when the trough meniscus is pushed up above the 
trough, the inverted meniscus goes through a flip-flop condition and be- 
comes S-shaped. The lower portion of the S hangs out over the outside of the 
trough and threatens to spill onto the heater. Under these threatening 
conditions the LSI moves out 1 to 2mm so that it is above point B in 
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Fig. 92. This condition of instability of the inverted meniscus is in a very 
narrou range of trough meniscus height, whereas the stable (concave) in- 
verted meniscus does not apparently change over a wide range of trough 
meniscus positions. Menisci observed within the SCIM coater at 10, 15, 
and 20° tilt angles are in good agreement with the calculated shapes to 
be discussed. 


TROUGH 

MENISCUS 


INVERTED 


20' TILT 



Figure 92. Meniscus shape during coating. The trough 
meniscus produces a liquid-solid inter- 
face at the lettered position. The LSI 
at "A" and "B" is <2mm apart. "B" can 
grow layers 100pm thicker than "A" but 
the inverted meniscus at "B" is unstable 
at 20° tilt. 


A coating problem which has been consistently experienced with the SCIM-II 
coater is related to a dendritic growth which occurs near the center of the 
substrate and continues longitudinally along the substrate. The cause of 
this dendritic growth is not understood, but when an argon gas flow is 
blown on the inverted meniscus, the dendrite disappears. Figure 93 is a 
photograph of a 10crr. x 100cm silicon coating where argon was periodically 
administered to the meniscus. As can be readily observed, whenever the gas 
flow was shut off, the dendritic growth reappeared. This dendritic growth 
was not, however, characteristic of layers coated in the SCIM- I coater. In 
consideration of this fact, existing differences in the two coaters were 
reviewed and found to be the following: 

1) The superior continuous coating features of SCIM-II cause it to be 
less gas-tight than SCIM-I. 


149 


ORIGINAL PAGE IS 
OF POOR QUALITY 




Figure 93 . Silicon coating where argon gas is periodically 
administered to the meniscus. 


2) The trough in SCIM-I remains horizontal irrespective of the coating 
angle used. SCIM-II initially was designed witn a horizontal trough, 
but was later tilted 10 ° toward the path of the substrate to enhance 
meniscus stability. 

3) The trough heater in SCIM-I is powered by a filtered (low-ripple) d-c 
power supply. The supply powering the SCIM-II trough heater, on the 
other hand, is merely SCR pulsed d-c and, because of large field 
variations, it imparts mechanical vibrations to the silicon melt. 

4) In SCIW-I the substrate preheater extends farther downstream over 

the trough than is the case in the SCIM-II coater. Here again, SCIM-II 
was initially constructed similarly, but was later modified in pur- 
suit of a more optimum longitudinal thermal profile. 

SCIM-II was made more gas-tight by adding gas Jets to entrance and exic 
gates through which the substrates pass. This modification makes it diffi- 
cult for air to enter the coater against the outward flow of these Jets. 

Our concern with excess oxygen entering the coater is based on our belief 
that it generates carbon monoxide which interacts with the melt to produce 
dendrites because of constitution supercooling. The addition of these gas 
Jets did not, however, improve growth conditions responsible for this 
dendritic growth. 

The orientation of the trough with respect to the coating angle was then 
altered, but dendritic growth continued. Lastly, near the end of the pro- 
gram, the substrate heater was again extended further downstream. In the 
runs that followed, considerable improvement in coating uniformity was 
experienced . 
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On the basis of a few runs only, it would appear that the longitudinal 
location of the substrate preheater is at least one of the sensitive para- 
meters contributing to dendritic growth. 

Growth Stability - During the course of this work we have learned much c )out 

the conditions for stability in supported growth. The reasons for stability 

in the dip-coating process are discussed in publication titled "Supported 

19 

Growth of Sheet Silicon from the Molt.” The early success of the project 
was due to the high stability of the meniscus during dip-coating. Because 
the liquid-solid contact angle is essentially the same for all (non- 
dendritic) growth conditions, the height of tne meniscus is constant. In 
SCIM-coating, the situation is more complex. 

The purpose of this subsection is to present some calculated meniscus 
shapes that agree quite well with the experimental observations, and to 
give theoretical arguments that provide a basis for selecting the coating 
angle of 20° as a preferred condition for SCIM-coating. 

The basic SCIU-coating configuration currently used is shown in Fig. 94. 

In this figure the substrate is inclined at 20° and the substrate thickness, 
trough size, and meniscus shape are shown to scale. In contrast to dip- 
coating, where only one meniscus shape is possible, in SCIU-coating an in- 
finite number of meniscus shapes art possible. The method of calculating 

meniscus shapes and the basic conditions for meniscus stability are des- 

28 

cribed in the Final Report to SERI, on Contract No. XS-9-8119-1. Some 
calculated meniscus shapes are shown in Fig. 95. In performing the calcula- 
tions, the meniscus pressure and the contact angle at the trough were used 
as parameters, but in actual experiment, the substrate angle and height 
are the parameters, and the contact angle can vary, as shown in Fig. 94. 

Figure 94 basically illustrates desirable conditions for uniform thickness. 
The condition is that the locus of the LSI be perpendicular to the sub- 
strate. Then, if the substrate height varies (due to nonflatness, for 
example), the position of the LSI with respect to the substrate will not 
vary. Thus, the rate of heat removal, and consequently the thickness, will 
be constant. 

There are also other conditions that are desirable for SCIM-coating. For 
example, it would be desirable to have a short meniscus to minimize edge 


LOCUS Of 
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Figure 94. SCIM-coating configuration shown with calculated men- 
iscus shapes for a pressure parameter, P - 0.5. The 
substrate is moving from left to right at an angle of 
20 to horizontal. The position of the liquid-solid 
interface is determined by the point at which the slope 
of the meniscus- is 20° + 11°, where 11° is the co.it act 
angle . 



Figure 95. Calculated meniscus shapes and loci of constant slope. 

The three cases show three meniscus pressures: (a) 

P ■ 0.5; (b) P * 0; (c) P * -0.5. The curves labelled 
from 0° to 30 are loci of constant meniscus slope. The 
distance scales are dimensionless, with capillarity 
length as the scaling factor (5.4mra for silicon). The 
pressure, P, is the melt reservoir level with respect to 
the lip of the trough. 
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effects due to curvature of the meniscus at the edge of the substrate. From 
Fig. 95 it is clear that '.he t rtest menisci are obtained by operating at 
negative pressure. A complete schematic with negative pressure (i.e., 

P - -0.5) is shown in Fig. 96. It is clear that a disadvantage of the nega- 
tive pressure is that the trough meniscus is too low to ensure positive 
contact at the back meniscus. Figure 94, on the other hand, shows that with 
sufficient positive pressure, the trough meniscus is higher than the sub- 
strate, thus ensuring positive contact at the back meniscus. 



Figure 96. SCIM-coating configuration shown with calculated 
meniscus shapes for a pressure parameter P ■ -0.5. 

For this case, the trough meniscus is too low 
to contact the substrate. Thus, to achieve initial 
wetting, the pressux-e must be increased. It can then 
be decreased to achieve negative pressure. 

From the above considerations, it is clear that the most stable growth con- 
dition is that in which the loci of constant meniscus slope have the least 
curvature, hence the widest range of stable substrate heights. This condi- 
tion occurs for the case of P * +0.5 with slope * 30° in Fig. 95. This is 
the special case shown in Fig. 94, corresponding to a substrate slope of 
30° - 11° - 19°. Thus, this is a condition for maximum growth stability, and 
is advantageous for both the front and back menisci. It is clear from 
Fig. 94 that the trough must be tilted backwards, as shown, to allow coating 
to take place at the optimum substrate height. 
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In summary, the effects of meniscus shape and meniscus stability are very 
important in SCIM-coat ing. Because of their complexity, calculations are 
quite helpful Id specifying trough design and conditions for stable growth. 

Physical Proportion of SCIM Material - Structural characterization of SCIM- 

28 

coated layers was performed under SERI support. In summary , SCIM-coated 
layers appear to >e structurally identical to dip-coated layers, i.e , t.ie 
predominant structure is heavily twinned with {ill} twin planes perpendi- 
cular to tho surface and parallel to the <211> growth direction; the sur- 
faces are close to the jlio} planes. The grains, however, are generally 
larger in the SCIM layers but, as with dip-coated layers, they too are 
generally perpendicular to the layer's surface. 

Electrical Properties of SCIM Material - A complete characterization of the 
electrical properties of SCIM material was never performed. This work 
could no|t be satisfactorily performed until the SCIM coating machines were 
fitted with purified carbon parts. Purified carbon parts could no., on the 
other hand, be installed until all necessary coater modifications were 
nade to ensure satisfactory coatings. 

Preliminary diffusion length (L r ) measurements on SCIM samples did, however, 
reveal values of about 20um, which is comparable to the 15 to 40um values 
measured lr> dip-coated material by the same techniques (see Appendix A). 

Solar Cell Effort 


Fabrication Procedure - The SOC research and development program began as a 
silicon sheet materials investigation. One method for assessing the quality 
of the silicon sheet produced was to evaluate the performance of solar 
cells fabricated from this material. Because of the uniqueness of SOC, con- 
ventional fabrication procedures could not be used. That is, because the 
insulating substrate prevented direct contact to the back surface of the 
cell, special fabrication procedures had to be developed. These various 
procedures are described in SOC Annual Reports Nos. 1 through 5 and are 
summarized here. 

First Year - Photodiodes were fabricated on dip-coated SOC material. The 
best diode had an active-area conversion efficiency, n, of 2.5% (AMI) with 
no antireflection (AR) coating. Other characteristics of this diode are 
given in Table 19, along with the corresponding characteristics of the best 
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SOC cells fabricated during the second through fifth years. For the 2.5% 
efficient photodiode, the Junction was formed using a ;, 2 0 5 diffusion and a 
mesa etch process. Electroless nickel was used as the electrical contact 
to the p-type base region and silver paste was used to contact the 
diffused layer. 


Table 19. Best cell summary data. 
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Second Year - Interdigitated SOC cells were fabrication on solid (non- 

slotted) SOC substrates. As shown in Table 19, the best cell had a totnl- 

2 

area conversion efficiency of 3.86%, for a total of 1.73cm . The junction 
in this cell was formed by a P 2 °5 diffusion. The interdigitated cell pattern 
was defined by using photoresist and an Fe 2 0 3 mask. The cell area was de- 
fined by using plasma etching. The base contact was made by using a plati- 
num-silicon contact, followed by nickel and indium solder. The coll 
patterns are shown in Fig. 97. 

Third Year - Slotted substrate SOC cells were developed during this year. 

The best cell had a conversion efficiency of 7.5% (AMI, AH), for a cell 
2 

area of 4.08cm . Also developed were new diffusion and AH coating 
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Figure 97. SOC solar cells, (a) An array of diodes. 

The central square is a base contact, as well 
as the- surrounding outer region, (b) An 
i nterdi gi t a t ed contact solar cell with lcm~ 
active area, tc) A solar cell with 10cm 2 
active area. 
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technique. For the diffusion process, which used phosphine as the source 
gas, a more uniform temperature profile was achieved in the furnace and a 
baffle was added to increase phosphine gas uniformity. The diffusion furnace 
gas flow rates were also optimized. For the AR coating process a new 
vacuum system was developed. This (diffusion-pumped) system used a silicon- 
oxide evaporation source from Vacuum Atmosphere Company. 

During the third year, it was also recognized that a serious series re- 
sirtance could be Introduced into the cell using slot openings for back 
electrical contact purposes. An analysis of parasitic resistance was made 
and is discussed in Appendix B of this report. 

Fourth Year - The efficiencies of SOC cells were improved by optimizing 
the base layer impurity concentration. The results of the study are in- 
dicated by solid lines in Fig. 98. The dashed lines represent earlier data 

reported in SOC Quarterly Report No. 12. The more recent data indicate that 

16 3 

cell performance can be improved if the Ng, is greater than 10 /era . The 
recent results differ from the previous results primarily at the low and 
high base doping concentrations. That is, several new SOC cells were fabri- 
cated with concentrations of 2.3 x 10^ and 4.6 x lO^/cm"*. It is assigned 
that the recent results are more accurate than the previous results. In 
the remainder of this subsection we will refer only to the solid curves in 
Fig. 98. 

Figure 98(a) shows that the open-circuit voltage, V qc , increases linearly 

as a function of log N Q . For a given value of Ng, the center of the vertical 

line indicates the average value of V and the length of the line is 2o, 

oc 

where a is the standard deviation. Thus, for a base doping concentration 
16 3 

of 2.3 x 10 /cm , the average value was 0.536V, with a standard deviation 
of 0.005V. For the low and high values of Ng, two vertical lines are given. 
The lower line corresponds to the results of SOC Quarterly Report No. 12 
and the upper line to the recent results. The upper lines are sometimes 
slightly displaced to the right or left so that the various vertical lines 
do not overlap. The solid line for V qc indicates approximately how the 
average value of V qc should vary as a function of Ng, based on the recent 
measurement . 

As has been stated, Fig. 98(a) shows that V increases linearly as a 

Q C 

function of log Ng. This result is consistent with conventional solar cell 
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Figure 98. Performance of slotted SOC solar cells as a function 
of base doping concentration, for cells with AR coat- 
ing. Dashed lines correspond to results given in SOC 
Quarterly Report No. 12. Solid lines correspond to 
results given in SOC Annual Report No. 4. 
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theory. Figure 98(b) shows that the fill factor, FF, also increases linearly 

as a function of log Ng. The increase is due primarily to a corresponding 

decrease in serie? resistance. Figure 98(c) shows that, as N n increases, the 

16 

short-circuit current density, J gc « is relatively constant until N Q - 10 / 

ora** and then decreases for higher values of Ng. Finally, Fig. 98(d) shows 
that the conversion efficiency, n, increases monotonically as Ng increases 
for the range of values measured. 

During the fourth year, photodiodes were fabricated from layers of various 
thicknesses in an effort to determine how cell performance is affected 
by film thickness. The results are shown in Fig. 99. It should be noted 
that all data are without an AR coating. The results indicate that the 
efficiency increases by about 20% as the silicon thickness increases from 
lOOum to 200um. 

The best cell fabricated during the fourth year had a total area conversion 

2 

efficiency, n, of 10% (AMI, AR) , fo** a cell area of 4.1cra . 

Fifth Year - Cell performance was improved significantly by using a 
slow cooldown following the phosphorus diffusion. The diffusion temperatu’'** 
is 850°C for 80 minutes. The diffusion cycle is divided into a 10-minute 
preheat time, a 60-minute source time, and a 10-minute postheat. The slow 
cooldown procedure was adopted as a routine practice. 

To improve cell performance, a number of process variations were investi- 
gated and the principal ones are listed in Table 20. The major results 
of these experiments were discussed in SOC Quarterly Reports Nos. 14 and 15. 
These results are summarized below, along with some more recent results: 

1) Melt-Current Experiment - Eight SOC cells were fabricated to deter- 
mine the effect of an electric current through the melt during dip- 
coating. It was hoped that this current might reduce melt contamina- 
tion by influencing the transport of impurities from the substrate 
into the melt. Based on the results of the eight cells, an electric 
current does not seem to have an effect on melt contamination. 

2) No-Doping Experiment - Four SOC cells were fabricated to determine 
the effects of progressive melt contamination on the sheet resistance 
of the dip-coated SOC material. For the ninth and twelfth substrates 
dipped in Run No. 228, the sheet resistance was 770 and 550 ohms/ 
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Figure 99. Photodiode characteristics as a function 
of silicon thickness. The photodiodes 
were fabricated on unslotted SOC materials 
grown at various speeds. Squares are ex- 
perimental davt; points. 


Table 20. Experiments for increasing cell performance. 
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Experiment 
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1 

Melt current 
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No effect 

2 
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Graph of p vs N 

3 

Aluminum doping 
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sc 

4 

Heavy doping 

Determine effect of 
heavy doping 

Reduces J 

sc 

5 

Carbon substrate 

Measure melt quality 

Poor structure 

6 

Ion- imp lan tat ion 

Increase J 

sc 

Preliminary results 

7 

Hydrogenat ion 

Increase J 

sc 
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8 

Two-step 

Increase J 

sc 

Some improvements 

9 

Slow-cool 

Increase J 

sc 

Significant improvements 

10 

Etching 

Increase J 

sc 

No improvements 
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square, respectively. Cells fabricated using these SOC substrates had 
efficiencies in the range from 4.3 to 6.1% (AMI, no AR). This indi- 
cates that the mullite progressively contaminated the melt with a 
p-type dopant, shown to be aluminum in earlier experiments. Such 
contamination is a basic problem for the dip-coating technique. A 
similar problem should not exist for the SCIM-coating technique, 
since the mullite will not be in direct contact with the silicon melt. 

3) Aluminum-Doping Experiment - Eleven SOC cells were fabricated to 
determine if aluminum doping could produce higher values of J for 

o bC 

a given doping concentration. The highest J is 14.68 mA/cnr 1 (AMI, 

® C i g q 

no AR) for an aluminum doping concentration of 2.2 x 10 /cm . This 

2 

value is low compared with an average value of 17 raA/cra for boron- 
doped samples. Thus, aluminum doping does not seem to produce high 
values of J 

sc 

4) Heavy-Doping Experiment - Eight polar cells were fabricated using 

18 3 

SOC material with a boron-doping concentration of 10 /era , as 

16 3 

compared with the normal concentration of 5 x 10 /cm . Our objective 
was to extend the experimental results given in Fig. 4 of SOC 
Quarterly Report No. 14 in order to determine the upper limit for 
the base-doping concentration. The fabricated cells were of poor 
quality and only the J results were of value. The most interesting 

SC ry 

result was that the average J was 19 mA/cm , as compared with 
2 

21 mA/cm , the value associated with the linear regression line given 

o 

in Fig. 4 of SOC Quarterly Report No. 14. The 19 raA/cra value in- 
dicates that Auger recombination begins to limit the performance of 

18 3 

SOC cells for base-doping concentrations greater than 10 /cm , as 
expected. 

5) Carbon-Substrate Experiment - Several carbon substrates were dipped 
first in a run in order to learn more about melt contamination from 
the ceramic substrate. We had planned to make photodiodes on the 
resulting silicon layers and to compare the performance of these 
diodes with the performance of corresponding diodes on SOC material. 
Unfortunately, the grain structure of the silicon layers on the carbon 
substrates was rather poor and, as a result, we did not complete this 
experiment . 

6) Ion-Implant Experiment - Eight SOC cells were fabricated using ion 
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implantation and thermal anneal to form the n + emitter region. Arsenic 
was used as the dopant. For one set of SOC cells the annealing tem- 
perature was 850°C for 60 minutes. The emitter sheet resistance was 
approximately 32 ohms/square and the efficiency of one ion-implant/ 
anneal cell was 5.6%, as compared with the 6.7% efficiency of a 
corresponding cell used for a control. More work needs to be done to 
determine the proper annealing conditions. 

7) Hydrogenation Experiment - Results for this experiment are discussed 
in the Evaluation of SOC Material subsection of this report. 

8) Two-Step Experiment - Eighteen solar cells were fabricated to deter- 
mine the effects of a two-step diffusion process. The best two-step 

cell had a total-area conversion efficiency of 10.26% (AMI, AR) for 

2 

a cell area of 5cm . For one group of cells, the two-step process was 

29 

similar to the process described by DiStefano and Cuomo, as well 
as by Lindmayer. 30 This two-step process basically is as follows: 
First, a low-temperature phosphorus diffusion is used to diffuse 
phosphorus atoms down grain boundaries using a spin-on diffusion 
source and diffusion conditions of 600°C for 60 minutes. Second, the 
initial diffusion source is then removed and the normal high-tempera- 
ture phosphorus diffusion is used to form the p-n junction. For a 
second group of cells the initial spin-on source was used for both 
the low- and high-temperature diffusion. For a third group the order 
of the two diffusions was interchanged. The latter process seems 
to produce the highest efficiencies. 

9) Slow-Cooldown Experiment - Twenty-seven SOC cells were fabricated 
using a slow cooldown after the phosphorus diffusion. Typically, the 
cooling rate was approximately 3°C/min from the diffusion temperature 
of 850°C to 700°C. The slow cooldown seems to produce a significant 
increase in cell performance. 

10) Etching Experiment - Nineteen cells were fabricated from SOC material 
that was etched before the phosphorus diffusion. The main idea is to 
determine if cell performance can be improved by removing impurities 
that segregate to the silicon surface during growth. Various etches 
have been used, including KOH etch, Sirtl etch, CP4 etch, and plasma 


29 

DiStefano and Cuomo, Appl. Phys. Lett. 30, 351 (1977). 

30 

Lindmayer, 13th IEEE Photovoltaic Specialist’s Conf. , p. 1096. 
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etch. In all cases, we have removed approximately, 3ura of material 
from the surface of the SOC material. None of the etches produced a 
noticeable Increase in cell performance. 

During 1980 all SOC cells were fabricated using a well-established and re- 
producible process. For a group of 74 cells - called the 1980 baseline 
cells - the conversion efficiency was 9.6% with a standard deviation of 
0.5%. This process is described below: 

1) Initial Clean - The as received SOC material is first dipped into 
straight HF for 30 seconds and rinsoci in DI water. The silicon sur- 
face is scribbed using soap and water and rinsed in DI water. The 
SOC material is dipped in sulfuric peroxide ( , 3:1) for 5 
minutes and rinsed in DI water. It is then dipped in Aqua Regia 
(HC1/HN0 3 , 3:1) for 5 minutes and rinsed in DI water. Finally, it is 
dipped in straight HF, rinsed in DI water, and blown dry with dry Ng. 

2) Phosphorus Diffusion - The SOC material is baked at 120°C (tlO C) 
overnight and then placed in a phosphine diffusion furnace for 80 
minutes at 850°C. After the diffusion, the SOC material is cooled 
down at a rate of 3°C/minute from 850°C to 700°C and then withdrawn 
from the furnace. The sheet resistance is about 45 ohms/square. 

3) Back Contact - Photoresist is applied to the front surface and the 
silicon is etched on the back side in the slot regions, to remove 
completely the back surface n + region formed during the diffusion. 
About 70oX of platinum is sputter-deposited on the back surface and 
then sintered at 600°C for 30 minutes to form a PtSi contact. About 
2000A of nickel is sputter-deposited to provide a good surface for 
soldering . 

4) Front Contact - The front metal pattern is defined by using photo- 
lithography. About 80oX of titanium is sputter-deposited, followed by 
2G00X of nickel. Finally, the whole substrate is dipped in tin- 
indium solder at 170°C. 

5) Mesa Formation - The cell area is defined by using photolithography 

and by using plasma etching to form a mesa cell structure, with an 
2 

area of about 5cra . 
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6) AR Coating - About 600X of SiO is evaporated on the surface of the 
cell to complete the SOC cell structure. 

Solar Cell Performance - The year-by-year best performance level of cells 

produced since the inception of this program is given in Table 19. Near 

the end of this program (December, 1980), the best dip-coated cell had a 

total-area conversion efficiency, n, of 10.54% (AMI, AR) and the best SCIM- 

coated cell had an n of 7.6% (AMI, AR). Each cell has a total area of 
2 

5cra and uses a slotted substrate. The current-voltage characteristics of 
these two cells are given in Figs. 100 and 101. 

During 1980, we fabricated approximately 250 SOC cells from dip-coated 
material. Of these cells, 74 were fabricated under nearly the same condi- 
tion and the principle characteristics of these ce).ls - which are called 
the 1980 baseline cells - are .listed in Tables 21 and 22. The 10.54% cell 
noted above is not included in these tables, since this cell was fabri- 
cated using a nonstandard fabrication process which consisted of using an 
H 2 /Ar atmosphere during the growth of the SOC material. Distribution graphs 
for the 1980 baseline cells are given in Figs. 102 and 103, corresponding 
to conditions before and after AR coating. As can be seen, the average 
cell efficiency is 9.6% (AMI, AR), with a standard deviation of 0.5%. This 
deviation indicates that our standard cell fabrication process can produce 
very reproducible results. 


To increase our understanding of the SOC cell, empirical and theoretical 
modeling was performed. The results of this modeling are summarized here: 

In Quarterly Report No. 14 we presented a figure showing cell performance 
as a function of the base doping concentration for the 1979 baseline cells. 
This figure, which is based on 136 SOC cells, is reproduced here as Fig. 
104. The experimental results are indicated by the vertical liDes; the cal- 
culated results, by the solid curves. These impirical curves can be 
calculated as follows: First, linear regression is used to determine a re- 
lationship between J and N n . The result is 

sc B 


J sc - 1.48 log 


9.7 x 10 


31 


(61) 


where J is short-circuit current 
sc 

base layer doping concentration in 


2 

density in units of mA/cra and N 0 is 

-3 ** 

units of cm . Second, linear regression 
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Figure 100. Current-voltage characteristics 
of a slotted cell fabricated 
from dip-coated material. 


Figure 101. Current-voltage characteristics 
of a slotted SOC cell fabri- 
cated from SCIII-coated material. 






Table 21. Conversion efficiencies of the 1980 baseline 
SOC cells. 
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0.06 

0.7 

30.0 

6.81 

9.41 

18 

248 


4-211 

3.C 

O' 6 

0.06 

0.9 

30.0 

7.02 

9.94 

\9 

24? 


7- 1 1 1 

5.0 

016 

0.06 

0 • / 

20.0 

7.07 

10.07 

20 

24? 


7-211 

5.0 

016 

0.06 

0.8 

20.0 

5.24 

8.39 

21 

231 


3-111 

3.0 

016 

0.06 

0.9 

15.0 

7.06 

9.90 

22 

231 


3-211 

3.0 

016 

0.06 

0.9 

13.0 

6.28 

9.12 

23 

231 


8-111 

3.0 

016 

0.06 

0.8 

21 0 

6.48 

9.26 

24 

231 


8-211 

5.0 

016 

0.06 

0.7 

21 .0 

6.81 

9.78 

23 

232 


3-111 

3.0 

016 

0.06 

0.8 

14.0 

6.89 

9.56 

24 

252 


3-211 

3.0 

016 

0.06 

0.7 

14.0 

6.61 

9.78 

27 

232 


4-111 

3.0 

016 

0.06 

0.8 

19.0 

6.92 

9.87 

28 

257 


4-211 

3.0 

016 

0.06 

0.9 

19.0 

7.06 

9.96 

2? 

252 


6-111 

5.0 

016 

0.06 

0.7 

21.0 

7.11 

9.80 

30 

232 


6-211 

3.0 

016 

0.06 

0.7 

21.0 

7.26 

10.12 

31 

232 


7-111 

5.0 

016 

0.06 

0.6 

23.0 

7.00 

9.96 

32 

232 


7-211 

3.0 

016 

0.06 

0.7 

23.0 

6.59 

9.48 

33 

232 


?-11l 

5.0 

016 

0.06 

0.8 

19.0 

7,05 

9.90 

34 

232 


9-21 1 

3.0 

016 

0.06 

0.9 

19.0 

7.16 

10.22 

33 

233 

- 

3-ill 

3.0 

016 

0.06 

0.7 

18.0 

7.22 

10.42 

36 

233 


3-211 

5.0 

016 

0.06 

0.7 

18.0 

7.17 

10.38 

37 

233 


4-111 

3.0 

016 

0.06 

0.8 

22.0 

7.19 

9.97 

38 

253 


4-211 

5.0 

016 

0.06 

1.0 

22.0 

7 • 22 

9.82 

3? 

234 


1-111 

5.0 

016 

0.06 

1.4 

21.0 

6.62 

9.15 

40 

254 


1-211 

3.0 

016 

0.06 

0.6 

21.0 

6.92 

9.97 

41 

234 


6-111 

3.0 

016 

0.06 

0.7 

21.0 

6.78 

10.19 

42 

254 


6-211 

3.0 

016 

0.06 

0.7 

21.0 

6.98 

10.11 

43 

234 


9-111 

3.0 

016 

0.06 

0.6 

26.0 

7.08 

10.03 

44 

234 


9-211 

3.0 

016 

0.06 

0.6 

26.0 

6.85 

9.96 

45 

255 


2-111 

3.0 

016 

0.06 

0.6 

21 .0 

6.90 

9.58 

46 

233 


2-211 

3.0 

016 

0.06 

0.6 

21.0 

7.09 

9.82 

47 

233 


3-111 

3.0 

016 

0.06 

0.6 

22.0 

?.oe 

9.87 

48 

253 

- 

5-211 

5.0 

016 

0.06 

0.7 

22.0 

7.58 

10.21 

4? 

233 

- 

7-111 

3.0 

016 

0.06 

0.6 

27.0 

6.97 

9.91 
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Table 21. Conversion efficiencies of the 1980 baseline 
SOC cells (concluded). 


Total Area 
Convaraion 
la»« Effic Uncut 

C* 11 »••• Dtp DA Ihiat lafora Aftar 

Nuabar Doping Spaad Product RatUtanca AR Alt 


(atoM/cc) (ca/tac) tohat-tqca) (ohat/M) (X) (X) 


itaa 

SO 

233 - 7-211 

3.0a«0l4 

0.04 

0.7 

27.0 

4.92 

10.03 

31 

234 - 4-111 

S.OatOlO 

0.04 

1.3 

32.0 

4.92 

9.77 

32 

254 - 6-211 

3.064016 

0.06 

1.4 

32.0 

7.0C 

9.90 

S3 

237 -13-111 

3.064014 

0.06 

0.9 

29.0 

6.70 

9.44 

34 

237 -13-211 

3.064014 

0.04 

1.4 

29.0 

6.80 

9.33 

33 

259 - 7-111 

3.064016 

0.06 

0.9 

22.0 

6.94 

9.52 

34 

258 - 7-211 

3.064016 

0.04 

0.7 

22.0 

7.14 

10.13 

57 

239 - 8-111 

3.064016 

0.04 

1.1 

35.0 

6.78 

9.76 

39 

258 -12-111 

3.064016 

0.06 

0.9 

34.0 

6.67 

9.38 

59 

238 -12-211 

5.064016 

0.06 

1.0 

34.0 

4.37 

9.13 

60 

473 - 7-211 

3.064014 

0.06 

0.9 

30.0 

7.14 

9.82 

61 

4?4 - 8-1 1 1 

3.064016 

0.06 

0.8 

23.0 

7.10 

9.88 

62 

474 - 8-211 

5.06401 6 

0.06 

0.9 

23.0 

6.87 

9.38 

63 

474 -IC-11I 

3.064016 

0.04 

1.0 

27.0 

6.92 

9.60 

64 

474 -10-211 

5.064016 

0.06 

1.2 

27.0 

6.73 

9.32 

65 

474 -11-111 

5.0t4016 

0.06 

0.7 

21.0 

7.02 

10.20 

66 

474 -11-211 

3.064016 

0.06 

0.8 

31.0 

6.79 

9.64 

67 

473 - 5-111 

3.064014 

0.06 

1.5 

31.0 

6.97 

9.73 

68 

473 - 3-211 

5.0a >014 

O'. 06 

1.2 

31.0 

7.36 

to. 10 

69 

477 - 7-111 

3.064016 

0.06 

0.9 

46.0 

6.89 

9.47 

70 

477 - 7-211 

3.064016 

0.04 

0.9 

44.0 

6.88 

9.33 

71 

477 -11-111 

3.064016 

0.04 

0.9 

37.0 

6.85 

9.73 

72 

477 -11-211 

5.064016 

0.06 

0.9 

37.0 

7.21 

9.86 

73 

478 -11-111 

3.064016 

0.06 

0.9 

29.0 

6.39 

9.72 

74 

479 -11-211 

3.064016 

0.06 

0.8 

29.0 

4.98 

9.86 
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Table 22. Characteristics of recent AR-coated 
SOC cells. 


fotal-Araa 



Cell 

lie 

Ooc 

Fill 

Jtc 

Elf. 


t 

UA) 

10) 

Factor 

(«A/*<jca) 11) 

lt«n 

1 

22* - 3*lll-0 

•0.40 

0.551 

0.728 

16.08 

6.43 

1 

22* - 3-111-1 

114.80 

0.570 

0.709 

22.9* 

9.28 048) 

2 

22* - 3-21 ' -0 

83.70 

0.554 

0.74 ; 

14.74 

6.93 

2 

22* - 3-21 1 - 1 

117.30 

0.570 

0.719 

23.4* 

9.62 (AR) 

3 

238 - 8-111-0 

82. *0 

0.53* 

0.737 

14.52 

6.95 

3 

238 - 8-1 11-1 

114.10 

0.S7C 

0.733 

22.82 

9.3* (Aft > 

4 

238 - 8-2M-0 

82.10 

0.54* 

0.6*4 

16.42 

3.93 

4 

238 - 8-2. 1-1 

1 .4.20 

0.5*0 

0.6*5 

22.84 

8.31 (Aft) 

V 

J 

241 - 4-1 1 1-0 

76.40 

0.554 

0.737 

13.68 

6.40 

3 

241 - 4-111-1 

107.40 

0.56* 

0.721 

21.48 

8.77 (Aft) 

* 

241 - 4-211-0 

77.70 

0.53* 

0.649 

13.54 

3.37 

6 

241 - 4-211-1 

110.00 

0.551 

0.*74 

22.00 

8 . , 7 ( Aft ) 

7 

244 - 3-111-0 

81.50 

0.552 

0.673 

16.30 

3.73 

1 

244 - 3-111-1 

114.20 

0.559 

0.656 

22.84 

6.40 (Aft) 

9 

244 - 3-211-0 

82.30 

0.5a3 

0,772 

16.46 

7.13 

P 

244 - J-2I1-I 

114.90 

0.5*5 

0.731 

22.98 

9.75 (Aft) 

9 

244 -11-111-0 

81.60 

0.357 

0.6*3 

16.32 

6.03 

9 

244 -11-111-1 

110.70 

0.568 

0.670 

22.14 

8.43 (Aft) 

10 

244 -1 1- 21 1-0 

83.00 

0. 5o2 

0.739 

16.60 

7.08 

10 

244 -11-211-1 

113.60 

0.572 

0.733 

22.72 

9.79 (Aft) 

1 t 

244 -12-111-0 

79.60 

0.560 

0.752 

13.92 

4.70 

1 1 

244 -12-111-1 

117.20 

0.565 

0.71 • 

23.44 

9.42 (Aft) 

12 

244 -12-211-0 

81.10 

0.332 

0.638 

16.22 

3.71 

12 

244 -12-211-1 

117.90 

0.360 

0.644 

23.56 

8.50 (Afti 

13 

248 - l-l 1 1-0 

81.00 

0.546 

0.449 

14.20 

5.92 

13 

248 - 1 -I 1 1-1 

111.20 

0.561 

0.690 

22.24 

8.61 (Aft) 

14 

248 - 1-211-0 

82.40 

0.550 

0.753 

16.48 

6.82 

14 

248 - 1-211-1 

113.90 

0.5*7 

0.743 

22.78 

9.62 i Aft ) 

13 

2 48 - 2-M1-0 

79.40 

0.556 

0.760 

15.88 

6.71 

13 

248 - 2-111-1 

113.50 

0.5*8 

0.718 

<> 3 • 1 0 

9.42 (Aft) 

16 

248 - 2-211-0 

79.60 

0.526 

0.469 

13.92 

3.93 

16 

248 - 2-211-1 

116.30 

0.548 

0.31 : 

23.30 

4 . aO i Aft ) 

1? 

248 - 4-1 1 1-0 

81.30 

0.330 

0.761 

16.2* 

6.81 

17 

248 - 4-111-1 

111.50 

0.3*3 

0.750 

22.30 

9.41 <AR> 

18 

248 - 4-211-0 

83.53 

0.532 

0.741 

16.70 

7.02 

18 

248 - 4-211-1 

115.80 

0.56* 

0.738 

23.16 

9.94 (Aft) 

1? 

249 - 7-111-0 

84.10 

0.360 

0.733 

14.82 

7.07 

19 

249 - 7-111-1 

121.10 

0.373 

0.723 

24.22 

10.07 (Aft) 
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Table 22. Characteristics of receri AR-coated 
SOC cells (continued). 


i t #11 

Cell 

Nun6*r 

lit 

(M«> 

Voc 

(V) 

fill 

factor 

Fotal-Art* 

Jit iff. 

<ft4/S0Cft> (!) 

20 

20 

244 

244 

• 7-211-0 

* 7-211-1 

78.30 

120.00 

0.340 

0.340 

0.418 

0.424 

13.70 

24.00 

3.24 

8.39 (68) 

21 

21 

23* 

23. 

- 3-111-0 

- J-lll-l 

84.40 

120.20 

0.344 

0.373 

0.739 

0.714 

14.88 

24.04 

7.04 

9.90 (48) 

22 

22 

251 

231 

- 3-211-0 

- 3-211-1 

84.80 

121.10 

0.340 

0.370 

0.441 

0.441 

14.94 

24.22 

4.28 

9.J2 (48) 

23 

23 

231 

231 

- 8-111-0 
- 8- 111-1 

80.00 

120.40 

0.364 

0.348 

0.718 

0.477 

14.00 

24.08 

6.48 

9.26 (41" 

24 

24 

231 

231 

- 8-211-0 
- 8-211-1 

81.60 

122.70 

0.343 

0.347 

0.741 

0.703 

14.32 

24.34 

4.81 

9.78 (48) 

23 

23 

232 

232 

- 3-111-0 

- 3-1 1 1 -1 

80.20 

117.30 

0.364 

0.372 

0.739 

0.713 

16.04 

23.44 

4.89 

9.34 (48) 

20 

24 

232 

232 

- 3-211-0 

- 3-211-1 

80.10 

120.20 

0.362 

0.370 

0.734 

0.714 

14.02 

24.04 

4.61 

9.78 (48) 

2? 

27 

232 

232 

- 4-111-0 

- 4-111-1 

82.30 

118.40 

0.349 . 
0.373 

0.739 

0.725 

16.46 

23.68 

4.92 

9.87 (48) 

28 

28 

232 

232 

- 4-211-0 

- 4-211-1 

82.90 

120.10 

0.363 

0.371 

0.734 

0.726 

14.48 

24.02 

7.06 

9 96 (48) 

2? 

2? 

232 

252 

- 6-111-0 
- 6-IM-1 

83.30 

119.30 

0.369 

C.374 

0.746 

0.708 

14.70 

23.90 

7.11 

9.80 (4R) 

30 

30 

232 

232 

- 6-211-0 
- 4-211-1 

83.80 

119.40 

0.367 

0.578 

0.764 

0.733 

14.76 

23.86 

7.26 

10.12 (48) 

31 

31 

232 

232 

- 7-111-0 

- 7-111-1 

80.70 

120.50 

0.366 

0.377 

0.747 

0.716 

16.14 

24.10 

7.00 

9.96 <4R> 

32 

32 

252 • 
232 • 

• 7-211-0 
■ 7-211-1 

80.20 

120.00 

0.336 

0.370 

0.736 

0.693 

16.04 

24.00 

6.39 

9.48 (48) 

33 

33 

232 ■ 
232 - 

4-111-0 

4-111-1 

82.60 

116.20 

V.370 

0.580 

0.749 

0.727 

16.32 

23.24 

7.05 

9.80 (48) 

34 

34 

- € 

232 - 

4-211-0 

4-211-1 

83.60 

122.30 

0.570 

0.578 

0.751 

0.723 

lo.72 

24.46 

7.lo 

10.22 ( AR ) 

35 

33 

253 - 
253 - 

3-111-0 

3-111-1 

84,80 

123.60 

0.365 

0.373 

0.734 

0.736 

16.94 

24.72 

7.22 

10.42 (48) 

36 

36 

253 - 
253 - 

3-21 i -0 
3-21 i-l 

84.80 

123.10 

0.564 

0.372 

0.750 

0.723 

16.96 

23.02 

7.17 

10. J6 (4R1 

32 

37 

253 - 
253 - 

4-1 ! 1-0 
4-111-1 

84.30 

119.10 

0.563 

0.578 

0.757 

0.724 

16.86 

23.82 

7.19 

9.97 (48) 
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Table 22. Characteristics of recent AR-coated 
SOC cells (continued). 


1 t»M 

Cell 

Nu«t*r 

lie 
( mA ) 

Voc 

<V> 

Pill 

Factor 

lotal-Arca 
J«c £ff. 

(AA/sqcnl (X) 

38 

38 

233 - 4-211-0 
253 - 4-211-1 

84.80 

121.10 

0.560 

0.574 

0.760 

0.706 

16.96 

24.22 

7.22 

9.82 (AR) 

39 

39 

234 - 1-111-0 
254 - 1-HI-l 

80.30 

118.20 

0.560 

0.569 

0.736 

0.680 

16.06 

23.64 

6.62 

9.13 (AR) 

40 

40 

254 - 1-211-0 
254 - 1-211-1 

79.00 

119.10 

0.363 

0.570 

0.778 

0.734 

15.80 

23.82 

6.92 

9.97 (Aft ) 

41 

4! 

254 - 6-111-0 
254 - 6-111-1 

83.40 

119.40 

0.564 

0.372 

0.721 

0.746 

16.68 

23.88 

6.78 

10.19 (AR) 

42 

42 

254 - 6-211-0 
254 - 6-21 1 -1 

81.40 

117.40 

0,362 

0.567 

0.763 

0.760 

16.28 

23.48 

o.9B 

10.11 (AR) 

43 

43 

254 - 9-111-0 
254 - 9-111-1 

82.70 

117.90 

0.565 

0.580 

0.738 

0.733 

16.54 

23.58 

7.08 

10.03 ( AR) 

44 

44 

254 - 9-2H-0 
234 - 9-211-1 

80.90 

119.30 

0.561 

0.373 

0.753 

0.726 

16.18 

23.86 

6.85 

9.96 <AR> 

45 

45 

255 - 2-111-0 
255 - 2-111-1 

80.20 

112.80 

0.564 

0.573 

0.763 

0.741 

16.04 

22.56 

>.90 

9.58 (AR) 

46 

46 

235 - 2-211-0 
255 - 2-211-1 

82.80 

114.70 

0.565 

0.373 

0.738 

0.747 

16.56 

22.94 

7.09 

9.82 (AR) 

47 

47 

235 - 5-111-0 
255 - 3-111-1 

62.30 

119.50 

0.569 

0.573 

0.756 

0.718 

16.46 

23.90 

7.08 

9.87 (AR) 

48 

48 

253 - 3-211-0 
255 - 5-211-1 

86.00 

120.60 

0.571 

0.57? 

0.772 

0.731 

17.2C 

24.12 

7.58 

10.21 (AR) 

49 

49 

235 - 7-111-0 
253 - 7-111-1 

81.70 

116.50 

0.565 

0.371 

0.753 

0.745 

16.34 

23.30 

6.97 

9.91 (AR) 

50 

30 

235 - 7-211-0 
235 - 7-211-1 

82.10 

116.70 

0.363 

0.573 

0.748 

0.750 

16.42 

23.34 

6.92 

10.03 (AR) 

5 1 
51 

256 - 6-111-0 
256 - 6-111-1 

86.00 

126.10 

0.561 

0.574 

0.717 

0.675 

17.20 

25.22 

6.92 

9.77 (AR) 

52 

52 

256 - 6-2H -A 
256 - 6-211-1 

87.80 

125.70 

0.557 

0.573 

0.716 

0.687 

17 56 
25.14 

7.00 

9.90 (AR) 

S3 

53 

257 -13-111-0 
257 -13-111-1 

83.50 

116.10 

0.551 

0.563 

0,732 

0.737 

16.70 

23.72 

6.70 

9.64 (AR) 

54 

54 

237 -13-211-0 
237 -13-21 1-1 

85.40 

117.30 

0.551 

0.363 

0.723 

0.708 

17.08 

23.46 

6.80 

9.35 (AR) 

35 

35 

253 - 7-111-0 
238 - 7-111-1 

83.00 

117.00 

0.560 

0.569 

0.746 

0.715 

16.60 

23.40 

6.94 

9.32 (AR) 


70 



Table 22. Characteristics of recent AR-coated 
SOC cells (concluded). 


Tot»l-Ar*» 




Ctll 

lie 

Voc 

Fill 

JfC 

Ef f 



Nunbtr 

\ «A ) 

(V) 

F»ctor 

UA/tqCA) 

(X) 


lt»* 

• as 








56 

258 

- 7-211 0 

84.00 

0.561 

0.760 

16.80 

7.16 


56 

258 

- 7-211 -1 

117.10 

0.57C 

0.759 

23.42 

10.13 

UR) 

57 

258 

- 8-111-0 

83.50 

0.560 

0.725 

16.70 

6.78 


57 

258 

- 8-111-1 

119.40 

0.571 

0.716 


9.76 

UR) 

58 

258 

-12-111-0 

85.00 

0.555 

0.707 

17.00 

6.67 

• 

58 

258 

-12-111-1 

119.30 

0.567 

0.708 

23.86 

9.58 

UR) 

59 

258 

-12-211-0 

84.00 

0.549 

0.691 

16.80 

6.3? 


59 

258 

-12-211-1 

118.00 

0.559 

0.692 

23.60 

9.13 

UR) 

60 

47J 

- 7-211-0 

85.10 

0.563 

0.745 

17.02 

7.14 


60 

473 

- 7-211-1 

121.20 

0.574 

0.706 

24.24 

9.82 

UR) 

61 

474 

- 8-111-0 

83. aO 

0.565 

0.752 

16 

7.10 


61 

474 

- 8-111-1 

120.80 

0.569 

0.719 

24. 1.. 

y .88 

UR) 

62 

474 

- 8-211-0 

81.30 

0.561 

0.753 

16.26 

6.87 


62 

474 

- 8-211-1 

119.70 

0.565 

0.708 

23.94 

9.58 

UR) 

63 

474 

-10-111-0 

82.10 

0.560 

0.752 

16.42 

6.92 


63 

474 

- 1 0- 1 1 1 - 1 

119.70 

0.568 

0.706 

23.94 

9.60 

(AR) 

64 

474 

-10-211-0 

82.50 

0.559 

0.730 

16.50 

6.73 


64 

474 

-10-211-1 

120.30 

0.569 

0.695 

24.06 

9.52 

UR) 

65 

474 

-1 1-1 1 1-0 

82.10 

0.562 

0.761 

16.42 

7.02 


65 

474 

-11-111-1 

116.50 

0.573 

0.764 

23.30 

10.20 

(AR) 

66 

474 

-11-211-0 

80.10 

0.559 

0.758 

16.02 

6.79 


66 

474 

-11-211-1 

114.30 

0.570 

0.749 

22.86 

9.64 

(AR) 

67 

475 

- 5-111-0 

87.10 

0.564 

0.709 

17.42 

6.97 


67 

475 

- 5-111-1 

119.50 

0.575 

0.708 

23.90 

9.73 

(AR) 

68 

475 

- 5-211-0 

86.50 

0.568 

0.749 

17.30 

7.36 


68 

475 

- 5-211-1 

121.10 

0.574 

0.726 

24.22 

10.10 

(AR) 

69 

477 

- 7-111-0 

82.30 

0.558 

0.750 

16.46 

6.89 


69 

477 

- 7-111-1 

1)3.60 

0.564 

0.739 

22.72 

9.47 

(AR) 

?0 

477 

- 7-211-0 

84.50 

0.558 

0.729 

16.90 

6.88 


70 

477 

- 7-211-1 

114.30 

0.563 

0.725 

22.86 

9.33 

ur ; 

71 

477 

-11-111-0 

82.80 

0.557 

0.743 

16.56 

6.85 


71 

477 

-tl-lll-l 

1 '8.70 

0.571 

0.719 

73.74 

9.75 

(AR) 

7 2 

477 

11-211-0 

85.90 

0.557 

0.753 

17.18 

7.21 


7? 

477 

-11-711-1 

121.00 

0.570 

0.715 

24.20 

9.36 

(AR) 

73 

478 

-11-111-0 

81.90 

0.555 

0.725 

16.38 

6.59 


73 

478 

-1 1-1 1 1-1 

116.20 

0.570 

0.734 

23.24 

9,72 

(AR) 

74 

476 

1 

1 

?\J 

1 

O 

82.90 

0.559 

0.753 

16.58 

6.98 


74 

470 

-11-211-1 

117.00 

0.568 

0.742 

23.-0 

9.86 

(AR) 
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SHORT-CIRCUIT CURRENT DENSITY lmA/cm J ) CONVERSION EFFICIENCY <%l 

Figure 103. Distribution graphs for the 1980 baseline 
SOC cells for conditions after AR coating. 
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is used to determine an empirical relationship between Jq^ and Ng. For the 
data given in Table 4 of SOC Monthly Report No. 34, the result is 


'01 


9.6 

7^.54 

1 B 


(62) 






Figure 104. Cell performance as a function of base doping con- 
centration, for the 1979 baseline cells. The 
characteristics of these cells were reported in 
SOC Annual Report No. 4 and in SOC Quarterly 
Report No. 14. Vertical lines are experiment re- 
sults; solid curves are calculated results. 

2 

where Jq^ is reverse-saturation current density in units of mA/cm . Third, 
the theoretical expression for the specific series resistance, R s A, as given 
in SOC Quarterly Report No. 12, is used, with an empirical constant of 1.75 
added to approximate the measured results. The resulting expressiou is 
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(63) 


R A - 2.25 + 
s 


3.35 x 10 


15 


TT 


B 


where R s A is in units of ohms-cm . For Eq. (63 > a growth velocity of 3.6 
cm/min is assumed. Finally, the current density is assumed to be given by 

[q(V 

rr l 


J * J 


01 


ex r 


"FT 


JR s A) 




' / 

I- 1 

+ J 02 

exp 1 


q( v - JR a A) 

— jjpp 


! v - JR S A) 

— rr~ 

i sh 


- J 


sc 


(64) 


-4 2 

where is 10 mA/cm and n is 2. Equations (61) through (64) give the 
empirical curves of Fig. 104. Clearly, the empirical curves approximate the 
experimental results. 


Now we turn from empirical modeling to theoretical modeling. Our objective 
is to understand more completely the performance limitations of SOC cells. 

In particular, we want to determine a relationship between the short-circuit 
current, J gc> and the minority-carrier diffusion length, L nB . assumiug 
AMI conditions. 


To begin, we consider only generation in the base region and monochromatic 
light. These restrictions will be removed as we progress. Moreover, we will 
assume that the geometry is one-dimensional and that the effects of the 
grain boundaries can be neglected. In other words, we initially consider an 
SOC cell structure that approximates a one-dimensional, single -crystal cell 
structure having a base diffusion length corresponding to that of SOC mate- 
rial. The basic cell configuration is shown in Fig. 105, along with a list 
of the various parameter values. A list of notation is given in Table 23. 
With the noted assumptions, conventional semiconductor theory yields the 
following five equations for low-level conditions: 

n 2 q(- v 1E ) 

p(0) = NT- exp ’ E T- " 


n i 2 C 1(" V 1E ) 
n(0) = tj — exp — r- T ‘ i — 

lN AB 


(65) 


( 66 ) 
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Electron diffusion length In base region, cm. 
Donor concentration in the emitter region, cm 

Acceptor concentration 1. the base region, cm 

- 3 

Intrinsic concentration, cm 

_ 3 

Electron concentration, cm 


Open-circuit voltage. 

Electrostatic potential drop across the emi t ter-base 
space-charge region minus the corresponding drop at 
equilibrium, V. 

Distance from the boundary between the space charge 
region and the base region, cm. 

Thickness the base region, cm. 

Absorption :oef f ic ient . cm 






(67) 


V JE + K S X D 


In • qA p(0) + I (0) 
D L pE n 


( 68 ) 


oiL 


I n (0) =* qAF( 1 - R) — 


nB 


< a L nB “ X > 


aL 


Gosh(X Q /L nB ) - exp( -aXg ) 


nB 


sinh(X Q /L nB ) 


qA,-SS n(0) 
L nB 


(69) 


For monochromatic light having a wavelength of 0.882um, these equations 
given the current-voltage characteristics shown in Fig. 106. (Also shown are 
the corresponding current-voltage characteristics for AMI light conditions, 
which will be discussed below.) For monochromatic light, the minority 
carrier profile in the base region is given by 



Figure 106. Calculated current-voltage char- 
acteristics for monochromatic light 
of wavelength of Q.882um and for 
AMI light . 
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aiil 

n ab 


! n<°) 

^bAb^AII 

- Y sinh (-ax) 


aL p YJ sir.h 


"nB 


♦ sS^ + Y 

W AB 


cosh 


aB 


(70) 


where 


Y 


° Fl \ - B > 

< A nB - 


(71) 


Equation (70) is plotted in Fig. 107 for three values of wavelength. The 
figures demonstrate the following: As the wavelength increases, the slope 
at x - 0 decreases, indicating a decrease in collection efficiency. The 
latter decrease is approximately given by 



Figure 107. Minority-carrier profiles in the 
base region for three values of 
wavelength: 0.823, 0.882, and 0.982um. 
The profiles were calculated at the 
maximum power point. 
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(72) 


^col 


L nB a 


1 ♦ L 


nB 


a 


which is a well-known approximation. In particular, Eq . (72) is used to 
evaluate the diffusion lengtn from LBIC measurements of n cQl . 

Now we turn to AMI conditions. Our primary objective is to calculate a 

theoretical upper limit for J gc and to compare this value with the measured 

values of J for SOC cells. For AMI conditions, ve must solve the following 
sc 

standard differential equation by numerical methods: 


d n n(x ) 
' ,2 
dx L nB 



0 


(73) 


where 


G(x) - jXaP(X) exp [-ax] dX 

J o 


(74) 


Here P(X) is the AMI solar power spectral density given as a function of 
the photon wavelength, . For Eq. (74) we have assumed that the AR coating 
is perfect (i.e., R = 0) and that the thickness of the emitter region is 
negligible. The boundary conditions for Eq. (74) are 

n(x) =» n(0) at x * 0 (75) 

and 

n(x)-0atx a X B (76) 

Using numerical methods and Eqs. (69) through (72), we can determine a 
relationship between 1^(0) and n(0), which is the generalization of the 
monochromatic result expressed by Eq. (69). The generalized relationship 
can then be used with Eqs. (65) through (69) to give the AMI current-voltage 
curve shown in Fig. 1C 6. Figure 108 shows the corresponding minority-carrier 
profiles for J gc> P max , and V qc 


178 



Figure .108. Minority-carrier profiles for J , 

P , and V . The profiles correspond 
max oc * 

to the AMI curve given in Fig. 106. 

2 

For all calculations, the incident power density, P in> is 88.92 mW/cm , 

which is the accepted value for AMI sunlight. This value is about 10% lower 

than the 100-mW/cra 2 value generally used for cell measurements with an 

ELH lamp. For the figure, the metallization coverage was neglected, which 

means the J is about 10% higher than possible with a typical metallization 
sc , 

pattern. The low value of P in approximately compensates for the neglect o. 

the metallization coverage. Thus, AMI results shown in Fig. 106 approxi- 
mately represent measurement conditions using an ELH lamp with an incident 

2 

power density of 100 mW/cm . 

So for all calculations, we have assumed a diffusion length of 50pm and 
negligible series resistance. Figure 109 shows the calculated cell charac- 
teristics for values of L nB in the range from 10um to lOOuro, where -he 
solid lines correspond to a specific series resistance of 1 ohm-cr, ' .nd 
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whore the dotted lines correspond to negligible series resistance. For the 
1 ohra-cm case, the calculated maximum efficiency is 11.2% for L 0 - 25uro 
13.3% for L R g ■ 50um, and 14.5% for L Q g ■ lOOum. The measured characteristics 
of the SOC cell shown in Fig. 110 are indicated by squares in Fig. 109. 

For this cell there is fair agreement between theory and experiment. 


io A 



§«. 

< 
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: 0 . 4 . 

0.2 • 

0 . l 
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60 100 
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Figure 109. Solar cell characteristics as a function 
of base diffusion length for a bare 
thickness of 200um. Dotted lines assume 
no series resistance, solid lines include 
l-ohm/cm^ specific resistance. The squares 
in the latter correspond to cell No. 
255-5-211-1. 
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COST ANALYSIS 


Int rodue'.lon 

The cost of silicon-coated ceramic has been determined at least annually by 
a sequence of different methods. Prior to 1978, our cost estimates were 

31 

derived according to procedures practiced in one of Honeywell's factories. 

32 

When the SAMIS methodology and the IPEC formula became available, our 

own cost procedures were no longer used. Early in 1978, a direct comparison 

was made using the Honeywell procedure versus the IPEG formula, with the 

same direct-cost input. The two prices differed by less than 5%, which is 

within the accuracy of the data. Our most recent estimates, including 

those contained in this report, are according to the revised IPEG-2 
•»3 34 

formula. A later revision, only recently "vailable, has not been used 
in our calculations. 

In over 4 years and several iterations of price estimations, the projected 

production costs for the SOC technology have not changed significantly in 

equivalent dollars. For example, an added-vaiue requirement of $13.01/(11^ 

(1980 dollars) was calculated in 1977, 31 compared with the IPEG-2 estimate 
2 

of $12. 99/m reported in this document . A calculation by SAMIS-III based 

2 

on 1979 data yielded an added value of $13.70/ra . Some major changes have 
occurred in the projected technology, particularly in the transition from 
dip-coating to SCIM-coat ing . It was reasonable to expect a lower labor 
intensity for SCIM, which is a continuous process, than for dip-coating, 
which is a complicated hatch operation. However, the labor saving was 
offset when the orginally projected pulling speed of 1.0 cm/sec war. re- 
vised to 0.25 cm/sec as a more realistic goal. 


31 ERDA/JPL 954356-77/3, Annual Report No. 2, published September 30, 1977. 

TO 

R. W. Aster and P.. G. Chamberlain, "IPEG: A precursor and an Adjunct to 
SAMIS-III, Version One," JPL Do ment 5101-33, September 10, 1977. 

R. W. Aster, "IPEG-2: Improved Price Estimation Guidelines," in JPL Docu- 
ment 5101-142 (December 1979), pp. 355-357. 

34 

P. J. Firnett, "Improved IPEG Computer Program User's Guide," JPL Document 
5101-156, Rev. A., November 17, 1980. 


This subsection presents the cost analysis for SOC, reflecting (a) 1986 tar- 
get technology, and (b) the technology as it exists today. The analysis of 
even today's technology is hypothetical, since the following assumed condi- 
tions have not yet been proven 

e Continuous automated operation at stated duty cycle 
# High-volume substrate availability at quoted price 
e 92% yield 

e Simultaneous coating of two substrates 

The required prices are expressed <.n both dollars per square meter and cents 

per peak watt. For the latter, we assume a net subsequent processing yield 

(cell and module product icn ) 3S of (0.95) (0.995) • 0.945. A final graph 

Indicates sheet cost versus efficiency that must be satisfied by our process 

in order to achieve a 70-cents/W module price. 

P 

The major econo: .ic assumptions will be more clearly understood in terms of 
the summary process description, which follows next. 

Process Summary 

Slotted mullite substrates, measuring 12.5cm by 50cm and 1mm thick, are 

automatically conveyed from the substrate unpacking area to the CARBCOAT 

equipment. A suspension of pulverized carbon in isopropyl alcohol is 

applied by roller to one surface of the substrates at a linear speed of 

3 cm/sec. Eight substrates ( 100cm total width) are processed simultaneously 

o 

for a per-unit throughput of 1.8 m /min. Compressed air is used to expedite 
drying. The process is highly automated, requiring one production operator 
per eight units. All material handling is automatic. 


35 

R. W. Aster, "Price Allocation Guidelines," JPL Document 5101-68, Rev. A, 
January 15, 1980. 
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The 12.5cm x 50cm carbonized substrates are automatically coveyed from the 
C 4RBC0AT line to the SILCOAT furnaces. The substrates approach the furnaces 
side-by-side In pairs. A pair of substrates enters a furnace by way of the 
heating zone, which is roughly 40cm long, so that the carbonized surfaces 
are facing downwards. The substrates are hea'.ted from room temperature to 
roughly 1420°C in 2.7 minutes, the time it would take for a point on the 
substrate to transverse the heating zone at 15 cm/rnin. The coating chamber 
immediately following is maintained near 1420°C, and molten silicon is 
applied from a transversely oriented trough to the underside of the sub- 
strates. Next is the cooling zone, which is also about 40cm long. Both 
heating and cooling zones are carefully designed to minimize residual 
stresses in the coated substrate. For this purpose, it is desirable to 
maintain constant temperature gradients in the moving substrate in the 
elastic-plastic transition between 500°C and 1000°C. A separate thermal re- 
quirement is imposed immediately following the silicon solidification 
position,' at the beginning of the cooling zone. Normally, active cooling 
(forced convection plus radiation) would be required here to accommodate the 
latent heat of fusion at the 15-cm/min linear growth speed; however, 
passive (radiative) cooling would suffice with an elongated solidification 
front as in cold-substrate or hot-substrate growth techniques. The pre- 
sent cost analysis assumes passive cooling. A slight positive pressure of 
inert atmosphere (argon) is maintained in all three zones. 

A nominal 7-year life is assumed on the cabinet, substrate drive and support 
structures, power supplies, and electrical and cooling water connections. 

Electrical furnace elements, insulation packages, crucibles, and troughs 
are replaced on a regular weekly shutdown schedule. Preventive maintenance 
on substrate drive is assumed to take place during these scheduled shut- 
downs. Scheduled shutdowns are 12 hours per week, allowing for cooling and 
startup. Only the crucibles are replaced weekly. 

The process will be automated to an extent requiring instrument monitoring 
but no material handling. One operator will monitor 12 furnaces, or 24 sub- 
strates in process. The minimal vertical height requirements will permit 
~ating furnaces to be stacked vertically in groups of six. 

From SILCOAT, the cooled, coated substrates are conveyed to INSPECT stations, 

where they are subjected in sequence to two tests called SCAN and RESIST. 

2 

The equipment for each test has an output of 12 pieces (0.75m ) per minute. 

I 
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The SCAN equipment monitors the coated surfaces for nonspecular reflection 
into a cylindrical collector, which would be an indication of a coating gap 
or some other gross fault. The raster scan used a 2.5mm spot, with trans- 
verse deflection provided by rotating mirrors. The sweep speed is 625 cm/sec, 
or 50 sec/sv/eep across the panel. Longitudinal displacement of successive 
sweeps is provided by the 12.5-cm/sec advance. A total of 5 seconds is 
required: 4 seconds scan plus 1 second between panels. The RESIST equipment 
makes resistance measurements by means of a four-point probe. The test is 
performed at four locations on each panel, 1 second each, plus 1 second 
between panels, totaling 5 seconds per panel. Direct labor is calculated on 
a piece rate basis and averages out to 0.648 production machine operator 
and 0.936 inspector for each pair of test machines. 

P rincipal Assumptions for Target Technology * 

2 " 2 

• Substrates are available at $5. 78/m , slotted, in 5M m /yr quantities. 
This is according to an IPEG-2 estimate made by a potential vendor 

in December, 1980. 

2 

• The graphite-in-propanol suspension is $6. 96/gal, or $0.054/ra as 
applied, in accordance with the actual experimental data. 

• Carbon coating apparatus is $10,200 per unit net, adjusted for sal- 
vage value and removal/installation cost. The useful life is 7 years. 

2 

• Carbon coating output is 1.8m /minute. 

• Silicon-coating units are $50,800 net, adjusted for salvage value and 
removal/ installation costs. Useful life is 7 years (See Process 
Summary) . 

• Crucible replacement cost is $80 per week per unit. 

• Furnace heating elements are replaced at 6-month intervals at a cost 
of $1200 per unit. 

• Insulation packages are replaced annually at a cost of $240 per unit. 

3 

• Argon gas consumption is 0.117 ft /rain., which is about one-sixth of 
the experimental rate. The reduced figure is the estimated seepage 
from redesigned entrance and exit ports. 

• • Electrical consumption in SILCOAT is calculated to be 3.67kW: 3.39kW 
required to heat the substrates and silicon and melt the silicon at 
the stated throughput; plus 0.28kW radiative losses. 


♦Items marked with two bullets differ from Progress-to-Date technology. 
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• Polysilicon cost is assumed to be $14/kg. 

• SILCOAT duty cycle is 85%. Of the downtime, 7.1% is for scheduled 
maintenance (Process Summary) and 7.9% is unscheduled maintenance. 

2 

• • Silicon-coating rate is 0.0368 m /min based on two panels coated 

simultaneously (25cm total width) at an advance speed of 15 cm/min. 
There is a 1cm gap between panels in the advance direction. 

• • Production area is reduced by stacking SILCOAT units in vertical 

groups of six. 

• • One production machine operator monitors 12 machines (two stacks). 

• Inspection requires SCAN and RESIST at $50,800 and $6,100 each, 
respectively. A useful life of 10 years is assumed for each machine. 

2 

• Inspect output is 0.75 m /min for each test. Production operator load 
is based on 0.9 standard hour per 1000 pieces; inspector load is 
based on 1.3 standard hours per 1000 pieces. 

• Inspect yield is 92%. This is the net residue from CARBC0AT , 

SILCOAT, and INSPECT processes; i.e.: 

(tested and approved S0C panels) r (uncrated substrates) » 0.92 
IP^G-2 Estimates fo” Target Technology 

The entries in Table 24 are direct costs for an annual production of 2.5 
million square meters of sheet SOC, not adjusted for yield in cell fabrica- 
tion and subsequent module manufacturing steps. Projected low-cost technology 
is postulated, as summarized by the listed assumptions. 


Table 24. Cost data for target technology. 


Task 

EQPT 

SQFT 

DLAB 

MATS 

UTIL 

Carbon coating 

$ 34,358 

649 

$ 61,805 

$ 687 , 800( a ) 

15 , 706 , 500(b) 

$ 14,131 

Silicon coating 

8,405,947 

3969 

851,827 

1,378, 350(c) 
8.881.464(d) 

252,987 

Inspect 

501 , 500 

2396 

862,491 

— 

7,345 

Total 

$8,941, 805 

7014 

$1,776,123 

$16. 531, 650(e) 
$25,413, 114( f ) 

$274 , 463 


(a) Carbon. ( b ) substrates, (c) argon, crucibles, furnaces. Insulation, 

(b) polysilicon, (e) excluding silicon, (f) including silicon. 
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Manufacturing prices estimates are based on the IPEG-2 formula: 


Price 


0.57 

or 

0.52 


| |eqp*i| + iob| 


109 |SQFT 


4-2.8 DLAB 4-1.2 


MAT 4- UTIL) 


/QUAN (77) 


In our case, vha £QPT coefficient, 0.52, applies to INSPECT, where the anti- 
cipated equipment life is 10 years; the coefficient 0.57 applies to the 
other two tasks, for which equipment life is 7 years. We have also applied 
the DLAB coefficient of 2.8 because our direct labor costs do not include 
fringe benefits. Reduced to price per square meter, the project technology 
estimates are: 

$12. 99/m 2 added value 
$17.25/m2, including silicon 


For an assumed 11% encapsulated cell efficiency and subsequent processing 

35 

yields (cell and module fabrication) totaling 94.5% , :he respective costs 

on a peak watt basis are: 


12.5 cents/W added value 

P 

16.6 cents/W including silicon 

P 

The price estimation formula can also be used to calculate price sensitivity 
data, or price changes resulting from a change in any cost input parameter. 
For illustration purposes, it is convenient to denote the cost input changes 
in terms of percent of nominal values. For example, a poly silicon price of 
$16/kg would be expressed as: 

1-g i £ - 100 x - 14.3% (78) 

since the nominal price is $14/kg. The price sensitivity is given relative 
to ceramic (CER) and silicon (SIL) costs, silicon coating time (SCT), coating 
machine cost (CMC), coating operator load (COL), and argon cost (ARG): 


$ - 17.25 4- 

0.0754 

% CER 

4- 

0.0426 

% SIL 


4- 

0.0353 

% SCT 

4- 

0.0192 

“’CMC 


4- 

0.00954 

% C0L 

4- 

0.00083 

T 

l AR G 

(79) 
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The sensitivit ies are the coefficients in Eq. 79, or the slopes of the 

corresponding gi phs in Fig. 111. These graphs, of course, intersect at 0%, 

o 

giving the nominal sheet value, $17. 25/m . The largest theoretically 
possible cost saving is represented by negative 100%, meaning that a parti- 
cular parameter has zero cost value. Hence, the intercepts on the left-hand 
vertical axis give the hypothetical sheet price in the zero-coat limit 
(e.g., the "added-value" is seen as the silicon graph intercept) of $12.99. 



CHANGE IN COMPONENT COST 
(% OP PROJECTED AMOUNT) 


Figure 111. Sensitivity profiles for total sheet 
value according to projected tech- 
nology. Abscissa represents relative 
(%) change in component cost (e.g., 
+100 and -50 denote doubling halving, 
respectively ) . 


Prcgress-To-Date Technology 

Another major category of economic analysis deals with progress-to-date 
technology instead of projected low-cost technology. A number of parameters 
have been changed to reflect our perception of progress-to-date, parti- 
cularly with respect to silicon coating throughput and operator load. 
Throughput is based cn a pulling speed of 0.0625 cra/sec, operator load is 
0.5, and machines are placed two to a stack instead of six. Compared with 
projected technology, direct silicon coatings costs are thereby magnified 
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by the following factors: 


• Coating machines x4 

• Materials (except Si) x4 

• Operators x24 

• Floor space xl2 


Table 25 giver airect costs for progress-to-date technology, with the as- 
sumptions according to Table 24 except as modified above. 


Table 25. Cost data for progress-to-date technology. 


Task 

EQP". 

SQFT 

DLAB 

MATS 

UTIL 

Carbon coating 

$ 34,358 

649 

$ 61,804 

$ 687.800(a) 

15, 706, 500(b) 

$ 14,131 

Silicon coating 

33,623,788 

47,628 

20,443,848 

5,513,400(0 
8, 881, 464(d) 

232,987 

Inspect 

501,500 

2,396 

862,491 

— 

7,345 

Total 

$34,150,646 

50,673 

$21,368,144 

$21, 907, 700(e) 
$30,789, 164( f) 

$274,463 


Ca) through (f): See Taule 24. 


Based on the values of Table 25, the IPEG2 formula yields: 

2 

$44.57/ra added value 
2 

$48. 83/m , including silicon 

For price per peak watt, we assume an encapsulated cell efficiency of 8.4% 
and a subsequent net processing yield of 94.5%: 

56.1 cents/W„ added value 
P 

61.5 cents/W^, including silicon 
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Price Allocat ion 


A commercial readiness goal of 70 cents per peak watt has been set for the 
year 1986, to be met by one or more low-cost, flat-plate, solar voltaic 
concepts that have demonstrated technological readiness by 1982. The JPL 
Price Allocation Guidelines 35 were developed as a tool for estimating 
price goals for low-cost sheet silicon technologies, such as Cz ingot and 
slicing, HEM ingot and slicing, EFG, dendritic web, and SOC. The basic 
assemption is that a iven set of pric'} allocations will have been met by 
the other ma1>r (nonsheet) task categories. These allocations are stated 
in 1980 dollars per square meter as follows: 


• 

Cell fabrication 

$21/ra 2 

of 

cells 

• 

Encapsulation materials 

$14/m 2 

of 

cells 

• 

Module assembly 

$14/m 2 

of 

module 


Assuming the remaining prices are for sheet and silicon, an allocation for 
SOC can be calculated from: 


$/W (SOC) * $/W total - $/W cell - $/W p encaps materials - $/W mod (80) 
assembly 

where 

$/W (SOC) « $/W p sheet + $/W p silicon (81) 

and $/ n (TOTAL) is the total module price, in dollars per peak watt. The 
price per peak watt can be expressed in terms of the price per m 2 , the 
silicon thickness, t (mils), silicon yield, Y gi , cell yield, Y n , module 
yieid, Y mfg , packing efficiency, n p , and encapsulated cell efficiency, n . 


• Silicon: 


$/» - $/kg |t • u • D ) / , I • r. e * Y • T cell • Y mfg ) 


Sheet : 


$/W ■ $/ru 2 sheet/ I • n • Y 

P 0 


e coll raf g i 


• Cell : 

$/» p - $/« 2 cell/ | I . n„ . Y„ fg | 


(84) 


t 

* 

* 

♦ 

l 

/ 

T 

4 


* 

\ 
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(85) 


• Encapsulation Materials : 

$/W p - $/rn module/ (i . . n p J 

• Module Assembly : 

$/Wp « $/m^ module/ | I • n 0 • n ) (86) 

where 

Insolation (I) - 1000 W/m 2 

Silicon density (D) - 2330 kg/m 2 

Unit conversion (u) - 2.f4 < 10 -5 m/mil. 

For SOC, { . -0.80 and t ■ 4 mils. With the suggested values, 5 Y 

si **“ ’cell 

” °' 95 ' Y m fg “ 0-995, r.nd silicon cost of $14/kg, Eq. 80 becomes 

(0,00106 $/m 2 sheet + 0.00438)/n e - $/W total - (0.0211 0.0147 ♦ (87) 

0.0147)/n P 

The left hand side of Eq. 87 is price per peak watt of sheet plus silicon. 
This equation is solved for the allocated sheet price, giving: 

$/m 2 sheet - 943(n e )($/W p total) - 51.8 (88) 

Figure 112 is a graph of $/m 2 sheet as a function of n e for cost goals of 
70, 60, and 50 cents per peak watt. For n - 0.11, cost allocations are 
$20. 82/m and $10. 44/m for 70 cents and 60 cents per peak watt, respec- 
tively. Equation 88 can also be used to calculate the required encapsulated 
cell efficiency to meet a given $/W p totai, assuming $12.99/ra 2 as found 
from our IPEG-2 analysis of the projected technology: 


70 

cents/W 

P 

requires 

n e " 

9.8% 

30 

cents/W p 

requires 

n e ' 

11.5% 

50 

cents/W p 

requires 

n e ■ 

13.7% 


Summary 

2 

The projected added-value cost of SOC is $12. 99/m , with the substrate 
material accounting for 58% of this amount and the actual silicon coating 
process contributing 29%. The remaining 13% is attributable to carbon 
coating and inspection. Of the silicon coating costs (29%), the relative 
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amount by the SAMIS categories 

are : 

• 

Capital Equipnont 

50% 

• 

Labor 

25% 

• 

Materials (non-silicon) 

17% 

• 

Manufacturing Space 

5% 

• 

Electricity 

3% 


All of these except electricity (i.e., 28% of the added value) scale in- 
versely as the silicon-coating throughput, which is thus a major cost 
driver after substrate cost. (See also Fig. 111.) 

The $1 2 . 99/ra 2 figure is consistent with the 70 cents/W target, assuming 
an encapsulated cell efficiency of at least 9.8%. For 60 cents or 50 

certs/Wp , the efficiency would be required to be at least 11.5% or 13.7% 
respectively . 



Figure 112. Allocated sheet price vs. encapsulated cell ef- 
ficiency. Numerical values are shown for n "0.11. 
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CONCLUSIONS 


From the work performed during this RM) program, we conclude that; 

• Silicon sheets of solar-cell quality can be produced by the SOC 
process, whereby thin layers of silicon are coated onto one surface 
of large-area ceramic (mulllte) substrates. 

2 

• Large-area (1000cm ) mulllte substrates, which contain slot openings 
for electrical contact purposes, can be inexpensively produced in 
large quantities. 

• Such substrates can be compositionAlly modified to obtain a thermal 
expansion coefficient which is similar to that of silicon. The slight 
dissimilarity leaves the silicon coating In a state of compression. 

• The SOC material can be satisfactorily produced by either a dip- 
coating or a continuous SCIM-coating process. 

• If the coating throughput rate of the SOC silicon sheet process is 

2 

increased to 350 cm /min, the process falls within DOE's 1986 
economic goals. 

• Carbon can be inexpensively applied to one surface of the substrate 
using a colloidal carbon suspension. 

• Thin SOC layers, by virtue of unidirectional solidification, contain 
columnar silicon grain as large as 5rara in width »r.d several cm long. 

• The single grains are heavily twinned, but the twin boundaries as 
well as the normal boundaries are, in general, perpendicular to the 
surface of the layer. 

• The SOC coatings are predominately {33l} surfaces with a growth direc- 
tion of <211>. This texture permits {ill} twin planes to occur per- 
pendicular to the substrate and propagate as the grain grows. This 

26 

texture is similar to that observed in EFG ribbon silicon. 

• The silicon-ceramic bond is mechanical in nature and has proven to be 
highly dependable in practice. 
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• Slot openings in the substrate which are perpendicular to the growth 
direction can be satisfactorily covered with silicon with very little 
perturbat ■" of the grain structure. 

• The thickn^s* of tile layer is inversely related to the melt tempera- 
ture and is also proportional to the Inverse of the coating speed 

squared wnen the heat of fusion is removed from the growing layer by 

radiation alone. 

• Experimentally, it was demonstrated that with additional -onvective 
cooling, coatxng rates commensurate with the project's economic goals 
are achievable. 

e The mullite substrates are slightly soluble in molten silicon, causing 
the melt to become progressively contaminated with the dtp-coating 
process. This is minimized with the SCIM technique by reducing the 
contact time and area the substate has with th«. melt. 

• The surfaces of SOC coatings are smooth and su 1 ny and need no further 
preparation other than a standard cleanup procedure prior to forming 
the p-n Junction. 

• SOC solar cells can be reproducibly produced with an average conver- 
sion efficiency of 9.6% (AMI, AR). 

• The s lot openings in the cerauic substrates must cover approximately 
50% of the total substrate area to avoid the introduction of serious 
series resistance in the finished solar cell. 

e The highest conversion efficiency of a dip-coated SOC ceil was 10.54% 
(AMI, AR). 

• The highest conversion efficiency of a SCIM-coated SOC cell was 7.6% 
(AMI, AR). 

• The highest efficiencies were produced with a PH 3 diffusion at 850°C, 
following by a slow cooldown at about 3°C/min. 

• Theoretical modeling of SOC solar cells shows that present cell per- 
formance is limited by diffusion length, and that 13% SOC cells are 
possible with a diffusion length of 50um. 


k 
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• The average diffusion length (L) in SOC material is 15 to 25uro. The 
diffusion length within good grains is three tiroes the average dif- 
fusion length, L. Closely spaced grain boundaries, high dislocation 
densities, and subsurface grain boundaries have been identified as 
structural causes of the reduction in L. 

e Cost estimates of the SOC procoss are $13/rr. 2 added value, or $17. 25/m 2 
including silicon at $14/kg. Based on 11% module efficiency, this 
corresponds to a sheet cost of 16.6 cents/W p . The analysis shows that 
the cost of ceramic is the largest cost driver in the SOC process. 
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APPENDIX A 

EVALUATION OF SILICON-ON-CERAMIC MATERIAL FOR LOW-COST SOLAR CELLS* 


SUMMARY 

In silicon-on-ceramic (SOC) solar cells, the short-circuit current is less 
than that of single-crystal solar cells because of (1) recombination at the 
grain boundaries, and (2) recombination within grains. These two effects 
can be separated by light-beam- induced-current (LBIC) measurements of the 
short-circuit current using a finely focused scanned beam from a mono- 
chromator. At any point on the surface of a cell, an effective minority- 
carrier diffusion length can be determined from the spectral response at 
long wavelengths (0.8 to l.Oum). At each grain boundary, an effective 
surface recombination velocity can be determined. The LBIC technique was 
used to evaluate the effectiveness of hydrogen plasma treatment for reducing 
grain boundary recombination in SOC material. 

INTRODUCTION 

The silicon-on-ceramic (SOC) process produces thin layers of silicon on 

A- 1 

inexpensive ceramic substrates. Layers r.re produced by unidirectional 
solidification of molten silicon, and are polycrystalline with grain 
widths on the order of a millimeter and lengths of several centimeters. 

SOC layers have been produced by a dip-coating technique as well as by a 

A_ 2 

"Silicon-Coating with Inverted-Meniscus" (SCIM) technique. Solar cells 
made from SOC material have demonstrated 10.5% conversion efficiencies 
(AMI), using ELH lamp simulator. 

To improve the efficiency of SOC solar cells, it is important to assess the 
relative importance of grain boundaries and impurities. To separate these 
two effects, light-beam- induced current (LBIC) measurements of the short- 
circuit current have been made using a micron-sized light beam from a mono- 
chromator. The LBIC measurements are closely related to solar cell per- 
formance, since the technique identifies the various spatial and spectral 
contributions to solar cell current. 


•Paper by J. D. Zook, Honeywell Corporate Technology Center , Bloomington , MN. 
Published in Proc. of 1980 Photovoltaic Solar Energy iont ., Cannes, France 
(1980), p. 565T 
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PROCEDURE 


The experimental arrangement is shown in Fig. A-l. The exit slit of a mono- 
chromator is focused on a pinhole aperture which in turn is focused on the 
sample by a 15X, N.A./0.28 reflecting objective. For maximum resolution, 
the spot width on the sample can be reduced if the image of the exit slit 
does not fill the pinhole. Viewing of the sample is done with a pellicle 
beamsplitter and viewing eyepiece. Scanning is done with a motor-driven 
stage on which the sample is mounted. The absolute spectral quantum effi- 
ciency, Q x , is determined from the ratio of the sample current to that 
of a calibrated, highly uniform silicon photodiode inserted in front of 
the sample. 



Figure A-l. Schematic design of the LB1C system. The chopper is 
optional and was not used for the results reported 
in this paper. 


From a knowledge of Q x at any given wavelength in the near-IR, the value of 
minority-carrier diffusion length, L, can be determined. It is convenient 
to use the peaks in the spectrum of the xenon lamp (see Fig. A-2). For cells 


202 





that are thick compared with the penetration of light, A-3 the quantum 
efficiency is given by:'* - '* 

(1 - R) u x L 

Q i * 

where R is reflectivity and is the absorption coefficient. Unfortunately, 
R varies from sample to sample because of variations in the oxide thickness. 
In order to deduce R from the data, the values of L determined at each 
wavelength are averaged, and the standard deviation of 6L is determined. 

The value ol 1 - R is then chosen to give the lowest relative error, 6L/L. 
This procedure gives reasonable values of R for antireflection-coated (AR) 
cells as well as for uncoated cells. 



WAVELENGTH (pm) 

1 . 0.992320 

2 . 0.979970 

3 . 0.941201 1 
0.961338 J °' 946 

4 . 0.916286 

5 . 0.904646 

6 . 0.896226 

7 . 0.881941 

8 . 0.823163 

9 . 0.784202 


TiRUM READING (ARBITRARY UNITS) 


Figure A-2. Spectrum of the 150W xenon lamp used 

in the experiment. The wavelengths are 
handbook values for low-pressure 
xenon discharges. 

The value of L depends on the values assumed for a.. The most reliable pub- 

A 

lished data in the region of interest (0.76 to l.Oym) appear to be those 
of Runyan A-4 and of Weakliera and Redfield A- ^, whose results do not. differ 
significantly. Average values of at the strongest peaks of the xenon 
lamp are given in Table A-l. 
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Table A-l. Inverse absorption lengths. 


1( um) 

um) 

0.992 

125 

0.980 

100 

0.946 

59 

0.916 

38 

0.882 

26 

0.823 

14 


At the long wavelengths, a^L<<l in Eq.(A-l)and the LBIC response provide 
a direct mapping of the effective diffusion length. By integrating the 
LBIC response for a cross-grain scan (across the narrow dimensions of the 
grains), an average value, L, can be determined for the whole cell. 


RESULTS 


LBIC scans of an SOC cell at five wavelengths are shown in Fig. A-3. The 
scans include four well-isolated grain boundaries. The light spot size was 
about 5um x 12pm, with the narrow dimension perpendicular to the grain 
boundaries to give maximum resolution. Interestingly, the depths of the 



i • 098»n< BASEL INC 


SCAN DISTANCE, a 

Figure A-3. Typical LBIC scans on SOC material. The 

dashed lines at grain boundary 2 indicate 
the limits used in deriving the inte- 
grated current loss at the grain boundary. 
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minima at the grain boundaries are not sansitive to the spot size or shape, 
as evidenced by the fact that the scan is virtually the same even when the 
spot is significantly out of focus. The slgnal-to-noise ratio in the d-c 
system is greater than 100:1 inside the grains. The intensity of the light 
spot is between one to two suns for the wavelengths of Table A-l. The value 
of L « 39um within a grain was deduced by the procedure given above. 

An important issue in modeling grain boundary recombination is the question 
of linearity of the boundary conditions at the grain boundaries. Because 
the potential barrier heights at grain boundaries depend on the density of 
filled traps, it is reasonable to expect that grain boundary recombination 
velocity, s, might depend on the light intensity . An experiment was 
performed to evaluate the effect of light intensity on grain boundary re- 
combination. The beam was attenuated with neutral-density filters, and the 
quantum efficiency was measured within grains and at grain boundaries. The 
result showed no change in Q within grains over the range 0.01 to 10 suns, 
and about a <.0% increase in Q at grain boundaries over the same range. 

These results imply that s is essentially constant, so that the principle 
of superposition should be valid. Thus, the short-circuit current of the 
cell is given by the focused spot (LBIC) response, integrated over the cell 
p.rea, and over the solar spectrum. 

Since s is constant, the analytical theory of LBIC reponse ~ is applicable 
to SOC material. The normalized value of s is given by S ■ sL/D, where D 
is the minority-carrier diffusivity. The parameter S can be evaluated by 

A _ 7 

integrating the total loss in current due to a grain boundary. ‘ For 
grain boundary 2 in Fig. A-3, this is the area within the dashed lines. The 
results of the procedure are shown in Fig. A-4. Each data point in Fig. A-4 
represents a different wavelength, and each data point is consistent with 

5 

a value of S 3 25, corresponding to s 3 1 x 10 cm/sec, which is typical 
of SOC material. 

LBIC scans of a number of SOC cells show that there is a characteristic 

"good-grain" diffusion length, L max , which is dependent on doping and cell 

processing. We assume that L provides a measure of the chemical defects 

max 

density (density of intragrain impurities with deep- level electron states). 

There is also a characteristic "worst case" grain boundary diffusion 

length. L. min> which is typical of many of the grain boundaries. The value 

of L . is lower than the value of L at isolated grain boundaries (as in 
min 
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FI*. A- 3 ) , and is typically a factor of 10 or 20 less than the value of 
L 

max 


1/d ♦ A) 

1.0 0.0 0.8 0.7 0.0 0.5 0.4 0.3 0 2 0.1 0 



Figure A-4. Effective width, w, of a grain boundary as a func- 
tion of absorption coefficient, following Ref. A-7. 

The effective width times the length of the boundary 
gives the equivalent loss in cell area for the short- 
circuit current. 

Effective diffusion lengths are strongly influenced by structural defects. 

Low values of L are usually observed in regions having closely spaced grain 
boundaries that are readily visible on th9 top surface. In most cells, 
there are also regions which have an LBIC response lower than that which 
would be expected by examination of the top surface. In order to understand 
the reasons for such poor response, samples have been cut, polished, and 
etched to reveal the subsurface structure. In some cases, this structure 

reduces the value of L. More frequently, however, the low values of L are 

5 2 

associated with high dislocation densities (>10 /cm ) and dislocation 
clusters at stacking faults. 

The average diffusion length, L, evaluated from the integrated LBIC response, 

is typically one-third of L_ . From the theoretical relationship between 

max r 

L and J , it follows that J is reduced by 11 to 16% due to structural 
sc sc 

defects. Typically, this reduction is greater than the reduction due to the 

difference between L and the values of L in semiconductor-grade Czocb- 

max 
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ralaki material. Thus, structural defects appear to be more important than 
chemical defects in limiting the performance of present SOC solar cells. 


Th* LBIC technique has also been used to evaluate grain boundary passivation 

in SOC material. The hydrogen passivation treatment was performed by C. 

A-8 

Seagar and D. Ginley at Sandia Laboratories- The LBIC scans before (0) 
and After (A) hydrogen treatment are showr in Fig. A-5. The overall 
increase in Q is due to an lnadvertant oxide that was present after hydro- 
genation. It is clear that different degrees of passivation occur at dif- 
ferent grain boundaries but that most grain boundaries are significantly 



A 



Figure A-5. Effect of hydrogen plasma threatrrent ir. LBIC re- 
sponse at three wavelengths. In the photo at 
the top, the horizontal lines are growth stria- 
tions; vertical lines are grain boundaries. Scans 
before AR coating are denoted by B; after, by A. 


207 


improved. It is interesting to note that the nonuniformity within a grain is 
also reduced by hydrogen treatment. Correlations of LBIC response within 
grains are due to variations in dislocation densities. Thus, hydrogenation 
apparently also reduces recombination at dislocations. 

In summary, the LBIC technique is a useful tool for evaluating diffusion 
lengths within grains, grain boundary recombination velocities, the effects 
of processing on cell performance, and the relative importance of structural 
defects vs. chemical defects in limiting performance of polycrystalline 
silicon solar cells. 
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APPENDIX B 

PARASITIC RESISTANCE ANALYSIS 


INTRODUCTION 

The integration of our sheet silicon with a nonconducting substrate presents 
a special problem in making the electrical contacts to tne base layer. 

These connections could be made at the front surface, in which case the 
active photovoltaic regions would have to be confined to small mesa struc- 
tures. The relatively high r.on-active surface fraction would seriously 
reduce the module efficiency. Back surface contacts are possible only 
through slots or holes in the substrate. Since these contacts are necessarily 
stripes or dots rather than continuous area, the photocurront must travers*- 
a longer average path from the Junction to the back contacts. The series 
resistance of the base layer therefore can not be neglected. The purpose 
of this appendix is to (a) predict the consequences of parasitic series re- 
sistance with respect to the cell efficiency, and (b) estimate a tolerable 
limit on electrode spacing based on calculated series resistances in the 
solar cell. 

The basic electrical structure of the ai licon-on-ceramic (SOC) concept is 
illustrated in cross section in Fig. B-l. For definiteness, an n f on p 
structure is assumed, giving the cell polarity shown in the figure. The 
illuminated front surface is provided with regularly spaced striped elec- 
trodes which are modeled as a top metallic layer contacting the n* diffused 

layer through an intermediate metal-semiconductor interface. The interface 

2 

is characterized by a specific contact resistivity, p ' , in ohms-cm . Photo- 
generated minority carriers which diffuse to, and transverse the junction, 
are responsible for the photocurrent.. This photocurrent enters the structure 
through the interconnected (negative) front contacts, sprer-io through the 
diffused layer, crosses the Junction, and converges through the base layer 
to the (positive) back contacts. A carbonized layer on the ceramic sub- 
strate is approximately in the plane of the front surface of the substrate, 
and connections are brought out through slots in the ceramic. This is at 
variance with an earlier representation (SOC Annual Report No. 2), where it 
was assumed that the silicon would fill the slots up to the back surface 
of the substrate and the electrical contacts would be in the plane of that 
surface. However, our experience with both dip-coated and SCIM-coated sub- 
strates has shown the present version to be correct. The spacings of the 
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two se: i of contacts are not necessarily the same. It is presumed through- »> 

out the report that the back contact spacing, b, may exceed the trout con- 
tact spacing, b', and the ratio b/b‘ can be any positive number. A unit 
of the periodic structure will be taken as the b x w rectangular area 
bounded by the lines centered on adjacent back electrodes, and having arbi- 
trary width, w, perpendiculai to the plane of Fig. B-l. Current flow is 

» * 

assumed everywhere parallel to the plane of Fig. B-l. This amounts to 
Ignoring the collecting effect of the contact bars which connect the 
stripes.®”* 



I I 


Figure B-l. One unit of the periodic structure of an 
SOC cell. Features (not drawn to scale). 

(a) contact metal; (b) metal-semiconductor 
interface; (c) diffused layer; (d) junc- 
tion; (e) base layer; (f) carbon layer; 

(g) ceramic substrate; (h) external connection. 

Parasitic resistances occur at contact interfaces and in the front and back 
(diffused and base) layers. These not only reduce the maximum power avail- 
able to the load, but shift the maximum power point toward both lower 
voltage and lower current density. This may be illustrated qualitatively by 
means of Fig. B-2. The upper curves gives the hypothetical ideal illuminated 
cell characteristic (no parasitic losses). In this case, the load voltage 
is the same as the voltage, \' D> across the junction, here assumed to be 
uniform. Maximum power density (product VpJ) is labelled W Q *. The lower 
curve represents a particular non-ideal characteristic due to a specific 

B_ i 

R. J. Handy, "Theoretical Analysis of the Series Resistance of a Solar 
Cell," Solid State Electronics 10, 765 (1967). 
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parasitic resistance, Rp. The available load voltage differs from V D by the 
series IR drop; i,e., V A - V D - J bw R p . Maximum power density on this curve 
is labelled W*. 



V 

Figure B-2. Determination of optimum solar cell operat- 
ing point according to lumped model for 
parasitic resistances. 


The foregoing illustration makes use of two assumptions: that the internal 
resistances can be lumped and that V Q and J are constant along the 
junction. These assumptions cannot be strictly true due to the obviously 
distributed nature of the device. The parasitic IR effect varies from a 
maximum at the middle of a unit to near zero at the ends so that V p is not 
exactly constant. Moreover the bulk and contact components of the resistance 
are not strictly separate because of so-called current crowding in tbe 
semiconductor layers near the contacts. Fortunately, if certain reasonable 
auxiliary conditions are satisfied, the assumptions are justified to a 
good approximation. Specifically, the dissipated power should be relatively 
small, and the contact spacings b and b' should be large compared with the 
dif f used-layer thickness and base-layer thickness, respectively, and the 
metallized areas should be small fractions of the total area. 
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In the subsections that follow, ar, approximate lumped-resistance model is 
developed and some of its implications with respect to cell performance 
are discussed. 

THE SOLAR CELL JUNCTION PARAMETERS 

The illuminated solar coll characteristic (Fig. B-2) is assumed to 
following the single-exponential model 

.1 - J Q - J L (e AV D - 1) (B-l ) 

locally at any given point in the junction plane*. J Q is the illuminated 

short-circuit current density and J. the reverse-saturation current density. 

1 1 

The parameter A has an ideal theoretical maximum of q/kT ■ 38.7 volt but 
in silicon solar cells is always considerable smaller. J^, J^, and A are 
strongly dependent on the material and junction fabrication technology. The 
choice of numerical values is therefore somewhat arbitrary and subject to 
speculation. This should be borne in mind wheu interpreting the calculated 
results included in this report. The output may depend critically on the 
parameters. The present choice is based on short-circuit density - 0.03 
a/cm“, open-circuit voltage - 0.55V, and A ■ 30. This value of J Q corres- 
ponds to AMI sunlight and is not corrected for reflection losses. is 
then found to be 

AV a o 

.»! - JJ (e oc - l) - 2.05 x 10' a a/cm (B-2) 

These parameters yield a power density maximum - 12.9 mW/cm^ and fill 
factor - 78%. 


B-2 

H. J. Hovel, S olar Cell Semiconductors and Semlmetals . Vol. II, K. K. 
Willardson ancTX! C7 Beer , Eds . , Academic iPress, M.Y7 , Chap. 3. 

♦The deplet ion- layer width may not be negligible. However, the word "plane" 
here implies that current passes through the depletion layer perpendicular 
to its boundaries and with negligible ohmic effects. 
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POWER DISSIPATION 


The ideal performance figures are degraded in practice by series and shunt 
parasitic resistances and by front electrode masking. The shunt component 
is not included here, since its causes B-3 are related more to diode fabri- 
cation than to device geometry considerations and should therefore respond 
to improved processing techniques. Electrode masking is omitted also; if 
the masked area is not too great, a multiplicative masking factor on J Q 
should serve this purpose within the accuracy of the model. 

Dissipation occurs serially through front contact interfaces (R c '), front 

silicon layer (R^). back silicon and carbon layers (R^), anc * back contact 

interfaces (R c ). The lumped equivalent circuit for a single-cell unit is 

shown in Fig. B-3, and total parasitic resistance is IL ■ R ' + R, ' t R, 

f c L L 



V 


D 


Figure B-3. Lumped components of parasitic resistance 
of a cell unit. 


B-3„ . 

Hovel, op. 


cit . , 


p. 50. 
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The operating characteristic J vs . V A for a given Rp 
B-2 ) may be expressed parametrically in terms of the 


D ' 


as 


AV r 


J o - 


J 1 (e 


- 1 ) 


(lower curve of Fig. 
Junction voltage, 


( B-3 ) 


V A " V D " bwJ V 0 i V D * V oc 


(B-4) 


The desired operating point is obtained by maximizing the product JV A with 
respect to V D . This gives the following condition on V D> which is easily 
solved numerically: 

° - 0 ( B-5 ) 

The power density decreases as the cell width b increases for two reasons. 
First, for a fixed Rp, the cell current increases with b, causing the IR p 
drop to increase. Second, Rp itself increases with b. The first effect 
is demonstrated in Fig. B-4, where maximum available power density is 


i 

r av d 

av d 

1 

0 

•“3 

e D - 1 

1 + 2J^Ae • bwRp 





AV 


- av d J l e 



bwRp 


Figure B-4. Maximum available power density (relative) 
as a function of cell area times R_ 
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2 

plotted as a fraction of the ideal 12.9 mW/cm (Rp ■ 0). From this plot, it 
can be seen that bwRp must be restricted to 0.83 ohra-cra^ or less, to keep 
the relative efficiency from falling below 95%. More specific inferences 
may be drawn from this arbitrary criterion after the components of Rp have 
been quantitatively identified. 

COMPONENTS OF PARASITIC RESISTANCE 

Front Contact Resistance 


For convenience, it will be assumed that b/b' is an integer a 1 , so that the 

front area contains b/b" (see Fig. B-5) subunits per cell unit. From the 

symmetry of a subunit, assuming uniform J, it is obvious that the cell 

current bwJ is shared equally by 2b/b' half-contacts each with area equal 

B— 4 B— 5 

to d'w/2. According to the extended transmission line model (ETLM), * 
the series resistance due to one of these strips is approximately 

(h'/w) • R d ' • Yn' + 0.2 coth |(dV2h')/Vn > + 0.2^ (B-6) 

"• - o c '/[<i>-) 2 R 0 '] 



Figure B-5. Front subcell geometry for the case 
b/b" * 2 (distorted scale). 


B-4 

H. H. Berger, "Models for Contacts to Planar Devices," Solid-State 
Electronics 15 , 145 (1972). 

B— 5 

H. H. Berger, "Contact Resistance and Contact Resistivity," J. Electro- 
chem. Soc . 119 , 507 (April, 1972). 

t 
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where R ' is the sheet resistance (ohms/square) of the n + layer and p " is 
U 2 ^ 

the specific contact resistivity (ohms-cm ) of the contact interface. Since 

the 2b/b" half-contacts are in electrical parallel, the net front contact 

resistance for the cell unit is 

R ' - (b"/2b)(h'/w) • Rq' • Vn' 0.2 • coth { ( d"/2h ' Wn ' + 0.2^ (B 7) 
Front Silicon Layer Resistance 

The subunit of area wb' is modeled here as a sheet with current Jwb'/2 
entering at each edge and exiting with uniform density, J, through the 
surface (Junction). 

If positive current is from left to right and x is the coordinate measured 
from edge to edge (0 £ x £ b"), the sheet current in the present approxi- 
mation is 

I(x)»wJ(b'/2-x) (B-8) 


and the voltage relative to the edges, not including the voltage drop in 
the contacts, is 

V(x) » - ± R Q j* I(z)dz » - JR Q '(b'x - *' 2 )/2 (B-9) 

The minus sign is required because positive current flows away from the 
edges. The power dissipated in the subunit is 

p' =»/ b wJV(x) dx - J 2 R Q 'wb' 3 /12 (B-10) 

-'o I ! 


The front layer resistance is power (for the whole bxw unit) divided by 
the square of the current : 


(b/b ' )P'/ ( wbJ )' 


14 


b b' 
w E 


(B-ll ) 


Back Layer Resistance 


Because of the direction of current from the junction toward the back con- 
tacts, the voltage in the back layer increases toward the center of the unit. 
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By a one-dimensional model analogous to the one in the preceding subsection, 
the back surface voltage relative to the silicon side of the back contacts 
is 

V(x) - JR 0 (bx - x 2 )/2 (B-12) 

where R Q is the sheet resistance of the base layer. 

Dissipated power is 

P >■ J 2 R q wb 3 /12 (B-13 ) 

and back layer resistance is 

h -(ts)(v) r c < b - 14 > 

Back Contact Resistance 


The analysis of front contact resistance applies also to the back contacts. 
There are in this case, however, only two half-contacts per unit cell, each 
with area equal to dw/2. In terms of the base layer thickness, h, sheet 
resistance, R , and specific contact resistivity, p c> the back contact re- 
sistance is 

R c “ (i )(§) R o^n ♦ 0.2 coth {(d/2h)/Vn + 0 . 2 } (B-15 ) 

where n - p Q / (h 2 R Q ). 

ESTIMATION OF CONTACT SPACING 

According to Fig. B-4, the area resistance product, bwRp, should be re- 
stricted to the order of 1 ohm-cm 2 . Using Eqs. (B-7), (B-ll), (B-14), and 
(B-15), one can calculate bwRp as a function of b and then estimate maximum 
allowable b. This will depend on base layer resistivity through both base 
sheet resistance and back contact resistance, if aluminum contacts are used 

on p-type silicon. In this case, p depends upon the silicon bulk resist!- 
B— 5 ^ 

vity. Three different base-layer resistivities are considered: p D ■ 0.3, 

O 

1.0, and 3.0 ohras-cm. The n+ layer is assumed to be 0.3pm thick with a sheet 

resistance of 50 ohms/square. The front contact length, d', is 0.3mm, and 
- 5 2 

P c = 10 ohm-cm . The base-layer thickness is 125pm, and the back contact 

length, d, is 0.8rara in accordance with the substrate slot width. Back con- 
B— 5 “4 2 

tact resistivity is 5 x 10 pg ohm-cm , where p Q is the bulk resistivity 
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of the p-type silicon. From Eqs. (B-7), (B-ll), (B-14), and (B-15), we 
calculate 

R c ' - 0.01122/w 
R^* ■ 4.17 b/w 

R l - 0.0833 R q b/w (U-16) 

R c • 0.00306 R 0 /w 

P D 

where R Q « - 24, 80, and 240 ohms/square. 

The total bwRp* bw (R^ ♦ R L " + R L + R c ) is plotted in Fig. B-6. For bwR p 
not to exceed 1.0 ohm-cm 2 , the slot separation can be as much as 4ram for 
^ □ = ^4 , but should be held to 2 or 3mm for more lightly doped base material. 



0 12 3 4 

b (mm) 

Figure B-6. bwR a product vs. b for solar cells made 
on slotted ceramic substrates. 
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CIRCULAR CONTACTS 


A possible alternative to striped electrodes is a two-dimensional square 
lattice of small circular electrodes as shown in Fig. B-7. Various methods 
might be used to produce such a pattern. One method of making small, 
roughly circular holes in the ceramic, while at the same time creating 
troughs in the back ceramic surface through which contacts to the holes 
can be made, is illustrated in Fig. B-8. 


• • - 


O 

O 

o 

o 


o 

o 


o ! o ! o 

• ' 

o.oo 

I I 

-* t--- 

o ! o ! o 

. i 

— f — t- — -» 


o 

o 


° i ° ; 

o i o ; 

j. i 



Figure B-7. Partial lattice of electrodes with square basic 
units (left). In the derivation of Eq , (B-17), 
the square units are approximated in circles 
( right ) . 





Figure B-8. Holes at the intersection of perpendicular 
sets of troughs in the front and back 
ceramic surfaces. 
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Neglecting variation in the thickness of the silicon layer that might result 

from such a geometry, we have approximated the square fields of Fig. B-7 

as roughly equivalent circular fields in order to estimate the series 

resistance. If d and d, are the diameters of contacts and circular fields, 
o l 

respectively, it can be shown that the HA product for an effective sheet 


resistance, R n , is: 
HA 


R a d^/8 [in(d 1 /d 0 ) - 1/4 1 - d 1 2 /d 0 2 ) |l - Sd^/d^ 

L ' 


(B-17) 


and this equation is plotted in Fig. B-9. One may note, for example, that 
RA - 1.5 ohras-cm 2 if d Q - 0.3mm, d t - 3mm, and R n - 50 ohms/square. 



Figure B-9. RA/Rq from Eq.(B-16) for the circular 
base field geometry. 

The contribution to RA due to contact resistance can be shown to be compara- 

B-4 

tively small. For this term, we use the transmission line model for 
circular geometries: 
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contact 


(B-18) 


j 


(RA) 



~o] 





whore J and J are the Bessel functions of orders 0 and l, respectively, 
j -yrT. o is the specific constant resistivity, and a - M a /p c' For R 0 
- 50 ohms/square and p c - 0.0005 ohm-cm~ , one finds ( R A) oontact 3 ° 2 ohm- cm 

The RA product, including contact resistance, is therefore approximately 
1.7 ohms-cm 2 . comparing favorably with the striped geometry of equivalent 
(0.3cm) spacing. For both types of geometry, the RA product will fall 
rapidly as the spacing is reduced. 
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APPENDIX C 

NOVEL CELL DEVELOPMENT 


Because silicon -on-ceramic (SOC) material is on an insulating substrate, it 
is difficult to fabricate an electrical contact to the base region of the 
cell. At the start of this work, the base contact was made by using an 
interdigital electrode pattern on the top surface. This method was ob- 
viously undesirable since the pattern covered too much area. Later in the 
program, a novel high-voltage structure (suggested by Dr. G. Turner of 
JPL ) was Investigated as an alternate technique for achieving contact to 
the base layer as well as providing a higher terminal voltage. 

The high-voltage structure is fabricated as follows: Carbon stripes are 
deposited on a ceramic substrate. The substrate is then dip-coated with 
silicon so that the carbon stripes are perpendicular to the surface of the 
melt. The resulting silicon stripes, after diffusing a p-n junction, can 
be series-connected to form a high-voltage structure. The active area of 
this structure is only one-half of the total area. To offset this area loss, 
an integral optical coupler (IOC) is added as illustrated in Fig. C-l. 

LOAD 



Figure C-l. Combined integral optical coupler 
high-voltage cell device. 
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This IOC is a cylindrical, trough-like, dielectric-compound , parabolic con- M 

centrator having a concentration ratio of 3. As compared with the inter- 
digital electrode cell, the high-voltage structure allows efficient col- 
lection of incident light even though all metal contacts are located on the ** 

top surface of the structure. The incident light, of course, is concen- 
trated primarily on the active area. Because the short-circuit current, J , 
of the cell increased by a factor of 2 when the optical coupler was added, 
the cell was designated the 2X cell." 

2X cells fabricated from striped SOC material had conversion efficiencies 

from 4 to 6% (AMI, no AR coating). The major problem was an undesirably 

high series resistance. mk 

In parallel with this novel cell study, a slotted-subslrate method for con- 
tacting the back surface of the cell was developed. This method is dis- 
cussed in detail in the Technical Discussion section of this report. 

In this appendix, the series resistance of the 2X structure and tha slotted- 
substrate structure will be calculated. Cross-sectional views of tnese de- 
vices are shown in Fig. C-2. 

For the slot ted-substrate cell, the parasitic series resistance, can 

be divided into three components: the front contact resistance, R q the 
silicon-layer resistance, ' ; and the back-layer resistance, R, . The back 
contact resistance can be neglected. The front-contact resistance, R ', is 
assumed to be given by 

V ’ V /l (C_1) 


where l is the peripheral length of the emitter metallization pattern and 
R,' is the specific contact resistance in ohras-cra. An accurate value of 
R ^ is not known; we assume here that R^ “ - 0.3 ohm/cm. Thus, R * - 0.013 
ohm, since i «* 24cm. The front silicon-layer resistance, R^', is given by 



h R n 



where is the front-layer sheet resistance in ohms/square, b' is the unit 
cell width, k' is the width of a contact finger, N' is the numbsr of unit 
cells per solar cell, and w is the width of the solar cell For a typical 
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200um 



(b) 2X SOC cell. 

Figure C-2. Cross-sectional drawings of a slotted-substrate 
SOC cell and a 2X-S0C cell. 

SCC ,cell, Rg - 50 ohms/square, b' - 0.34cm, k/ ■ 0.03cra, N' - 7.06, and 
w - 1.7cm, given that R L ' « 0.108 ohm. For the back-lay r resistance, R^, 
we use the following equation: 


R o 

TiTBJw 


b 2 

Y1 



(C-3) 


where R Q is the sheet resistance of the back layer, b is the center-to-center 
distance between 3lots, d is the slot width, h is the thickness of the 
back layer, N is the number of slots per solar cell, and w is ‘„he width of 
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!>. solar cell. Thus, (Nb>. Is the total are. of the solar cell. By making 

t:/™" 8 ' ' ,S ““ reWr ‘ te EQ ' <C " 3) t,r " 5 ° f base concentra- 

“Hh e B ' 1 »*T th SPeed ’ V Flr ' t ’ " — " that b “ 8 « -«* resistivity, 
quals R 0 h, is inversely proportional to N Q . Thus, 


R Q h - K/N q 


( C-4 ) 


where K is an empirical constant. For SOC material with a base donin* 
concentration^ N Q ^ of 2.3 x 10 16 /cm 3 , ft 8 
K 


23a lo 10 /™ 0 ° * 1U /Cm ' fl a h 13 about 1 ohm- cm, given that 

2-3 a 10 /cm . Second. we assume that the silicon thickness. h, Is 


icon cnicKness, h 

nversel, proportional to the square of the growth velocity, v. Thus, 


h - K*/v 2 


(C-5 ) 


where K* is an empirical constant, 
aoout 200ym. Thus, K* is about 7.2 
b » 0. 25cra, d - 0.1cm, and (Nb)w - 
and (C-5) combine to give 


For a growth speed of 0.06 cra/sec, h is 
x 10 cm . For a typical SOC cell, 
4.08cm , so that Eqs . (C-3). (C-4), 


1.49 x 10 13 


tt: 


B 


v 2 + 6 ■ 912 x 10 


- 6 ! 


(C-6 ) 


Because the totaJ parasitic 
specific series resistance, 


resistance, Rp, equals R' + R' + R 
RpA, is given by 


the 


RpA => 0.494 + 


6.07 x 10 ] 


N. 


B 


„2 A 6.912 x 10- 

2 

v 


( C-7 ) 


since A = 4.08cm' 


Equation (C-7) yields the data in Table C-l, 


of R- aid R' WhiCh 18 illUStrated in Fi *- -3. The values 

slot? r > *" aSSU " ed t0 be dOUble the cor tesponding values for the 
slotted-substrate SOC cell. The back-layer resistance, R ls given by 


1 

lS 


u iTw 


+ i R ^ 

2 a Nw 


(C-8) 
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where R^ is the sheet resistance of the back layer, l is the width of the 
active region of the unit cell, N is the number of unit cells per solar 
coll, w is the width of the solar cell, and s is the lateral distance from 
the edge of the active region to the edge of the base contact. By Fig. C-l, 
R ■ 50 ohms/square , t - 0.31cm, N * 4, w * 1.7cm, and s - 0.09cm. Thus, 

Rj - 0.52 ohm, R p - 0.76 ohm, and the specific series resistance of the 2X 
cell is 3.2 oluas-cm“ . This value is too high for high-efficiency cells. 


Table C-l. Specific series resistance vs. base doping 
concentration and thickness 


V 

( cm/sec ) 

h 

C um ) 

2 

Specific Series Resistance (rt-cra ) 

Base Doping 

- 3 ) 

Concentration (cm 

2.3 x 10 16 

10 16 

4.7 x 10 15 

0.06 

200 

0.64 

0.83 

1.21 

0.07 

147 

0.66 

0.88 

1.31 

0.08 

113 

0.69 

0.95 

1.46 


INTEGRAL OPTICAL COUPLER 



Figure C-3. The 2X-eell. 
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